MODELO MOHR-COULOMB

En las últimas décadas, un área importante del estudio de la mecánica de suelos ha sido
conducida al desarrollo de modelos que se ajustan mejor a la realidad física de estos
geomateriales. Así, por ejemplo, varias investigaciones formularon con éxito modelos
basados en diferentes conceptos como la elasticidad no lineal, el endurecimiento plástico,
la teoría del estado crítico y la hipoplasticidad. No obstante, el desarrollo de la
investigación constitutiva implicó, al mismo tiempo, una mayor dificultad tanto teórica
como práctica, traducida en la determinación de un número significativo de parámetros,
la aplicación de procedimientos matemáticos y ensayos experimentales complejos, y la
consideración de características y variables dependientes a la naturaleza de cada tipo de
suelo. Esto, sumado a información de laboratorio limitada, hace de la metodología para la
determinación de los parámetros usados en el modelo, así como de la evaluación de su
capacidad y aplicabilidad frente a una tarea específica, partes fundamentales en todo
análisis de modelación [4]. De esta manera, a continuación se presenta un trabajo de
investigación para la determinación de los parámetros correspondientes a dos de los
modelos constitutivos más relevantes formulados en el marco de la elastoplasticidad: el
modelo Mohr-Coulomb y modelo de endurecimiento Hardening Soil; dicha metodología se
valida mediante la comparación entre las observaciones experimentales y la
implementación de los modelos en simulaciones numéricas. Los resultados permiten
identificar las posibilidades y condiciones de aplicación de los modelos considerados en la
representación del comportamiento mecánico de suelos arcillosos, y en el análisis y diseño
de estructuras geotécnicas mediante el uso de esta metodología en modelaciones
numéricas.

MODELO MOHR-COULOMB
El modelo constitutivo Mohr-Coulomb es considerado como una aproximación de primer
orden al comportamiento no lineal del suelo. Se trata de un modelo elastoplástico
perfecto (isotrópico) desarrollado a partir de la composición de la ley de Hooke y la forma
generaliza del criterio de falla Mohr-Coulomb. Puede simular el comportamiento de suelos
granulares sueltos o finos normalmente consolidados y se debe tener en cuenta que no
representa el comportamiento elastoplástico progresivo sino que es un modelo elástico y
luego plástico perfecto. Su formulación involucra dos elementos generales: la elasticidad
perfecta y la plasticidad asociada al desarrollo de deformaciones plásticas o irreversibles.
Para evaluar si la plasticidad ocurre o no en un análisis específico, el modelo Mohr-
Coulomb involucra un conjunto de funciones de ( f) que definen el límite entre el
comportamiento elástico y plástico del material. La representación gráfica de la funciones
(cono hexagonal) en el espacio de los esfuerzos principales constituye un contorno o
superficie de fluencia fija, como se muestra en la Figura 2a. De esta manera, para diversos
estados de esfuerzos representados dentro de la superficie, el comportamiento es
puramente elástico, y todas las deformaciones son reversibles. Cuando los esfuerzos
igualan o superan la frontera definida por esta superficie, se presentan deformaciones
tanto elásticas como plásticas. Este criterio de fluencia es una extensión de la ley de
fricción de Coulomb para un estado general de esfuerzos, y es definido a partir de seis
funciones formuladas en términos de esfuerzos principales.

El comportamiento elastoplástico perfecto en el modelo es formulado a partir de la

relación entre los dos tipos de deformaciones (tasas de deformación elásticas y plásticas)
mediante la aplicación de la ley de Hooke en su forma clásica.

Donde De es la matriz de rigidez (elástica) del material. No obstante, la teoría de

plasticidad asociada relativa a esta relación básica tiende a una sobreestimación del
fenómeno de dilatancia en el suelo, para lo cual el modelo Mohr-Coulomb en adición a la
función de fluencia f, incorpora una función de potencial plástico g, que en el caso
particular g ≠ f describe plasticidad no asociada. De esta manera, y en adición a las
funciones de fluencia, un grupo de funciones de potencial plástico g son definidas para el
modelo:
A partir de esta consideración, las deformaciones plásticas son expresadas como:

Donde el multiplicador plástico λ define el tipo de comportamiento del material (λ = 0,

para comportamiento netamente elástico; y λ > 0 en el caso de comportamiento plástico).
Las ecuaciones anteriores componen la relación fundamental del modelo Mohr-Coulomb
entre esfuerzos efectivos y deformaciones para elastoplasticidad.

