Modelo Constitutivo de Mohr

Coulomb

Integrantes:
Esteban Martín Caballero Campos
Miguel Enrique Paats Nicora

Índice
1.Conceptos Básicos .............................................................................................. 2
a)Idealización del material................................................................................. 2
b)Tensión, Tensor e invariantes ........................................................................... 2
c)El criterio de Fluencia y el potencial Plástico ................................................ 5
d)Relación elastoplástica tensión deformación .............................................. 6
2.Modelo Constitutivo de Mohr Coulomb ............................................................ 7
3.Bibliografía: ........................................................................................................... 9
 1. Conceptos Básicos

Los modelos constitutivos se pueden definir como procedimientos numéricos

que intentan describir y representar de manera aproximada el comportamiento el
comportamiento de los materiales. Algunos conceptos fundamentales de la teoría
de la mecánica de medios continuos (MMC) son presentadas para facilitar la
compresión de los conceptos posteriores.

a) Idealización del material

El modelo constitutivo de Mohr Coulomb idealiza al material como un material
isótropo, homogéneo y continuo. Con un comportamiento elastoplástico perfecto.
Como consecuencia la relación tensión deformación será lineal con una pendiente
E, que representa la región de las deformaciones elásticas (ξe) del material, hasta
que se alcance un valor de tensión fluencia iniciando la región de las deformaciones
plásticas (ξp), volviendo a ser una relación lineal pero esta vez con pendiente 0, ya
que esta tensión de fluencia no puede ser superada (no existe endurecimiento por
las deformaciones plásticas) y las deformaciones plásticas son irreversibles, ver figura
1. Por esto, el modelo de Mohr Coulomb es considerado como una aproximación de
primer orden del comportamiento real del geo material porque el comportamiento
tensión deformación es curvo, línea discontinua en la figura 1.

Figura 1 - Grafico Tensión-deformación. En línea llena el comportamiento elastoplástico

perfecto idealizado, en línea discontinua el comportamiento real observado en laboratorio.

b) Tensión, Tensor e invariantes

En matemáticas y en física, un tensor es cierta clase de entidad algebraica de
varios componentes que generaliza los conceptos de escalar, vector y matriz de una
manera que sea independiente de cualquier sistema de coordenadas elegido. El
tensor de tensiones de Cauchy consiste en nueve componentes σij que definen
completamente el estado de tensiones en un punto dado dentro de un material. El
tensor T relaciona la dirección de un vector unidireccional ei con las tensiones Tei
que se desarrollan en el punto considerado sobre el plano perpendicular al vector
ei. Gráficamente se ilustra en la figura 2.
 σ11 σ12 σ13
𝑇 = [σ21 σ22 σ23 ]
σ31 σ32 σ33
σ11 σ21 σ31 1 σ11
𝑇 = 𝑇 . 𝑒1 = [σ12
𝑒1 𝑇 σ22 σ32 ] . [0] = [σ12 ] Tensión sobre el plano perpendicular a X .
1

σ13 σ23 σ33 0 σ13

σ11 σ21 σ31 0 σ21
𝑇 = 𝑇 . 𝑒2 = [σ12
𝑒2 𝑇 σ22 σ32 ] . [1] = [σ22 ] Tensión sobre el plano perpendicular a X2.
σ13 σ23 σ33 0 σ23
σ11 σ21 σ31 0 σ21
𝑇 𝑒3 = 𝑇 𝑇 . 𝑒3 = [σ12 σ22 σ32 ] . [0] = [σ22 ] Tensión sobre el plano perpendicular a X3.
σ13 σ23 σ33 1 σ23

Figura 2 - Estado tensional de un punto en 3D. Fuente: Wikipedia

De acuerdo con la ley de conservación del momento angular, el equilibrio

requiere que la suma de los momentos con respecto a un punto arbitrario sea cero,
lo que conlleva a la conclusión de que el tensor de tensiones es una matriz simétrica,
teniendo entonces solo 6 componentes independientes y no nueve.

σ𝑖𝑗 = σ𝑗𝑖 ∀𝑖 ≠ 𝑗 ∴ 𝑇 𝑇 = 𝑇

Si bien las componentes σij de un tensor de tensiones T depende de la

orientación de los ejes de coordenadas optados para definir la matriz de tensiones T
en el punto de consideración. La matriz de tensiones T en si es una entidad
independiente del sistema de coordenadas optado.

Hay ciertos valores asociadas con cada tensor que son constantes,
independientemente al sistema de coordenadas utilizado para representar el tensor,
a estas constantes se le denominan invariantes del tensor.

