MODELO CONSTITUTIVO UBC SAND

El Modelo UBCSAND es un modelo bidimensional de plasticidad en tensiones efectivas para el

análisis de la licuefacción de suelos. El modelo predice básicamente el comportamiento tensión-
deformación de corte del suelo empleando la relación hiperbólica de Duncan-Chang y estima la
respuesta volumétrica asociada del esqueleto del suelo usando una regla de flujo que es una
función de la relación de esfuerzos. El modelo UBCSAND ha sido desarrollado por la University
of British Columbia por el Prof. P.M. Byrne (1998). La formulación del modelo UBCSAND
presenta diversas investigaciones realizadas por Puebla et al. (1997), Beaty & Byrne (1998),
Puebla (1999), Naesgaard et al. (2009.

COMPORTAMIENTO ELÁSTICO NO LINEAL

Para explicar y predecir el comportamiento elástico de materiales granulares constituidos por

partículas discretas, como la arena, la mecánica de suelos ha recurrido tradicionalmente a la
teoría de la elasticidad que rige la mecánica de los medios continuos. Por lo tanto, el
comportamiento elástico de la arena queda caracterizado mediante dos parámetros elásticos.
Aunque existen varias posibles elecciones de este par de parámetros elásticos, el modelo
UBCSAND emplea el módulo volumétrico elástico Ke y el módulo de Poisson ν.

El comportamiento elástico que ocurre dentro de la superficie de fluencia está gobernado por
el modelo de elasticidad no lineal descrito por Puebla (1999). Se asume un módulo volumétrico
elástico Ke dependiente de la presión y la relación viene dada por la ecuación:

𝑝
𝐾 =𝐾 𝑝
𝑝

Como el módulo de Poisson permanece constante, esto implica que el módulo de corte elástico
Ge depende de la presión, a través de Ke, tal como se observa en la siguiente ecuación:
3𝐾 (1 − 2𝜈)
𝐺 =
2(1 + 𝜈)
donde Pref es la presión de referencia, que habitualmente se toma como la presión atmosférica.
KBe es el módulo volumétrico elástico de referencia (esto es el módulo volumétrico de referencia
dividido entre la presión de referencia), p es la presión en el suelo, me es el parámetro utilizado
para definir el grado de dependencia de presión con la rigidez.

El módulo de corte elástico de referencia KGe (esto es el módulo de corte de referencia dividido
entre la presión de referencia), se puede obtener a partir de la siguiente ecuación:

𝑝
𝐺 =𝐾 𝑝
𝑝

donde ne es el parámetro utilizado para definir el grado de dependencia de presión con la

rigidez.

COMPORTAMIENTO PLÁSTICO

El comportamiento plástico de los materiales granulares constituidos por partículas discretas,

como la arena, es se ha modelizado tradicionalmente haciendo uso de la teoría de la plasticidad
que rige la mecánica de los medios continuos. De tal forma que, macroscópicamente el material
granular puede ser considerado como un continuo en el cual es posible definir un estado de
esfuerzos interno, además de unas deformaciones asociadas.

En el modelo UBCSAND se consideran de forma desacoplada el comportamiento plástico del

modelo, la resistencia al corte y la resistencia a compresión. Esto implica considerar en el modelo
doble endurecimiento, y que el comportamiento plástico es una combinación de la rotura por
cortante y una deformación irrecuperable.

Funciones de fluencia

Las funciones de fluencia del modelo UBCSAND definen tres superficies de fluencia:

- una superficie de rotura por cortante similar a la de Mohr-Coulomb.

- una capa esférica como superficie de rotura por compresión.
- una superficie de rotura por tracción (tension cut-off surface).

La superficie de rotura por cortante es una aproximación de la bien conocida superficie de Mohr-
Coulomb, donde los vértices son reemplazados por una superficie continua. Por otra parte, la
formulación del plano de la superficie de rotura por tracción es similar a la de los modelos de
Mohr-Coulomb y Drucker-Prager. Las superficies de rotura por cortante y por compresión son
definidas por las funciones f1 y f2, respectivamente, en el espacio p-q.

La función f1 viene dada por la siguiente ecuación:

𝑞 6 sin(𝜙)
𝑓 = − (𝑝 + ∆𝑝) = 0
𝑅 (𝜃) 3 − sin(𝜙)
donde φ es el ángulo de rozamiento movilizado en compresión triaxial, Δp es el cociente entre
la cohesión c y tan(φ), y que permite modelizar el comportamiento cohesivo del material, R1(θ)
permite considerar las diferencias existentes entre la compresión triaxial y la extensión triaxial
y que es función del ángulo de Lode, θ.

La capa elíptica f2 viene dada por la siguiente ecuación:

𝑞
𝑓 = (𝑝 + ∆𝑝) + 𝛼 −𝑝 =0
𝑅 (𝜃)
donde pc es la presión de preconsolidación y α determina la forma de la capa elíptica. Como en
la expresión anterior se emplea R2(θ) para modelizar las diferentes resistencias en compresión
triaxial y extensión triaxial. En la figura aparece el modelo en el espacio p-q:

Superficies de fluencia por cortante y compresión en el espacio p-q para θ=π/6

En el plano desviador, existen varias posibilidades para considerar la resistencia de las arenas.
Una de ellas es el criterio de rotura de Mohr-Coulomb, que ha demostrado ser preciso para
ensayos triaxiales en arenas (Groen, 1995), que es conceptualmente simple y muy adecuado
para cálculos manuales por el otro, aunque en los métodos de elementos finitos, las trayectorias
de tensiones no son ni conocidas de antemano ni constantes. La utilización del criterio de Mohr-
Coulomb en un espacio continuo tridimensional plantea algunas dificultades que se unen al
hecho de emplear múltiples superficies para el criterio de fluencia.

Por lo tanto, el objetivo es definir un modelo que sea continuo, y capaz de representar con
precisión la resistencia de la arena bajo varias trayectorias de esfuerzo triaxial, para ello se
adopta el modelo de Eekelen (1980), que proporciona una solución simple y fácilmente aplicable
a procedimientos de elementos finitos (Groen, 1995).

El modelo puede ser ajustado sobre la superficie de fluencia del modelo Mohr-Coulomb en el
plano desviador, mediante las funciones R1(θ) y R2(θ), donde θ es el ángulo de Lode que viene
dado por la siguiente ecuación:
27𝐽
cos(3𝜃) = −
2𝑞
De esta forma, el esfuerzo q puede ser asimilado al radio en el plano desviador y el ángulo de
Lode θ al ángulo, así se tiene un sistema coordenado polar en el plano desviador con
coordenadas (q,θ), tal como se observa en la figura.
Definición del ángulo de Lode y envolvente última respecto a la convexidad (adaptado de Groen,
1995)

El parámetro α determina la forma de la capa elíptica, y viene dado por la siguiente ecuación:

𝛼 = (1 − 𝛽)
donde α, β y n son constantes.

Tal y como se ha indicado anteriormente el ajuste al modelo Mohr-Coulomb se realiza mediante

las funciones R1(θ) y R2(θ). La primera viene dada por la ecuación:

1 − 𝛽 cos(3𝜃)
𝑅 (𝜃) =
1−𝛽

donde β1 viene dado por la expresión:

3 + sin(𝜙)
−1
3 − sin(𝜙)
𝛽 =
3 + sin(𝜙)
+1
3 − sin(𝜙)
De acuerdo con Eekelen (1980) el mejor valor de n es -0,229. Según Groen (1995), n=-0,229 toma
el mismo valor que en el caso de Eekelen y tiene un valor máximo β≤0,7925, que se corresponde
con un ajuste al modelo Mohr-Coulomb con un ángulo de fricción φ=46.55ᵒ. Otra posibilidad
sería tomar un valor de n=-0,25 y β≤0,75592, que corresponde a un ajuste al modelo Mohr-
Coulomb con un ángulo de fricción φ de46.99ᵒ Para β=0 la superficie de fluencia se reduce a la
superficie de fluencia de Drucker-Prager mostrada en la figura.

Superficie de fluencia en el plano desviador

La función R2(θ) modifica la forma de la capa de la superficie de fluencia de compresión y viene
dada por la siguiente expresión:

1 − 𝛽 cos(3𝜃)
𝑅 (𝜃) =
1−𝛽

Tomando un valor de n=-0,229 y siendo β2=0 se obtiene una capa esférica. Por lo tanto, los
factoresβ1 y β2 son los parámetros del contorno de fluencia.

Finalmente, la superficie de rotura por tracción se define mediante la función f3 mediante la

siguiente ecuación:

𝑓 =𝑝+𝑝 =0
donde pt es el esfuerzo máximo a tracción.

Regla de flujo

Una superficie de potencial plástico es aquella que va cambiando según la regla de flujo plástico,
que es la que fija la dirección del vector incremento de deformación plástica en el espacio de
direcciones principales, mediante una función que liga sus componentes con las tensiones.

En caso del modelo UBCSAND se consideran dos superficies de potencial plástico. Para el caso
de la rotura por cortante se supone que el flujo plástico es no asociativo y se determina mediante
la superficie g1:
6 sin(𝜓)
𝑔 =𝑞− (𝑝 + ∆𝑝) = 0
3 − sin(𝜓)
El ángulo de dilatancia ψ se relaciona con el ángulo de rozamiento φ mediante una regla basada
en la energía similar a la teoría de esfuerzo-dilatancia desarrollada por Rowe (1962) y linealizada
y simplificada de acuerdo con consideraciones de energía, tal como lo propuso Puebla (1999),
según la expresión:

sin(𝜓) = sin(𝜙) − sin(𝜙 )

donde φcv es el ángulo de rozamiento a volumen constante.

El flujo plástico en la capa de compresión se determina mediante la superficie de potencial

plástico g2:

𝑔 = (𝑝 + ∆𝑝) + 𝛼 𝑞 − 𝑝 = 0

Que supone flujo asociativo en el espacio p-q y flujo no asociativo en el plano desviador,
basándose en los mismos argumentos que en el caso de la regla de flujo para la superficie de
rotura. Por lo tanto αg=α, por ser flujo asociativo (ψ=ϕ).

Finalmente, para la superficie de rotura por tracción se considera una regla de flujo asociativo.

Mecanismo de endurecimiento monotónico

El modelo UBCSAND, considera el endurecimiento independiente de la superficie considerada.

Así el endurecimiento por cortante está delimitado por la línea de rotura (ángulo de rozamiento
último) tal como se muestra en la figura.
Mecanismos de endurecimiento monotónico de las superficies de corte y compresión (adaptado
de Rahman & Schreppers, 2014).

Para la superficie de rotura por cortante se supone que existe una relación hiperbólica entre la
relación de esfuerzos presente en el suelo η y la deformación plástica de corte equivalente tal y
como aparece en la figura.

Relación entre la relación de esfuerzos, η, y la deformación plástica de corte ks (adaptado de

Petalas & Galavi, 2012)

La formulación ha sido desarrollada de acuerdo con las consideraciones de Puebla et al. (1997)
y viene dada por la siguiente expresión:

Δ𝜂 = 𝐺 Δ𝑘

donde Δ𝜂 es el incremento de la relación de esfuerzos y Δks es el incremento en la deformación

plástica de corte equivalente y Gsp es el módulo de corte plástico tangente normalizado dado
por la expresión:

𝑝 𝜂
𝐺 =𝐾 1− 𝑅
𝑝 𝜂

Donde KGp es el módulo de corte plástico, p es la presión, pref es la presión de referencia, np es

el exponente del módulo de corte plástico, ηf es la relación de esfuerzos en rotura, y Rf es la
relación en rotura. La relación de esfuerzos presentes en el suelo viene dada por la siguiente
ecuación:
𝑞
𝜂=
𝑝
El incremento de la deformación plástica de corte equivalente Δks viene dada por la expresión:

∆𝑘 = ∆𝜀 − ∆𝜀

donde Δε1p y Δε3p son los incrementos de las deformaciones plásticas principales mayor y menor,
en relación con la superficie de rotura por cortante.

Para el endurecimiento de la superficie de rotura por compresión, la presión de

preconsolidación viene dada por una función de un parámetro interno k2 mediante la siguiente
expresión:

𝑝 = 𝑝 (𝑘 )
Estando k2 definido como:

∆𝑘 = −∆𝜀

La variación de pc está gobernada por la siguiente ley exponencial:

1+𝑒
𝑝 = 𝑝 𝑒𝑥𝑝 ∆𝑘
𝛾
en la que pc0 es la presión de preconsolidación inicial, y ϒ es el factor de saturación que puede
ser considerado como un parámetro del material.

Finalmente, para la superficie de rotura por tracción no presenta endurecimiento.

Mecanismo de endurecimiento cíclico

Para simular la respuesta bajo carga cíclica en un estado de esfuerzos tridimensional se emplea
un mecanismo de endurecimiento isotrópico-cinemático mixto basado en un criterio de
superficie de delimitación (bounding surface approach), de acuerdo con lo planteado por
Dafalias & Popov (1976), Shen et al. (1993), Bakmar et al. (2008). Sin embargo, en el modelo
UBCSAND se propone que este mecanismo de endurecimiento llegue a ser idéntico al
mecanismo de endurecimiento descrito por Puebla (1999) en un estado bidimensional de
esfuerzos. El tensor de tensiones σ en la superficie de rotura por cortante es reemplazado por ξ
que se define como:

𝜉 = 𝛔 − 𝑝𝛂
donde p es la presión, σd es la parte desviadora del tensor de tensiones y α es el tensor de
cambio de esfuerzo (back stress). De acuerdo con la ecuación anterior, la variación del cambio
de esfuerzo α provoca un efecto de endurecimiento cinemático que mueve la superficie de
rotura por cortante alrededor de su vértice, tal como se muestra en la figura.
Superficie de fluencia tipo cono definido en términos de α (tomado de Bakmar et al., 2008).

El módulo de corte plástico KGp es aplicable para un primer tiempo carga, es decir siempre que
el incremento de carga expanda la superficie de rotura por cortante fuera de la anterior relación
de esfuerzos máxima η. Por otra parte, cuando un incremento de carga ocurre dentro del límite
de la anterior relación de esfuerzos máxima y mueve la superficie de rotura por cortante junto
con él (es decir mecanismo de endurecimiento cinemático), se asume que la arena se comporta
plásticamente, pero con un módulo de corte plástico que es varias veces más rígido que para un
primer tiempo de carga tal como lo mencionado por Beaty & Byrne (2011). Este fenómeno se
conoce como la densificación del suelo.

La rigidización de KGp se considera como una función del número de ciclos experimentados por
la arena usando un enfoque similar al empleado en el código de elementos finitos descrito por
Petalas & Galavi (2012):
𝑛−1
𝐾 =𝐾 4+ 𝑓𝑎𝑐
2
donde KG0p es el módulo de corte plástico inicial, n es el número de semiciclos contados desde
el comienzo de la carga cíclica y facpre es un factor de calibración para la regla de esfuerzos de
densificación. KGp puede llegar a ser muy elevado para un valor elevado de n. Por lo tanto, KGp
se limita mediante KGmaxp de tal manera que:

𝐾 = 𝐾 60 0,003 + 100

Finalmente, cuando se alcanza la máxima relación de esfuerzos posible ηf, se produce la

licuefacción y entonces KGp se reduce a:

𝐾 =𝐾 𝑓𝑎𝑐

Donde facpos es un factor de post licuación.

COMPORTAMIENTO NO DRENADO

Si en un análisis geotécnico la permeabilidad de un suelo saturado es baja y la velocidad de la

aplicación de la carga es rápida, el comportamiento a corto plazo se caracteriza por no
producirse cambio de volumen, predominando la generación de las presiones intersticiales
frente a la disipación. Este comportamiento no drenado puede ser modelizado definiendo un
exceso de presión de poro del fluido (𝑝𝑒) en el material como:
 𝑝 = −𝐾 𝜀

siendo εv la deformación volumétrica y Kf el módulo de compresión no drenada que es igual al

módulo de compresión drenada K multiplicado por un factor de penalización.

𝐾 = 𝑓𝑎𝑐 𝐾

PARÁMETROS DEL MODELO UBC SAND

PARÁMETROS ASOCIADOS A LA RESPUESTA ELÁSTICA

Presión de referencia

El parámetro Pref es la presión de referencia. Es habitual tomar el valor de la presión atmosférica

(Pref=100KPa). Igualmente ha sido empleado por diversos investigadores tales como Puebla
(1999), Tsegaye (2010), Jalili et al. (2012), Petalas & Galavi (2012) y (2013), Ju & Vassalos (2015),
Borowiec & Stanuszek (2016), entre otros.

Exponentes elásticos

Los exponentes elásticos me y ne son constantes de ajuste. Los valores del exponente elástico
ne, varía entre 0,38 y 0,5. Según Puebla (1999), Byrne et al. (2004), Beaty & Byrne (2011) este
parámetro puede ser aproximadamente 0,5.

Por otro lado, en diversos trabajos de investigación, tales como Beaty & Byrne (1998), Tsegaye
(2010), Jalili et al. (2012), Petalas & Galavi (2012) y (2013), Ju & Vassalos (2015), Borowiec &
Stanuszek (2016), entre otros, se toma me=ne=0,5.

Coeficiente de Poisson

El coeficiente de Poisson de la arena ν permanece constante en la formulación del modelo

UBCSAND, lo cual implica que el módulo de corte Ge depende de la presión según la ley de
elasticidad no lineal del modelo.

Según estudios realizados por Hardin (1978) el módulo de Poisson, para la arena varía entre 0 y
0,2, recomendando un valor de 0,12. Este valor de 0.12 hace que el valor del módulo volumétrico
elástico Ke sea prácticamente igual al del módulo de corte elástico Ge, según la siguiente
ecuación:
𝐺 𝐾 3(1 − 2𝜈)
= =
𝐾 𝐾 2(1 + 𝜈)

En ausencia de datos, se puede suponer el mismo valor del módulo volumétrico elástico y del
módulo de corte elástico según Puebla (1999), lo cual se cumple para un coeficiente de Poisson
de 0,125. Este valor se ha utilizado en diferentes trabajos publicados por Puebla et al. (1997),
Beaty & Byrne (1998), Jalili et al. (2012), Petalas & Galavi (2012) y (2013), Ju & Vassalos (2015),
entre otros.

Módulo de corte elástico

Según estudios realizados por Hardin (1978) el módulo de corte elástico Ge es función del índice
de huecos e, y de una variable A que relaciona el tamaño y la forma de las partículas de la arena,
de la siguiente manera:

𝐺 = 𝐴 𝑓(𝑒)
Los valores para A, f(e) y el exponente elástico ne sugeridos por varios autores son los que
aparecen en la siguiente tabla.

Material A f(e) ne Referencia

Arena limpia Yu & Richart,
700 (2,17-e)2/(1+e) 0,5
redondeada 1984
Arena limpia Yu & Richart,
326 (2,97-e)2/(1+e) 0,5
angulosa 1984
Hardin &
Arena limpia 320 (2,97-e)2/(1+e) 0,5
Drnevich, 1972
Arena limpia 625 0,3+0,7e2 0,5 Hardin, 1978
Iwasaki &
Arena limpia 900 (2,17-e)2/(1+e) 0,38
Tatsuoka, 1977
Valores típicos de A, f(e) y ne (adaptado de Puebla, 1999).

El módulo de corte elástico Ge también puede ser aproximado a partir de valores del número de
golpes normalizado del ensayo de penetración estándar (Byrne et al., 1987), o de la prueba de
penetración de cono (Robertson et al., 1983). Los valores aproximados de Ge a partir de la
densidad relativa, Dr, o la resistencia a la penetración estándar normalizada, (N1)60, aparecen en
la siguiente tabla.

Estado Dr (%) (N1)60 Ge

Muy suelto 0-15 0-3 300
Suelto 15-35 3-8 300-900
Medio 35-65 8-25 900-2000
Denso 65-85 25-42 1500-2500
Muy denso 85-100 >42 >2000
Valores aproximados del módulo de corte elástico Ge (adaptado de Puebla, 1999).

Por otro lado, el módulo de corte elástico Ge tambien puede obtenerse directamente de los
ensayos de laboratorio usando “bender elements” o un aparato de columna resonante.

Módulo volumétrico elástico

El módulo volumétrico elástico de la arena KBe, se puede obtener a partir de la consolidación

isotrópica en un ensayo triaxial drenado (CD) o a partir de un ensayo triaxial no drenado (CU)
con medida de la deformación volumétrica en la fase de aplicación de la presión de
confinamiento.

Según se ha visto anteriormente el módulo volumétrico elástico viene dado por la siguiente
ecuación:
∆𝑉
∆𝑝 = 𝐾
𝑉
Por otro lado el módulo volumétrico elástico no lineal, dependiente de la presión según Puebla
(1999) se puede expresar:

𝑝
𝐾 =𝐾 𝑝
𝑝

Para me=0,5, operándo se tiene:

 ,
𝑝 ∆𝑉
∆𝑝 = 𝐾 𝑝
𝑝 𝑉

Reordenando e integrando se obtiene la ecuación:

,
𝑑𝑝 𝐾 𝑝
= 𝑑𝑉
𝑝 , 𝑉

KBe puede ser obtenido mediante la siguiente expresión (*):

, ,
2𝑉 𝑝 2 𝑝
𝐾 = =
𝑉 −𝑉 𝑝 𝜀 𝑝

Donde (V0-Vf) es el cambio de volumen y εvc es la deformación volumétrica unitaria, resultantes

de la aplicación de la presión de confinamiento.

Por lo tanto, por lo indicado por Puebla (1999), Tsegaye (2010), Petalas & Galavi (2012) y (2013),
el módulo volumétrico de la arena KBe se puede obtener a partir de un ensayo de compresión
triaxial (CD o CU) mediante una presión de confinamiento igual a la presión de referencia
(p=σ3=pref=100KPa). Por lo tanto:
2𝑉 2
𝐾 = =
𝑉 −𝑉 𝜀

Cuando se dispone de los resultados de un ensayo triaxial a diferente presión de confinamiento,

esta debe ser corregida para obtener el módulo volumétrico elástico de la arena KBe.

PARÁMETRO ASOCIADO A LA RESPUESTA INICIAL DE LA ARENA

El ángulo de rozamiento inicial φ0 es un parámetro que indica la posición inicial de la superficie

de fluencia de cortante. Este ángulo para un suelo sin cohesión está relacionado con el estado
de esfuerzos inicial de acuerdo con la siguiente ecuación:
𝜎′ − 𝜎′
sin(𝜙 ) =
𝜎′ + 𝜎′
Para un ensayo triaxial con consolidación isotrópica, es decir σ’10=σ’30, el valor del ángulo de
rozamiento inicial es φ0=0. Por otro lado, cuando en el estado de esfuerzos inicial, σ’10≠σ’30 el
punto de partida de la trayectoria de tensiones seguirá el patrón indicado en la figura.

Movimiento del punto de inicio de la trayectoria de tensiones (tomado de Vaid et al., 2001).
PARÁMETROS ASOCIADOS A LA RESPUESTA PLÁSTICA DE CORTE

Los parámetros asociados a la respuesta plástica de corte pueden ser divididos en dos grupos:
los relacionados con la rigidez y los relacionados con la resistencia. En el primer grupo están el
módulo de corte plástico KGp y el exponente plástico de corte np, mientras que en el segundo
están el ángulo de rozamiento en rotura φf, la relación de rotura Rf y el ángulo de rozamiento a
volumen constante φcv.

Parámetros relacionados con la rigidez

Exponente plástico de corte

El exponente plástico de corte np, es una constante de ajuste para el comportamiento de

endurecimiento por corte, este parámetro puede ser tomado igual al exponente elástico ne,
según lo sugerido por Vermeer (1980).

Módulo de corte plástico

El modelo constitutivo UBCSAND considera en su formulación deformaciones elásticas y

plásticas. Restando la componente elástica de la deformación total de corte, la componente
plástica se puede obtener a partir de su formulación, ya que el módulo de corte G, está
relacionado con los módulos de corte elástico y plástico, Ge y Gp, como sigue:
1 1 ∆𝜏
∆𝛾 = ∆𝛾 + ∆𝛾 = ∆𝜏 + =
𝐺 𝐺 𝐺
donde Δϒ, Δϒe y Δϒp son los incrementos de deformación por cortante total, elástica y plástica
respectivamente. En consecuencia, operando se obtiene:

𝐺𝐺
𝐺 =
𝐺 −𝐺

De acuerdo con el modelo original desarrollado por Puebla (1999), para deformaciones totales
bajas, es decir, ϒ<10-4, G es igual al módulo de corte inicial Gi y Gp es igual al módulo de corte
plástico inicial Gip, según se ve en la figura.

Módulo de corte plástico inicial Gip (tomado de Puebla, 1999).

Por lo tanto, el módulo de corte plástico KGp puede ser estimado a partir de Gip utilizando la
siguiente ecuación:

𝑝
𝐺 =𝐾 𝑝
𝑝

Según Puebla (1999) en ausencia de otros datos, KGp puede ser estimado dentro del rango de
⅓KGe a KGe.

Parámetros relacionados con la resistencia

Ángulo de rozamiento en rotura

El ángulo de rozamiento en rotura φf es el ángulo de rozamiento máximo posible de acuerdo

con la curva hiperbólica de endurecimiento.

El ángulo de rozamiento en rotura se puede obtener a partir de la pendiente de la línea de estado

crítico (CSL) empleando la ecuación:

6 sin 𝜙
𝑀=
3 − sin 𝜙

La línea de estado crítico se puede obtener experimentalmente a partir de un ensayo triaxial

drenado (CD) o a partir de un ensayo triaxial no drenado (CU) tal como se observa en las figuras.

Trayectoria de tensiones para el ensayo triaxial drenado

Trayectoria de tensiones para el ensayo triaxial no drenado.

Donde M es la pendiente de la línea de estado crítico.

Relación de rotura

La relación de rotura Rf es un parámetro de corte plástico relacionado con la resistencia última

de la arena. Este parámetro relaciona la relación de esfuerzos en rotura ηf y la relación de
esfuerzos último ηult a partir de la mejor hipérbola de ajuste de los datos de laboratorio, y se
define por:
𝜂
𝑅 =
𝜂
Dado que la relación tensión-deformación es hiperbólica, la curva de relación de esfuerzos se
acerca al valor último ηult asintóticamente. Según Puebla (1999), Rf podría estar en el rango de
0,5 hasta 0,98.

Ángulo de rozamiento a volumen constante

El ángulo de rozamiento a volumen constante φcv depende principalmente de la mineralogía

(Puebla, 1999). Este parámetro de corte plástico se puede obtener directamente a partir de
ensayos triaxiales drenados midiendo las deformaciones y los cambios de volumen en la fase de
corte, o a partir de ensayos no drenados en los que se midan las presiones de poro en lugar de
las deformaciones volumétricas. Este parámetro se puede estimar mediante la evaluación de la
relación de esfuerzos η correspondiente al punto en el que la pendiente de la curva εvp-ϒp o la
curva u-ϒp, llega a ser cero tal como lo indicado en las figuras.
Cálculo del ángulo de fricción a volumen constante φcv, a partir del ensayo de corte simple
drenado sobre arena Ottawa (adaptado de Park, 2005).

Cálculo del ángulo de fricción a volumen constante φcv, a partir del ensayo triaxial no drenado.
Data de (Phillips & Byrne, 1993).

Estrictamente, las curvas εvp-ϒp o u-ϒp, deben ser utilizadas para evaluar φcv. Sin embargo, ya que
en este nivel de deformación la componente elástica es mucho menor que la componente
plástica, el error incurrido por la utilización de la deformación total, puede no ser despreciable
según Puebla (1999).

PARÁMETROS ASOCIADOS A LA RESPUESTA PLÁSTICA DE COMPRESIÓN

Tres parámetros son necesarios para definir el mecanismo de la superficie esférica de fluencia,
la posición inicial de la superficie de fluencia pc, el módulo volumétrico plástico KBp y el
exponente plástico de compresión mp. Los parámetros KBp y mp deben ser obtenidos a partir
del análisis de ensayos de compresión isotrópica, ya que en este tipo de pruebas la respuesta
volumétrica es independiente de los efectos de corte según Puebla (1999).

Posición inicial de la superficie de fluencia de compresión

El parámetro pc define la posición inicial de la superficie de fluencia de compresión esférica y se

identifica con la presión de preconsolidación inicial tal como se muestra en la figura y toma un
valor: pc=pref=100KPa.
Posición inicial de la superficie de fluencia de compresión

Exponente plástico de compresión

La constante mp es el exponente del módulo volumétrico plástico. Tiene un rango de variación

entre 0<mp≤1. A partir de observaciones experimentales de las trayectorias de deformaciones
y tensiones en las pruebas realizadas a relación de esfuerzo constante, Vermeer (1980) sugirió
el uso de los mismos exponentes elásticos y plásticos mp=me=0,5.

Módulo volumétrico plástico

El módulo volumétrico plástico KBp es un parámetro de la ley de potencia de la curva de

endurecimiento de la superficie de fluencia esférica.
El módulo volumétrico plástico inicial Kip puede ser obtenido de una manera similar a Gip. El
módulo volumétrico inicial Ki se determina a partir de la pendiente de la curva εv-p’. Por lo tanto,
el módulo volumétrico plástico inicial Kip puede ser calculado a partir de la ecuación:

𝐾𝐾
𝐾 =
𝐾 −𝐾

Por lo tanto, KBp puede ser obtenido a partir de Kip, utilizando la siguiente ecuación:

𝑝
𝐾 =𝐾 𝑝
𝑝

PARÁMETROS ASOCIADOS AL COMPORTAMIENTO COHESIVO

Cambio de presión

El cambio de presión Δp es una constante que permite simular el comportamiento cohesivo del
material. Es un cambio de presión para la superficie de fluencia de cortante según se ve en la
figura.
Cambio de presión para la superficie de fluencia de corte Δ𝑝 (adaptado de Rahman &
Schreppers, 2014).

Para suelos no cohesivos este parámetro se toma igual a 0.

Este parámetro se puede obtener a partir del trazado de la línea de estado crítico (CSL). Se
obtiene experimentalmente el parámetro c* y seguidamente se obtiene la cohesión del suelo
empleando la ecuación:

3 − sin 𝜙
𝑐= 𝑐∗
6 cos 𝜙

Finalmente se usa la expresión:

Δ𝑝 = 𝑐 tan 𝜙

donde φf es el ángulo de rozamiento en rotura.

Parámetro de tracción

Pt es un parámetro que define una superficie de rotura por tracción, que varía entre 0≤Pt ≤ Δp.
Por lo tanto, para suelos no cohesivos este parámetro vale Pt=Δp=0.

PARÁMETROS ASOCIADOS AL COMPORTAMIENTO NO DRENADO

El modelo UBCSAND está basado en el comportamiento observado de la arena en ensayos de

laboratorio en condiciones de carga monotónica y cíclica. El modelo estima directamente la
respuesta del esqueleto del suelo para incrementos de carga. La respuesta del fluido de los poros
es acoplada a la respuesta del esqueleto del suelo mediante el módulo volumétrico del fluido.

Porosidad

La porosidad en un suelo es la relación del volumen de poros con respecto al volumen total. Este
parámetro puede ser determinado a partir de datos de laboratorio y está definido por la
siguiente expresión:
𝑉
𝑛=
𝑉
Donde Vv es el volumen de huecos y V el volumen total.
Parámetro B de Skempton

El parámetro B de Skempton se determina directamente en laboratorio y se define como, la

relación que existente entre el aumento de la presión de poros Δu3 y el aumento del esfuerzo
isotrópico de confinamientoΔσ3, según la siguiente ecuación:
∆𝑢
𝐵=
∆𝜎
Módulo volumétrico del fluido

Es posible predecir el comportamiento no drenado en el ensayo triaxial mediante la adición del

módulo volumétrico del fluido Kf a los demás parámetros obtenidos del modelo numérico del
comportamiento drenado, según lo señalado por Puebla (1999). El módulo volumétrico del
fluido depende del grado de saturación de la muestra. Por lo tanto, la predicción de la respuesta
no drenada requiere obtener el valor del parámetro B de Skempton para cada triaxial no
drenado. El cálculo del módulo se hace mediante la siguiente ecuación:
𝑛𝐾
𝐾 =
1
𝐵−1
siendo Ke es el módulo volumétrico elástico de referencia.

PARÁMETROS ASOCIADOS AL COMPORTAMIENTO CÍCLICO

El comportamiento cíclico puede ser caracterizado mediante los parámetros de pre y pos
licuación. Estos parámetros deben ser calibrados hasta que coincida con los datos
experimentales, tal como lo señalado por Rahman & Schreppers (2014).

Factor de pre licuación

facpre es el factor de densificación utilizado para la calibrar el comportamiento de pre licuación

en caso de carga cíclica. Es un multiplicador que controla el aumento del módulo de corte
plástico durante la carga secundaria y que puede tomar valores mayores de 0. Para el caso de
carga monotónica este parámetro toma el valor de 1. Para valores mayores de la unidad, el valor
del módulo de corte plástico KGp aumenta y el comportamiento se hace más rígido y para valores
menores de la unidad el valor del módulo de corte plástico se hace menor y el comportamiento
es menos rígido.

Factor de post licuación

Es un factor para el comportamiento de post licuación en caso de carga cíclica. Su rango varía
entre 0<facpos≤ 1. Para el caso de carga monotónica este parámetro toma valor unidad. Es un
factor que se puede usar para “ablandar” el comportamiento del suelo posteriormente a la
licuación. Si el parámetro es menor que uno, el comportamiento post licuación del suelo es
aproximado, por el contrario si el factor es igual a uno, entonces se toma un módulo de corte
plástico idéntico al utilizado para la carga primaria.

PARÁMETROS ADICIONALES

Para añadir un comportamiento cohesivo o adaptar la forma predeterminada de las superficies

de fluencia se puede especificar los siguientes parámetros adicionales.
Parámetros de contorno de fluencia

Los parámetros para el contorno de fluencia son 𝛽1 𝑦 𝛽2. El parámetro 𝛽1 es el parámetro de

ajuste para la superficie de fluencia de cortante en el plano desviador, que está por defecto
ajustado a Mohr-Coulomb. Para 𝛽1=0, la superficie se reduce a la superficie de fluencia de
Drucker-Prager. El parámetro 𝛽2 es el parámetro de ajuste para la capa de la superficie de
fluencia, por defecto igual a cero 𝛽2=0, que implica una capa esférica

Parámetro de rigidez adicional

Es el parámetro de rigidez adicional, que tiene un rango de variación (0≤𝑎𝑑𝑑𝑠𝑡𝑖≤1).