CÁLCULO ANALÍTICO ESPECTRO DE DISEÑO EN
SUPERFICIE
Ciudad objeto de estudio: Manila – Filipinas
Presentado por:
Juan Camilo Mejía Herrera
Sebastián Zapata Díaz
Profesor:
Juan Carlos Olarte
Octubre de 2023
TABLA DE CONTENIDO
1. PARÁMETROS DINÁMICOS........................................................................................4
2. SELECCIÓN DE SISMOS.............................................................................................5
3. ESTUDIO DE SITIO Y ESPECTRO DE DISEÑO.........................................................6
4. AMPLIFICACIÓN DEL MOVIMIENTO...........................................................................8
5. ANÁLISIS DE RESULTADOS.....................................................................................10
Referencias
Anexos
LISTA DE FIGURAS
Figura 1. Curvas de degradación de rigidez y amortiguamiento para la rigidez y
amortiguamiento en arena. Fuente: Barrientos (2014) …………………………………4 pág.
Figura 2. Curvas de degradación de rigidez y amortiguamiento para la rigidez y
amortiguamiento. Fuente: INGEOMINAS. (2005) ………………………………………4 pág.
Figura 3. Selección de sismos vs Espectros de Amenaza en Roca. Fuente: PEER…. 5
pág.
Figura 4. Espectro de respuesta en superficie para cada sismo y promedio. Fuente: Los
autores…………………………………………………………………………………………7 pág.
Figura 5. Espectro de diseño para caso específico de análisis. Fuente: Los autores…7
pág.
Figura 6. Función de Transferencia y RRS para sismo RSN 590. Fuente: Los autores 8
pág.
Figura 7. Función de Transferencia y RRS para sismo RSN 31. Fuente: Los autores9 pág.
Figura 8. Función de Transferencia y RRS para sismo RSN 1632. Fuente: Los autores.9
pág.
1. PARÁMETROS DINÁMICOS
El perfil de suelo analizado, posee dos tipos de suelo predominantes, el primero con un
espesor aproximado de 4 metros el cual corresponde a una arena y el segundo un
depósito aluvial fino. Para ambos, se presentan las respectivas curvas de degradación del
módulo de cortante y del amortiguamiento. [1] Barrientos (2014) establece unas curvas de
degradación para la rigidez y el amortiguamiento para arenas con condiciones muy
similares a la del estrato de análisis, estas se pueden apreciar en le Figura 1
Figura 1. Curvas de degradación de rigidez y amortiguamiento para la rigidez y
amortiguamiento en arena. Fuente: Barrientos (2014).
Para materiales finos (depósito aluvial), las curvas de degradación, se observan en la
Figura 2, la cuales fueron planteadas por [2] INGEOMINAS. (2005) para unas
características de suelo similares a la del perfil analizado.
Figura 2. Curvas de degradación de rigidez y amortiguamiento para la rigidez y
amortiguamiento. Fuente: INGEOMINAS. (2005).
2. SELECCIÓN DE SISMOS
Con base en el espectro de amenaza en roca por Norma consultado en la Asignación 1
de la materia, se consulta en la base de datos [3] PEER tres sismos de tal forma que
cubran de la mejor manera el contenido frecuencial del Espectro de Respuesta de la
Norma y que a su vez, los rangos de aceleración espectral se encuentren ajustados a los
valores del espectro de la Norma.
De acuerdo con el modelo de clasificación sísmica del [3] PEER, los tres sismos
escogidos para efectos del presente análisis, se identifican bajo los ID RSN 31, RSN 590
Y RSN 1632, los espectros de los sismos en superposición del espectro de diseño en roca
por norma de la ciudad de Manila – Filipinas, se pueden apreciar en la Figura 3.
Figura 3. Selección de sismos vs Espectros de Amenaza en Roca. Fuente: PEER.
3. ESTUDIO DE SITIO Y ESPECTRO DE DISEÑO
Para el modelado en el software DEEPSOIL, se tuvieron en la cuenta varias
consideraciones que parten de conceptos fundamentales de la ingeniería sísmica hasta la
escogencia de los diversos tipos de curvas de degradación y amortiguamiento que han
propuesto diferentes autores a través de los años.
A. PERFILES DE SUELO
Para el caso de estudio se tiene una columna de suelo la cual está dividida en
suelo arenosos y deposito aluvial fino, con un nivel freático presente a los 5.0m de
profundidad. Las características del perfil de suelo se pueden apreciar en el Perfil
8.
En el caso de la arena, se tiene un espesor de capa 4.0m, considerando la arena
como un material denso (ver propiedades).
Para el depósito aluvial fino, se tiene un espesor de capa 25m, que, de acuerdo
con sus características, se considera un suelo sobre consolidado, por ende, se
deduce que es un suelo que se ha asentado en el tiempo, aumentado la presión
efectiva entre sus partículas, generando de esta forma un suelo rígido.
La resistencia a cortante o “shear strength” τ para cada material fueron
seleccionados a propósito del estudio de suelos con curvas de degradación y
amortiguamiento, [2] INGEOMINAS, en el cual se encuentran perfiles de suelo con
algunas semejanzas a los perfiles de estudio que facilitan la viabilidad, como, por
ejemplo, además del tipo de material, la profundidad, índices de plasticidad,
densidades relativas, entre otros. Teniendo entonces para la arena τ =159.85 kPa .
, y para el depósito aluvial τ =188.29 kPa .
Al momento de realizar el modelado, el perfil de suelo se discretizo en múltiples
capas de 0.3m de altura cada una, esto con el fin de garantizar el buen tránsito de
la onda de corte mediante una sucesión de pasos pequeños, frente a las alturas
de perfiles que se tienen, simulando el método de elementos finitos en cuanto a la
eficiencia.
B. CURVAS DE DEGRADACIÓN Y AMORTIGUAMIENTO
Para suelo de arenoso se tomaron las curvas de Seed&Idriss, Barrientos (2014),
estipuladas para las arenas. En la gráfica se presentan tres tipos de curvas para
diferentes tipos de arenas, en este caso se escogió la curva de arena Upper
Bound, ya que corresponde a arenas densas o muy densas.
Para el depósito aluvial fino, se recurrió a la construcción de la curva por medio de
Darendeli, gracias a se tienen las características necesarias del suelo, como lo son
el índice de plasticidad y el OCR; adicionalmente para el cálculo del factor k 0, se
(1)
realizo con la siguiente ecuación:
k 0=[ 1−sen( φ) ]∗OCR
sen(φ)
donde φ es elángulo de fricción del material.
En valores típicos φ para materiales finos varía entre 25° y 35°, para este caso se
consideró φ=30°. Posteriormente se verifico la curva con las curvas de estándar de
Darendeli para finos con diferentes con índice de plasticidad 15 y diferentes OCR,
la curva de Idrris de 1990, considerando los resultados muy acordes ya que se
trabajan materiales muy parecidos al perfil de estudio.
Con todas las hipótesis planteadas, y una vez calculados los espectros para cada
uno de los tres (3) sismos ingresados, se procede a superponer los tres espectros
y uno adicional que corresponde al promedio de los valores de los tres espectros,
los cuales se evidencian en la Figura 4.
1.4
1.2
1
Aceleración (%g)
0.8
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Periodo (s)
RSN 31 RSN 590 RSN 1632 PROMEDIO
Figura 4. Espectro de respuesta en superficie para cada sismo y promedio.
Fuente: Los autores.
Con el valor promedio de los sismos y teniendo calculado el espectro en superficie
por norma, de acuerdo con NSCP (2015) [4] (ver anexo), se determina el espectro
para el caso específico estudiado, se decide que la meseta será uniforme, pero
aumentada para suplir la parte del espectro por norma que no está contemplado
en el espectro de diseño, además que el ajuste que se debe hacer es mínimo, lo
que no supone un gran cambio en el espectro inicial. Dichos espectros se
relacionan en la Figura 5.
1.4
1.2
1
Aceleración (%g)
0.8
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Periodo (s)
ESPECTRO NORMA EN SUPERFICIE PROMEDIO
DISEÑO ESPECIFICO
Figura 5. Espectro de diseño para caso específico de análisis. Fuente: Los
autores.
4. AMPLIFICACIÓN DEL MOVIMIENTO
Haciendo uso del software libre Transfer Function Tool, se procede a calcular la Función
de Transferencia para cada uno de los sismos extraídos de la base de datos PEER; se
crean los respectivos perfiles de suelo, ingresando los parámetros dinámicos calculados
en los puntos anteriores y los otorgados por la información de los perfiles de suelo. Para
cada caso, lo primordial es contrastar la Función de Transferencia calculada con la gráfica
RRS calculada en el punto anterior también para cada uno de los sismos. Como el
software no permite extraer archivo de puntos y comas para superponer en una sola
gráfica, se colocan ambas gráficas para ser analizadas
Las gráficas correspondientes al sismo RSN 590, se pueden apreciar en la Figura 6
RRS
12
10
Aceleración (%g)
8
6
4
2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Periodo (s)
RSN 590
Figura 6. Función de Transferencia y RRS para sismo RSN 590. Fuente: Los
autores.
De la misma forma, se presenta la misma información para los otros dos sismos,
en las Figuras 7 y 8.
RRS
12
Aceleración (%g)
10
8
6
4
2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Periodo (s)
RSN 31
Figura 7. Función de Transferencia y RRS para sismo RSN 31. Fuente: Los
RRS
12
Aceleración (%g)
10
8
6
4
2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Periodo (s)
RSN 1632
autores.
Figura 8. Función de Transferencia y RRS para sismo RSN 1632. Fuente: Los
autores.
5. ANÁLISIS DE RESULTADOS
Las curvas de degradación de la rigidez requieren de claridad en los conceptos
fundamentales y parámetros como las características físicas de los tipos de suelo,
densidades, índice de plasticidad, esfuerzos, peso específico, relación de vacíos, para
poder así realizar un análisis acertado en categorizar el suelo y construcción de estas
curvas. Las curvas de la arena y los finos presentan tempranas degradaciones, por las
pendientes que presenta, sin embargo, suelen tener gran energía para resistir
deformaciones medianas, así como sus curvas de amortiguación presentan esta misma
consistencia de tener poca disipación de energía en periodos cortos, pero una buena
energía de disipación para periodos largos.
La columna de suelo que fue el caso de estudio, cuenta con suelos que son consistentes
entre sí por sus propiedades, donde se tiene un perfil de suelo rígido en finos y denso en
las arenas, lo que proporciona que las ondas generadas por movimientos sísmicos viajen
por un suelo que conserva sus propiedades sirviendo de medio continuo. Las curvas
histeréticas nos muestran dos perfiles de suelo que podrían tender a la una falla súbita o
frágil debido a su rigidez, sin embargo, el nivel freático de esta columna de suelo se
encuentra a los 5 metros de profundidad, entre el nivel inferior de la arena y superior de
los finos, lo que podría, debido al movimiento, generar desconfinamiento de los
materiales, disminuyendo sus esfuerzos internos y permitiendo un poco de ductilidad a
medida que el tiempo pasa. Es decir, que sismos de cortas duraciones podrían ser fatales
para esta columna de suelo, mientras que los de larga duración podría permitir un
comportamiento controlado evitando, tal vez, fatalidades. (ver curva histerética anexos).
Al calcular las funciones de transferencia para el perfil de suelo indicado, se puede
apreciar que a pesar que las funciones se obtienen a partir de tres excitaciones sísmicas
muy distintas entre sí, los comportamientos son muy similares, puesto que se maneja la
misma morfología de la curva y las diferencias entre los picos no son tan relevantes, lo
anterior permite concluir que los parámetros dinámicos y mecánicos del suelo gobiernan
el proceso de cálculo para dichas funciones de transferencia, teniendo en cuenta la gran
incertidumbre que se tiene en este paso del análisis, puesto que los sismos ingresados a
las bases de datos, no necesariamente ocurrieron en la zona de análisis (Manila –
Filipinas) y los perfiles de suelos tampoco son representativos de la zona de estudio.
El espectro en superficie por norma, se ajustó bastante bien al espectro promedio de los
tres sismos analizados, la zona que se salió de la meseta no fue representativa de la
totalidad de rango frecuencial, generando que los sismos estuvieran bien enlazados a la
envolvente definida por la misma, sin embargo, este resultado aun sigue teniendo una alta
fuente de incertidumbre debido a la poca representatividad de los sismos y los perfiles de
suelos analizados
Un comportamiento similar en las tres funciones de transferencia, es que las mayores
amplificaciones, se concentran en las frecuencias mas bajas, situación que entra en
concordancia con los sismos, puesto que los mayores valores de aceleración también se
presentan en las zonas donde el periodo es pequeño.
Referencias
[1] Barrientos, A. (2014). CURVAS DE DEGRADACIÓN Y DAMPUNG EN
SUELOS CHILENOS [Tesis de pregrado, Universidad de los Andes. Repositorio
institucional.]
[2] INSTITUTO COLOMBIANO DE GEOLOGÍA Y MINERÍA – INGEOMINAS.
(2005). Análisis y modelación de la Respuesta Sísmica Local del subsuelo en
Santiago de Cali. Extraído de:
chrome-extension://efaidnbmnnnibpcajpcglclefindmkaj/[Link]
.co/B7/21003002515005/documento/pdf/[Link]
[3] PACIFIC EARTHQUAKE ENGINEERING RESEARCH CHENTER. (23 de
octubre de 2023). Shallow Crustal Earthquakes in Active Tectonic Regimes.
Extraido de: [Link]
[4] ASOCIATION OF STRUCTURAL ENGINEERS OF THE PHILIPPINES, INC.
(2015), National Structural Code of The Phiiippines. Extraído de:
[Link]
Anexo
Tabla: Parámetros sísmicos extraídos de Norma de Diseño Filipina
Zona sismica 4 Todo el país excepto Palawan
Factor de zona 0.4 Tabla 208-3
SD (Velocidad de onde de
Tipo de suelo corte de 180 a 360 m/s) Tabla 208-2
Coeficientes de sitio
Ca 0.44 Tabla 208-7
Cv 0.64 Tabla 208-8
Tipo de sismofuente A Tabla 208-4
Magnitud sísmica 7<M<8.4 Tabla 208-5
Cercanía de la sismofuente 2.1 a 5 km
Coeficiente de cercanía a la sismofuente
Na 1.2 Tabla 208.5
Nv 1.6 Tabla 208.6
Categoría de ocupación de V(Miscellaneous structures) Tabla 103-1
las estructuras
Factor de importancia 1 Tabla 208-1
Parámetro elástico de
2.9 Tabla 208-12
respuesta para tipo de
estructura (aporticada)
Figura anexo: Diseño típico de espectro de Respuesta de acuerdo con National Structural
Code of The Philippines (2015) y curva histerética análisis resutados.