The Effects of Induced Anisotropy on Seismic Stability of Clayey Slopes: Numerical Study

Resumen Técnico: Los Efectos de la Anisotropía Inducida en la Estabilidad Sísmica de Taludes Arcillosos

Planteamiento del Problema
Los depósitos de arcilla naturales exhiben frecuentemente una pronunciada anisotropía debido a la evolución de su estructura y cohesión durante la deposición, la consolidación unidimensional y la carga subsiguiente. En consecuencia, las propiedades del suelo, tales como la rigidez, el fluencia (yielding), la generación de presión de poros y la resistencia al corte no drenada, dependen tanto del nivel de esfuerzo como de la dirección de la carga. Esta característica es particularmente crítica en problemas de estabilidad de taludes sísmicos, donde las trayectorias de esfuerzos oscilan entre compresión, extensión y corte simple durante el movimiento telúrico. Los modelos constitutivos isotrópicos tradicionales, como el modelo Modified Cam Clay (MCC), asumen una fluencia isotrópica y no pueden representar directamente el endurecimiento rotacional, la variación de resistencia direccional o la anisotropía inducida. Omitir estos factores puede conducir a una simplificación excesiva del comportamiento del suelo bajo carga dinámica.

Metodología
Este estudio emplea un enfoque numérico para investigar los efectos de la anisotropía inducida en la respuesta sísmica de un talud arcilloso. La investigación compara dos marcos constitutivos:

1. Modelo Isotrópico: El modelo Modified Cam Clay (MCC).
2. Modelo Anisotrópico: El modelo Simple ANIsotropic CLAY (SANICLAY), el cual incorpora el endurecimiento rotacional para simular la evolución de la estructura de la arcilla.

El análisis se llevó a cabo utilizando el software FLAC 2D. Se implementó un modelo de usuario definido personalizado (UDM) utilizando el lenguaje Fish para simular el comportamiento constitutivo de SANICLAY, mientras que el modelo CAMCLAY estándar se utilizó para la comparación isotrópica. Para asegurar una comparación válida, los parámetros de SANICLAY fueron ajustados (estableciendo $\alpha=0$, $C_{san}=0$ y $x=1$) para transformar efectivamente el modelo en una formulación MCC isotrópica para el caso de control.

El modelo geométrico representa un talud arcilloso de 7 metros de altura con una pendiente de 7:4, asentado sobre roca madre a una profundidad de 4 metros. El talud fue sometido a un historial de aceleración dinámica aplicado en la base, caracterizado por una aceleración máxima de 2.85 m/s² durante una duración de 10 segundos (analizada durante 13 segundos). El estudio monitoreó desplazamientos laterales, tasas de deformación por corte (para identificar superficies de deslizamiento), redistribución de esfuerzos cortantes, módulos de compresibilidad y de corte, y la evolución de las superficies de fluencia en puntos específicos de "historia" a lo largo del talud.

Contribuciones Clave y Resultados
El análisis comparativo entre el modelo anisotrópico (SANICLAY) y el isotrópico (MCC) revela varias respuestas mecánicas distintas inducidas por la anisotropía:

• Desplazamiento: El modelo SANICLAY predice desplazamientos laterales significativamente mayores en comparación con el modelo MCC. Esta diferencia aumenta con el tiempo, siendo el desplazamiento máximo en la parte superior del talud en ambos casos, pero con una magnitud mayor en la simulación anisotrópica.
• Superficies de Deslizamiento: Los contornos de la tasa de deformación por corte máxima indican que las superficies de deslizamiento en el modelo anisotrópico son más extendidas que en el modelo isotrópico, lo que sugiere una zona de falla potencial más amplia.
• Esfuerzo Cortante: La anisotropía conduce a un aumento del esfuerzo cortante en todo el talud. El modelo SANICLAY muestra consistentemente valores absolutos de esfuerzo cortante mayores en comparación con el modelo MCC, observándose concentraciones de esfuerzo en la base y en el pie del talud.
• No Uniformidad de la Rigidez: La inclusión de la anisotropía resulta en distribuciones más no uniformes de los módulos de compresibilidad y de corte en varios puntos del talud. Los módulos en el modelo SANICLAY exhiben una mayor variabilidad espacial en comparación con el caso isotrópico.
• Evolución de la Superficie de Fluencia: El estudio rastrea los parámetros de la superficie de fluencia ($\alpha$ y $p_0$). Los resultados demuestran que la anisotropía inducida causa una reducción en el tamaño de la superficie de fluencia ($p_0$) en varios puntos del talud en comparación con el modelo isotrópico. Además, el parámetro anisotrópico $\alpha$ evoluciona durante la carga, representando la rotación y distorsión de la superficie de fluencia, un mecanismo ausente en el modelo MCC.

Significancia y Conclusión
El artículo concluye que omitir la anisotropía inducida constituye una simplificación excesiva del comportamiento del suelo, incluso bajo cargas sísmicas relativamente débiles. Los resultados numéricos indican que la anisotropía altera fundamentalmente las características de estabilidad sísmica de los taludes arcillosos al aumentar los desplazamientos, expandir las superficies de deslizamiento, elevar los esfuerzos cortantes y reducir el tamaño de la superficie de fluencia.

Los autores sostienen que para lograr una precisión apropiada en los análisis de ingeniería y aplicaciones que involucren taludes arcillosos bajo carga dinámica, es necesario considerar la anisotropía. El estudio valida el uso de modelos constitutivos anisotrópicos, específicamente el marco de trabajo SANICLAY, como un enfoque más físicamente significativo para capturar el complejo comportamiento de esfuerzo-deformación-resistencia de las arcillas naturales durante eventos sísmicos en comparación con los modelos isotrópicos tradicionales.