De esta forma, la formulación del modelo exige cinco parámetros básicos de entrada: el
módulo de Young E, la relación de Poisson υ, la resistencia por cohesión del suelo c, el
ángulo de fricción interna ϕ, y el ángulo de dilatancia ψ, (tabla 1). Dado que los ingenieros
están por lo general familiarizados con los cinco parámetros anteriores y ocasionalmente
disponen de información adicional acerca de otros parámetros del suelo, el modelo Mohr
Coulomb es el modelo básico más utilizado en la práctica geotécnica, a pesar de su
incapacidad para reproducir adecuadamente los cambios de rigidez del suelo y modelar
situaciones donde diferentes trayectorias de esfuerzos son experimentadas.
 OBTENCIÓN DE PARÁMETROS DEL MODELO MOHR-COULOMB
En el apartado 2 se trataron las principales características, variables y formulación básica
del modelo elastoplástico perfecto Mohr-Coulomb. A continuación se describe en detalle
la metodología desarrollada para la determinación de los parámetros del modelo.
Parámetros de resistencia
1 - Ángulo de fricción ϕ

El ángulo de fricción es uno de los parámetros más comunes en la mecánica de suelos y es

utilizado para numerosos procedimientos de diseño de fundaciones y otras estructuras
geotécnicas. En general, el ángulo ϕ hace referencia a la medida de la resistencia cortante
debida a la fricción mecánica directa entre las partículas de suelo [13]. El valor del ángulo
de fricción caracteriza notoriamente a cada tipo de suelo, y puede ser determinado a
partir de procedimientos de laboratorio y de campo, a través de correlaciones y gráficos; y
si no se cuenta con ninguna otra información puede obtenerse a partir de tablas. En este
sentido, la prueba triaxial es uno de los ensayos más adecuados para evaluarlo mediante
la representación gráfica del criterio de falla Morh-Coulomb. La determinación del ángulo
de fricción total y efectivo para el suelo en estudio se muestra en la figura 4 a partir de los
resultados de los ensayos triaxiales (CU) efectuados.
 2- Cohesión
La resistencia por cohesión es definida como la medida de las fuerzas que cementan las
partículas de suelo [13]. La cohesión del suelo puede ser establecida a partir de variados
procedimientos; sin embargo, a partir de la información disponible de los ensayos
triaxiales ejecutados y del trazado de la envolvente de falla del suelo, esta puede ser
determinada como el intercepto de la línea de falla con el plano vertical. El valor de la
cohesión es cero, un valor típico para arcillas normalmente consolidadas. No obstante, en
los cálculos computacionales es aconsejable valores ligeramente mayores para un mejor
desempeño en el cálculo numérico.
C. Ángulo de dilatancia ψ El ángulo de dilatancia puede tener una influencia significativa
en la modelación del comportamiento del suelo, especialmente en arenas densas y arcillas
altamente sobreconsolidadas. En el caso de las arenas este parámetro puede ser
determinado mediante la evaluación de las deformaciones unitarias volumétricas en un
ensayo triaxial drenado, considerando de manera particular la variación del parámetro
con el grado de sobreconsolidación (OCR) del suelo. Dado que de manera general la
dilatancia corresponde a un fenómeno particular de suelos granulares, y por tratarse de
especímenes de suelo con comportamiento normalmente consolidado es recomendable
usar el valor de cero grados para este parámetro.
3- Parámetros elásticos
A. Módulo de Young E
Plaxis usa el módulo de Young como un módulo de rigidez, y dado el comportamiento no
lineal típico de los suelos normalmente consolidados, no es conveniente usar el módulo
de rigidez tangente a la parte inicial de la curva esfuerzo deformación (zona elástica) sino
que se recomienda usar el módulo de rigidez correspondiente a la pendiente de la recta
entre el origen y el punto correspondiente al 50% de la resistencia máxima como se
presenta en la figura 5. Este módulo de Young es denotado como E50 y es el más indicado
para representar dicha respuesta. El valor de este parámetro se determinó gráficamente a
partir de la información experimental obtenida de las pruebas triaxiales.
Dado que el modelo Mohr-Coulomb no considera en su formulación la dependencia de
este parámetro de rigidez con el estado de esfuerzos del suelo, se consideró la medición
gráfica de los valores de este módulo para cada uno de los tres estados ensayados
(esfuerzos de confinamiento de 300, 350 y 400 kPa), y debido a que se está modelando el
comportamiento drenado del suelo, el módulo de Young o módulo de rigidez está en
términos efectivos.
B. Relación de Poisson υ
Se trata de uno de los parámetros más difíciles de evaluar en el contexto de la mecánica
de suelos, pero al mismo tiempo, de uno de los más relevantes en el análisis numérico por
elementos finitos y particularmente en el modelo Mohr-Coulomb. Normalmente, la
relación de Poisson efectiva puede ser evaluada a través de la realización de un ensayo
triaxial drenado, midiendo las distintas variables de deformación unitaria. Durante la
etapa de falla en los ensayos triaxiales CU no se presenta cambio volumétrico en el suelo,
por lo que la relación de Poisson evaluada en este tipo de ensayo corresponde al
parámetro no drenado υu (generalmente υu≈ 0,5). Dado que el modelo considerado en
este estudio requiere la introducción de este parámetro en términos efectivos, el uso del
valor obtenido del ensayo no drenado reproduciría de manera equivocada el
comportamiento del suelo.

