Los análisis sísmicos de túneles han sido tradicionalmente abordados mediante expresiones analíticas para geometrías sencillas que no incluyen las secuencias constructivas ni historiales de esfuerzos (Wang, 1993; Penzien y Wu, 1998; Penzien, 2000). Últimamente, algunos softwares de análisis geotécnico han entregado herramientas para la resolución de problemas complejos, permitiendo incorporar las variaciones en los historiales de tensiones, métodos constructivos, secuencias de excavación y solicitaciones sísmicas a través de registros de aceleraciones.
Este artículo presenta un estudio comparativo de 3 métodos de análisis sísmico para un túnel NATM construido en suelos finos del noroeste de Santiago. Se describe la metodología, consideraciones particulares y los parámetros empleados en cada caso. Se indican las complejidades y los tiempos computacionales requeridos para el desarrollo de cada metodología. Finalmente, se presenta un análisis comparativo de los resultados obtenidos: esfuerzos sísmicos en revestimiento del túnel y cálculo de espesor de revestimiento.
Geometría del túnel y propiedades del suelo de fundación
La geometría del túnel se muestra en la Figura 1. La secuencia constructiva considera 3 secciones principales:side drift I, sección central, side drift II y 9 subsecciones que se enumeran en la misma figura. La metodología utilizada es acorde a los principios del método NATM (New Austrian Tunnelling Method) y simula las secuencias de excavación en tres etapas constructivas: bóveda, banco y contrabóveda; con desfase entre etapas y la aplicación de revestimiento estructural. Entre cada frente de avance de sidedrift, hay un desfase de 10 m, así también de la pared central. El nivel de riel del túnel se encuentra a una profundidad de 22 m del nivel de terreno y la clave del túnel se encuentra a 16 m de profundidad. La sección del túnel abarca un área aproximada de 190 m2.
Figura 1: a) Geometría (dimensiones en cm) y b) secuencia constructiva
del túnel
Se considera un túnel construido en el sector de suelos finos del noroeste de Santiago, al cual se le han asignado las propiedades geotécnicas presentadas en la Tabla 1. El módulo de deformación ha sido considerado lineal aumentando en profundidad, también se han considerado distintos valores de cohesión y coeficiente de empuje en reposo in situ K0 para dos distintos estratos de suelo.
Tabla 1: Propiedades de los materiales (ARCADIS, 2014)
Z: profundidad medida desde la superficie en m
Dj: profundidad sello fundación en m, B: dimensión menor de
estructura en m
Solicitación sísmica
Con el fin de simular la solicitación sísmica, se utilizan dos procedimientos: desangulación sísmica y análisis dinámico con registro de aceleraciones. Para la desangulación sísmica, la metodología empleada se basa en las recomendaciones del Manual de Carreteras (2014), que se sustentan en la propuesta de Kuesel (1969) para el diseño sísmico del metro de San Francisco. En este estudio se ha considerado una desangulación θs de 1.1·10-3rad, obtenida de los valores tabulados en el Manual de Carreteras (2014) para un rango de compresión no confinada qu entre 20 y 40 kPa, para zona sísmica con ao= 0.4g.
El análisis dinámico se basó en uno de los registros de aceleraciones del terremoto de Chile, ocurrido el 27 de Febrero del 2010, que tuvo una magnitud momento Mw de 8.8. El sismo fue subductivo tipo thrust con epicentro marítimo frente a la localidad de Cobquecura, Región del Bío Bío (Saragoni y Ruíz, 2012). El registro de aceleraciones fue obtenido de la Red Nacional de Acelerógrafos de la Universidad de Chile (RENADIC). Corresponde a un registro de superficie con componente horizontal, obtenido en una estación ubicada en Maipú, sobre depósitos de ceniza volcánica denominados comúnmente como "Pumicita". Las principales características del registro de aceleraciones se indican en la Tabla 2.
Tabla 2: Principales características sísmicas registro aceleraciones
terremoto 2010, estación Maipú (Saragoni y Ruíz, 2012)
La Figura 2 presenta las componentes de aceleraciones, velocidades y desplazamientos del registro utilizado. El registro ha sido sometido a corrección de línea de base. La Figura 3 presenta los espectros de Fourier y pseudo-aceleración para un amortiguamiento del 5%.
Figura 2: Registros de aceleración, velocidad y desplazamiento en función del tiempo.
Sismo 27F2010, estación Maipú
Figura 3: Espectro de Fourier y espectro de respuesta de aceleraciones
(5% de amortiguamiento)