Por ejemplo, un vector es un tensor de rango uno. En tres dimensiones posee

tres componentes. El valor de estas tres componentes depende del sistema de
coordenadas utilizado para describir el vector, pero su magnitud es una cantidad
física (un escalar) y es independiente del sistema de Coordenadas cartesianas
escogido.

Similarmente, la matriz de tensiones T tiene tres variables independientes

asociadas con él. Un conjunto de dichas invariantes son los esfuerzos principales, que
son los eigenvalues λi de la matriz de tensiones. Los eigenvectors vi contienen la
dirección de los esfuerzos principales expresadas en el sistema de coordenadas
electo para escribir la matriz de tensiones T. Para determinar se los eigenvalues λi y
eigenvectors vi se debe resolver la siguiente ecuación:

(𝑇 − 𝜆𝑖 𝐼). 𝑣𝑖 = 0 donde I es la matriz de identidad

Pudiendo también resolverse el polinomio característico de la matriz T para
determinar los eigenvalues λi :
3 2 1 0
det(𝑇 − 𝜆𝑖 𝐼) = 0 ∴ −𝜆 + 𝐼1 𝜆 − 𝐼2 𝜆 + 𝐼3 𝜆
𝐼1 = 𝜎11 + 𝜎22 + 𝜎33 Primer Invariante del tensor de tensiones T.
𝜎11 𝜎12 𝜎11 𝜎13 𝜎22 𝜎23
𝐼2 = 𝑑𝑒𝑡 ([𝜎 𝜎22 ]) + 𝑑𝑒𝑡 ([ 𝜎31 𝜎33 ]) + 𝑑𝑒𝑡 ([𝜎32 𝜎33 ]) Segundo Invariante del tensor
21
de tensiones T.

𝐼3 = det (𝑇) Tercer Invariante del tensor de tensiones T.

La notación usual de las tensiones principales es 𝜎1 > 𝜎2 > 𝜎3 y es la que se

respetara a lo largo de este trabajo.

Todo tensor de tensiones puede descomponerse en dos matrices de tensiones,

una llamada tensiones desviadoras y otras denominada tensiones hidrostáticas.
𝑠11 𝜎12 𝜎13 𝑝 0 0
𝜎
𝑇 = 𝑆 + 𝑃 = [ 21 𝑠22 𝜎23 ] + [0 𝑝 0]
𝜎31 𝜎32 𝑠33 0 0 𝑝
1 1
Donde 𝑝 = (𝜎
3 11
+ 𝜎22 + 𝜎33 ) = 3 𝐼1

Al igual que la matriz de tensiones T las tensiones desviadoras S también

poseen sus invariantes J1, J2, J3 .

Para terminar, cabe resaltar que al momento de resolver un problema existen

2 sistemas de coordenadas de interés. El primero, es el sistema de coordenadas
globales por lo general cuenta con el eje vertical paralelo a la dirección de la
gravedad y dos ejes horizontales que se definen arbitrariamente para acomodar
mejor la geometría de la estructura estudiada. El segundo es un sistema local de
coordenados que tiene como eje de coordenadas, ejes paralelos a los vectores de
esfuerzo principal Vi en el punto estudiado. Es en este último sistema donde se definen
los modelos constitutivos.

Además, cabe resaltar que existe una dirección particular de interés que es el
eje de las tensiones hidrostáticas definido por la dirección que posee igual
espaciamiento angular con los vectores de las tensiones principales. Un plano que
contenga el eje hidrostático es conocido como plano meridional y un plano
perpendicular a dicho eje es conocido como plano de desviación.

Figura 3-Eje de Coordenadas locales referido a las tensiones principales y eje hidrostático.
Fuente: Continuum Mechanics Course (MMC) - ETSECCPB – UPC CH 08

c) El criterio de Fluencia y el potencial Plástico

El criterio de fluencia, la potencia plástica y la ley de endurecimiento
usualmente definen el modelo constitutivo.

Como se ha visto en anteriormente el termino elastoplástico se refiere a que el

mecanismo de estado de deformaciones está compuesto puramente por una
componente elástica y otra plástica. Para evaluar si una deformación plástica ocurre
o no en un punto del material, definimos el criterio de fluencia como una función de
las tensiones principales σ en el punto estudiado, y un parámetro de estado del
material α.

𝑓(𝜎, ∝) < 0 ∴ 𝐷𝑜𝑚𝑖𝑛𝑖𝑜 𝑑𝑒 𝑙𝑎𝑠 𝑑𝑒𝑓𝑜𝑟𝑚𝑎𝑐𝑖𝑜𝑛𝑒𝑠 𝐸𝑙𝑎𝑠𝑡𝑖𝑐𝑎𝑠