METODO DE CULLMAN
Los muros son obras destinadas a la contención de tierras en general. En particular
pueden contener granos, agua, etc. como resulta evidente en los muros que se
encuentran a la intemperie, la lluvia se filtra a través de la tierra y entonces el
muro pasa a sostener los efectos de empuje dados por la tierra y por el agua, por
lo que habrá de tener en cuenta este factor en cuanto a su cálculo. la utilización de
muros de contención es muy frecuente en todo tipo de obras.
Empujes: se denomina empuje a la acción que las tierras ejercen sobre el muro.
Los empujes pueden ser activos y pasivos.

Empuje activo (Ea): es el que ejerce la tierra que es sostenida por el muro y que
para dicho fin se construye este.

Empuje pasivo (Ep): el empuje pasivo contrarresta la acción del empuje activo, y
es el producido por un terreno que absorbe la acción producida por la estructura.
Este método puede utilizarse en muchos tipos de muros y con muchos tipos de
sobrecarga.

El método de cullman sigue el siguiente proceso:

a) define la línea del talud natural como la que partiendo del vértice B del trasdós
del muro, forma un Angulo φ (que es el de rozamiento interno del terreno) con la
horizontal.

b) define la línea de dirección como aquella que pasando por B forma un Angulo φ
+ δ con el paramento del muro.
El método de cullman dice que si a partir del punto B, que hemos considerado
como origen de coordenadas, llevamos sobre la línea de talud natural BD, la
magnitud del peso del prisma ABC a una determinada escala, nos dará el punto J .
Si ahora por JB se traza una paralela a la línea de dirección cortara a la línea BC en
el punto N. Este valor JN representa a la escala indicada para la fuerza el valor del
empuje activo producido por el prisma ABC.

BJ= valor del peso del prisma ABC.

JN= valor del empuje sobre el muro producido por el prisma ABC.

En vista de esto podemos indicar que como lo que aquí se pretende es el determinar el
empuje máximo (Ea) se consideraran tantos puntos C como sean necesarios para así
describir una curva en la que podamos determinar el Ea máximo.
EJEMPLO DE METODO DE CULLMAN:
Sea un muro de contención tal como el indicado en la figura anterior para contener unas
tierras de peso específico γ =2000 kg/m3. De φ= 45 grados y δ=10 grados. Haremos los
cálculos por metro lineal de muro. (Para 1 m de ancho unitario)
Hemos dividido el terreno en seis partes, cuando en realidad y si es que se pretende
mayor aproximación se puede dividir en tantas como se quiera, no obstante en líneas
generales una división de ocho partes γc suele dar buena aproximación.
Una vez determinados los pesos de los diferentes prismas de terreno se llevan a la escala
indicada de fuerzas y sobre la línea de talud natural los valores de las fuerzas, obteniendo
así los puntos, J1, J2, J3, J4, J5. Por cada uno de estos puntos se trazan paralelas a la línea
de dirección hasta que corten las líneas BC1, BC2, BC3, BC4, BC5 obteniendo así los puntos
N1, N2, N3, N4, N5. Estos últimos puntos se unen mediante una curva tal como se indica
en la figura. Construida la curva se traza una tangente a dicha curva que sea paralela a la
línea del talud natural tal como se indica. Por el punto de tangencia se traza ahora una
paralela a la línea de dirección obteniéndose de esta forma los puntos J" y N" (primas). El
empuje sobre el muro vendrá dado por el valor del segmento J"N" medido con la escala
de fuerzas. En el ejemplo citado J"N"= 2.1cm por lo que el valor del empuje sobre dicho
muro bale:
Empuje= 2.1 x 3 = 6.3 toneladas.