𝑓(𝜎, ∝) = 0 ∴ 𝐷𝑜𝑚𝑖𝑛𝑖𝑜 𝑑𝑒 𝑙𝑎𝑠 𝑑𝑒𝑓𝑜𝑟𝑚𝑎𝑐𝑖𝑜𝑛𝑒𝑠 𝑃𝑙𝑎𝑠𝑡𝑖𝑐𝑎𝑠
𝑓(𝜎, ∝) > 0 ∴ 𝐸𝑠𝑡𝑎𝑑𝑜 𝑖𝑛𝑎𝑑𝑚𝑖𝑠𝑖𝑏𝑙𝑒
𝜎1
𝜎
𝜎 = [ 2 ] 𝑣𝑒𝑐𝑡𝑜𝑟 𝑑𝑒 𝑙𝑎𝑠 𝑡𝑒𝑛𝑠𝑖𝑜𝑛𝑒𝑠 𝑝𝑟𝑖𝑛𝑐𝑖𝑝𝑎𝑙𝑒𝑠 𝑦 ∝ 𝑝𝑎𝑟𝑎𝑚𝑒𝑡𝑟𝑜 𝑑𝑒 𝑒𝑠𝑡𝑎𝑑𝑜 𝑑𝑒𝑙 𝑚𝑎𝑡𝑒𝑟𝑖𝑎𝑙
𝜎3

El criterio de falla puede ser visto como una hiper superficie en el espacio de
las tensiones principales, representando la superficie de falla cuando f = 0, y es
conveniente formulado a través de las invariantes del tensor de tensiones T. Cabe
destacar que la superficie de falla en si es una invariante ya que el material debe
fallar de igual manera independientemente del sistema de coordenadas optado.
 En análisis no lineales, un acercamiento incremental es usado para acomodar
las no linealidades. Cuando el incremento de la deformación resulta en la fluencia
del material f = 0, una deformación plástica ocurre y ya no existe una relación única
entre las tensiones y el incremento de las deformaciones. Los incrementos de
deformaciones son descompuestos en deformaciones puramente elásticas y
deformaciones puramente plásticas.
𝑑𝜀 = 𝑑𝜀 𝑒 + 𝑑𝜀 𝑝

Las deformaciones elásticas están gobernadas por la ley de Hooke, y las

deformaciones plásticas son descriptas por la regla de flujo.
𝜕𝑔(𝜎, 𝛽)
𝑑𝜀 𝑝 = 𝑑𝜆
𝜕𝜎
𝜎1
𝜎 = [𝜎2 ] 𝑣𝑒𝑐𝑡𝑜𝑟 𝑑𝑒 𝑙𝑎𝑠 𝑡𝑒𝑛𝑠𝑖𝑜𝑛𝑒𝑠 𝑝𝑟𝑖𝑛𝑐𝑖𝑝𝑎𝑙𝑒𝑠 𝑦 𝛽 𝑝𝑎𝑟𝑎𝑚𝑒𝑡𝑟𝑜 𝑑𝑒 𝑒𝑠𝑡𝑎𝑑𝑜 𝑑𝑒𝑙 𝑚𝑎𝑡𝑒𝑟𝑖𝑎𝑙 (𝑐𝑡𝑒)
𝜎3

𝑆𝑖 𝑔 = 𝑓 → 𝐷𝑒𝑓𝑜𝑟𝑚𝑎𝑐𝑖𝑜𝑛𝑒𝑠 𝑃𝑙𝑎𝑠𝑡𝑖𝑐𝑎𝑠 𝑜 𝑟𝑒𝑔𝑙𝑎 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜 𝐴𝑠𝑜𝑐𝑖𝑎𝑑𝑎𝑠

𝑆𝑖 𝑔 ≠ 𝑓 → 𝐷𝑒𝑓𝑜𝑟𝑚𝑎𝑐𝑖𝑜𝑛𝑒𝑠 𝑃𝑙𝑎𝑠𝑡𝑖𝑐𝑎𝑠 𝑜 𝑟𝑒𝑔𝑙𝑎 𝑑𝑒 𝑓𝑙𝑢𝑗𝑜 𝑛𝑜 𝐴𝑠𝑜𝑐𝑖𝑎𝑑𝑎𝑠

Donde dλ es llamada el multiplicador plástico, y, g es una función escalar

denominada potencial plástico. Como las deformaciones plásticas son irreversibles
el multiplicador plástico está restringido a valores mayores a cero d λ >0. El
multiplicador plástico dλ controla la magnitud de la deformación plásticas y es
dependiente de la ley de endurecimiento, para este caso es un valor constante.

En geotecnia el cambio de volumen durante la deformación plástica es un

fenómeno usual. Este fenómeno es gobernado a través de la función del potencial
g.
plástico En materiales granulares la dilatación del material es descripto por el
ángulo Ψ.

d) Relación elastoplástica tensión deformación

Para un material homogéneo e isótropo la relación elástica lineal esta dada
por la siguiente ecuación:

𝜎 = 𝐷𝜀 Ley de Hooke

𝜎1 1−𝑣 𝑣 𝑣 𝜀1
𝐸
⌈𝜎2 ⌉ = (1+𝑣)(1−2𝑣)
[ 𝑣 1−𝑣 𝑣 ] [𝜀2 ] E: Modulo de Young y v: módulo de Poisson
𝜎3 𝑣 𝑣 1 − 𝑣 𝜀3

Cuando la fluencia ocurre en punto del material, la ecuación de Hooke ya no

es aplicable, y se necesita una nueva relación entre el incremento de tensión y
deformación. Una relación para el criterio de fluencia con un incremento infinitesimal
en las tensiones puede ser dada como:
 𝜕
𝑓(𝜎 + 𝑑𝜎) = 𝑓(𝜎) + (∇𝑓)𝑇 𝑑𝜎 = 0 𝑑𝑜𝑛𝑑𝑒 ∇ =
𝜕𝜎
Como 𝑓(𝜎) = 0 por encontrarse en el régimen plástico, entonces debe
cumplirse que el segundo término también sea igual cero.
(∇𝑓)𝑇 𝑑𝜎 = 0 𝐶𝑜𝑛𝑑𝑖𝑐𝑖ó𝑛 𝑑𝑒 𝐶𝑜𝑛𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎

Un incremento infinitesimal de las tensiones puede ser expresado por la ley de

Hooke y la regla de flujo como:
𝜕
𝑑𝜎 = 𝐷𝑑𝜀 𝑒 = 𝐷(𝑑𝜀 − 𝑑𝜀 𝑝 ) = 𝐷(𝑑𝜀 − 𝑑𝜆∇𝑔) 𝑑𝑜𝑛𝑑𝑒 ∇ =
𝜕𝜎
Substituyendo esta nueva ecuación en la ecuación de consistencia se
obtiene:
(∇𝑓)𝑇 𝐷 𝑑𝜀
(∇𝑓)𝑇 (𝐷(𝑑𝜀 − 𝑑𝜆∇𝑔)) = 0 ↔ 𝑑𝜆 = 𝑝𝑜𝑟 𝑠𝑒𝑟 𝑑𝜆 = 𝑐𝑡𝑒
(∇𝑓)𝑇 𝐷 ∇𝑔

Colocando el valor dλ en la ecuación del incremento de tensiones

infinitesimal se obtiene la matriz de rigidez elastoplástica:
𝐷 ∇𝑔 (∇𝑓)𝑇 𝐷
𝑑𝜎 = 𝐷 𝑒𝑝 𝑑𝜀 = (𝐷 − ) 𝑑𝜀
(∇𝑓)𝑇 𝐷 ∇𝑔

En la practica la matriz Dep no suele ser consistente y conlleva una pobre

convergencia, por lo que es reemplazada por la matriz Depc, matriz de rigidez
elastoplástica consistente.

𝐷 𝑐 ∇𝑔 (∇𝑓)𝑇 𝐷 𝑐
𝐷 𝑒𝑝𝑐 = (𝐷 𝑐 − ) 𝑑𝑜𝑛𝑑𝑒 𝐷 𝑐 = (𝐼 + 𝑑𝜆𝐷∇2 𝑔)−1 𝐷
(∇𝑓)𝑇 𝐷𝑐 ∇𝑔

2. Modelo Constitutivo de Mohr Coulomb

En esta sección se definirán el criterio de fluencia y la función de potencial

plástico utilizado por el modelo constitutivo de Mohr Coulomb, así también con
algunas ventajas y desventajas del modelo. El criterio data de 1773 y ha sido
descripto extensivamente en la literatura. Es todavía uno de los criterios más utilizados
para describir materiales que resisten por mecanismo de fricción. ´

Si los resultados de ensayos triaxiales o corte directo, en suelos, son graficados

en términos de esfuerzos efectivos, es usual obtener un grafico aproximado al grafico
mostrado en la figura 4. Esta línea es el criterio de fluencia f de Mohr Coulomb.

𝑓(𝜎, ∝) = 𝜎1, − 𝜎3, − (𝜎1, + 𝜎3, ). sin 𝜑, − 2𝑐 , . cos 𝜑, Criterio de fluencia Mohr Coulomb

Esta ecuación también suele ser reescrita como:

1 + sin 𝜑,
𝑓(𝜎, ∝) = 𝑘𝜎1, − 𝜎3, − 𝜎𝑐, 𝑑𝑜𝑛𝑑𝑒 𝑘 = ,
𝑦 𝜎𝑐, = 2𝑐 , √𝑘
1 − sin 𝜑
 Figura 4-Superficie de falla obtenida en ensayos triaxiales.

Donde σ1 yσ3 son las tensiones efectivas principal mayor y menor

respectivamente. C´ y ϕ´ son la cohesión y la fricción efectivas respectivamente del
material estudiado.
𝑐,
𝛼 = [𝜑, ] = 𝑐𝑡𝑒 Parámetro de estado del criterio de fluencia de Mohr Coulomb

En el espacio de tensiones efectivas la función de fluencia genera una

pirámide hexagonal irregular como se muestra en la Figura 5.

Figura 5-Superficie de Falla en el espacio de las tensiones efectivas

Figura 6-Proyeccion de la superficie de falla en un plano desviador (a) y en un plano

meridiano (b)
 Para completar el modelo constitutivo de Mohr Coulomb se debe definir la
función de potencial plástico g. Para el criterio de Mohr Coulomb se puede asumir
una función de potencial plástico asociada es decir g = f. La desventaja de este
enfoque es que la magnitud de las deformaciones volumétricas plásticas predichas
por el modelo son mayores que las observadas, la otra es que el modelo predice que
el suelo en estado fluencia se dilatara permanentemente, pero la experiencia
demuestra que el suelo alcanza un estado de deformación a volumen constante
para grandes deformaciones.

El primero de los problemas puede ser solucionado introduciendo el

′
parámetro de estado 𝛽 = [ 𝑎 , ] en g, siendo Ψ’ el ángulo de dilatación y a’ el valor
𝛹
de la cohesión en estado plástico pudiendo ser considerado a’ = 0 o a’ = C’.

𝑔(𝜎, 𝛽) = 𝜎,1 − 𝜎,3 − (𝜎,1 + 𝜎,3 ). sin 𝛹, − 2𝑎, . cos 𝛹, Potencial Plástico en Mohr Coulomb

Aunque con el uso del criterio del potencial plástico no asociado 𝑔 ≠ 𝑓 se logre
salvar las grandes deformaciones volumétricas, el modelo aun predice el incremento
de volumen con el incremento en la deformación, sin importar cuanto el suelo
deslice. Esto es una predicción no realista y puede derivar en problemas de contorno.
Además, como se puede apreciar tanto en el criterio de fluencia como en el
potencial plástico no se tiene en consideración del efecto del esfuerzo principal
intermedio σ2, y posee una superficie de falla discontinua lo que conlleva a ciertos
problemas de convergencia en los problemas.

Sin embargo, su ventaja radica en que brinda una buena aproximación a la

mayoría de los problemas de geotecnia y necesita la determinación de pocos
parámetro para poder ser implementado, tan solo cinco (C’ ϕ´ , E, u, Ψ’).

3. Bibliografía:
➢ Cauchy stress tensor - Wikipedia
➢ Cálculo tensorial - Wikipedia, la enciclopedia libre
➢ Lode coordinates - Wikipedia
➢ Elasto-Plastic Constitutive Modelling of Geotechnical Material
Investigation of the Intermediate Principal Stress Influence and
Exploration of the General Parametric Yield Surface Format - Johan
Nielsen & Kenneth Sixhøj Jepsen - Aalborg University Esbjerg
➢ DETERMINACIÓN DE PARÁMETROS PARA LOS MODELOS
ELASTOPLÁSTICOS MOHR-COULOMB Y HARDENING SOIL EN SUELOS
ARCILLOSOS - Andrés Nieto Leal, Javier Fernando Camacho Tauta,
Edwin Fernando Ruiz Blanco - Revista Ingenierías Universidad de
Medellín
➢ Finite element analysis in geotechnical engineering Theory - David M.
Potts and Lidija Zdravkovic, 1999
➢ Continuum Mechanics Course (MMC) - ETSECCPB – UPC - CH 08
“Plasticity”