ESPECTRO DE DISEÑO Y EFECTOS DE SITIO

PARA SECTOR PUERTO MALLARINO.

Elaborado por:

ELKIN CAMPO MEDINA, JEYSSON STIVEN PORTILLO MELO

Profesor:

M. Sc. ALEJANDRO CRUZ ESCOBAR

Asignatura:

INGENIERÍA SÍSMICA

ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA

UNIVERSIDAD DEL VALLE
SANTIAGO DE CALI, COLOMBIA
JULIO DE 2019
__________________________________________________________________________________
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

INTRODUCCIÓN

Según el ​Servicio Geológico Colombiano y​ la ​Norma Sismo Resistente Colombiana del 2010
(NSR-10) la ciudad de Cali se encuentra en una zona de amenaza sísmica alta teniendo en cuenta este
factor y la importancia de la ciudad en la economía del país y la región además de su número de
habitantes y extensión geográfica fue necesario elaborar un estudio detallado el cual se realizó en
2004. Este estudio es la ​Microzonificación sísmica de Cali ​el cual dividió el territorio de la ciudad en
zonas regidas por características de comportamiento sísmico similar.
Para la ingeniería estructural es de vital importancia conocer el comportamiento del suelo ante un
evento sísmico y las implicaciones que podría tener la respuesta de éste en las estructuras construidas
en estas zonas.
En el presente trabajo se describe el procedimiento y elaboración de un espectro de diseño para una
zona en particular de la ciudad de Santiago de Cali. El sector para el cual se ha desarrollado el
espectro de diseño es Puerto Mallarino; este lugar se encuentra clasificado en la zona de llanura
aluvial y que está clasificada como zona geotécnica 8​ a y zona de respuesta sísmica 6 según la
Microzonificación Sísmica de Cali (​ MZSC).

PROCEDIMIENTO

Imagen # 1. Localización del sector de Puerto Mallarino

Fuente: Alcaldía de Santiago de Cali
A continuación, se muestra el perfil de suelo para la zona de estudio de acuerdo al ​Anexo 3a de la
Microzonificación sísmica de Cali.

Lugar / Perforación:​ Puerto Mallarino / P. Mallarino P10

Profundidad cuaternario / Terciario: ​1.1 Km / 2.0 Km

1
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Período Elástico:​ 1.7 Segundos

Zona Geotécnica:​ 8a. Llanura Aluvial.

Tabla # 1. Perfil de Suelo Puerto Mallarino Zona 8a. Llanura Aluvial.

Fuente: Anexo 3a. Microzonificación Sísmica de Cali.

En la ​Tabla # 1 se tuvo que hacer uso de las tablas: ​4.18 y 4.20 de la ​Microzonificación Sísmica de
Cali​ para obtener los valores de amortiguamiento para cada estrato.

Usando el programa ​EFESIO se ingresan los valores de cada estrato de suelo descrito en la Tabla #1.
Con los datos del perfil estratigráfico ingresados se procede a seleccionar un sismo en este caso se
determinó usar las aceleraciones del sismo de ​Pizarro-2004.​

Fecha - Hora Magnitud

(UTC) Departamento-Municipio Latitud Longitud Profundidad(Km​) MI
2004-11-15 Choco-Bajo Baudo 4.768 -77.724 Menos de 30 km 6.7
Tabla # 2 . Sismo de Pizarro
Fuente: Servicio Geológico Colombiano

Ingresado el sismo se procede a analizar los resultados de la función de transferencia proporcionada

por ​EFESIO donde se debe identificar una frecuencia fundamental. Con este valor podemos verificar
si efectivamente este perfil de suelo se ajusta de acuerdo a los datos experimentales de la zona de
estudio (microtrepidaciones). Para comprobar esto es necesario conocer el rango de periodo
fundamental de la zona de respuesta sísmica 6 (seis) y del sector de Puerto Mallarino, de acuerdo al

2
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

estudio de microtrepidaciones de la zona, y su respectivo rango de frecuencias los cuales se muestran

a continuación.

Zona 6 - Llanura Aluvial (Zona de Respuesta Sísmica)

Periodo Fundamental (s) 1.4-1.75
Frecuencia Fundamental
(Hz) 0.714-0.571
Zona 8a - Puerto Mallarino(Zona Geotécnica)
Periodo Fundamental (s) 1.7
Frecuencia Fundamental
(Hz) 0,588
Tabla # 3.Rangos de periodo y frecuencia para la zona de estudio
Fuente: anexo 3a e Informe 5.2 Microzonificación Sísmica de Cali.

Gráfica # 1. Función de transferencia en EFESIO (Sismo-Pizarro)

Fuente: propia

3
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

La frecuencia fundamental del suelo calculada por EFESIO, ​0.511 hz​, está por fuera del rango del
establecido por el estudio de microtrepidaciones. Lo anterior describe la no validez del perfil
estratigráfico del suelo por lo que este deberá ser ajustado. Variando la altura (h) y velocidad de suelo
(Vs) del estrato más profundo, se ajusta el modelo para que concuerde con los rangos ya definidos.

h (m) Vs (m/s)
Valores establecidos 145 410
Valores ajustados 115 410

Tabla # 4. Ajuste al modelo del perfil estratigráfico

Fuente: propia

Con este ajuste la frecuencia fundamental del suelo según ​EFESIO e​ s de ​0.585 hz que es muy cercana
a la frecuencia fundamental mostrada en la ​Tabla #3.

Ahora es necesario contrastar esta información con datos de sismos registrados en esta zona. Con el
sismo ya definido (​Pizarro-2004.​ ) se procede a evaluar las diferentes estaciones de acelerómetros
instaladas en Cali con el fin de elegir dos (2); una de estas estaciones tiene que estar en suelo rocoso
mientras que la otra debe estar en el lugar de estudio, o muy cerca de este, y sus aceleraciones deben
haber sido tomadas en suelo blando.

En este punto y con los datos de aceleración y tiempo de las estaciones escogidas se procede por
medio de MatLab a realizar la transformada rápida de Fourier tanto en los datos de entrada (Estación
en roca) como en los de salida (Estación en suelo blando) con el comando FFT.
Una vez realizada la transformada rápida de fourier para cada señal de aceleraciones se procede a
hacer una función de transferencia punto a punto.

4
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Gráfica # 2. Función de transferencia en MatLab (Sismo-Pizarro)

Fuente: propia

Se identifica que la frecuencia fundamental experimental difiere un poco con la frecuencia

fundamental identificada en el modelo de EFESIO. Es necesario entonces ajustar el perfil
estratigráfico de suelo para que este concuerde con la frecuencia fundamental encontrada
experimentalmente a través de la función de transferencia.

h (m) Vs (m/s)
Valores establecidos 145 410
Valores ajustados al análisis experimental 107 410
Tabla # 5. Ajuste al modelo del perfil estratigráfico
Fuente: propia.

Con el modelo ya ajustado es posible ingresar los sismos de diseño en EFESIO y de esta manera
obtener la respuesta en superficie del perfil para cada uno de estos eventos telúricos.

REGISTROS

Se utilizaron registros del ​Catálogo de Sismos ​correspondiente a cada país d​ e donde se obtuvo la
información necesaria para categorizar los diferentes sismos. Se usaron los parámetros mencionados
en el informe ​1.9 de la Microzonificación Sísmica de Cali​ para cada fuente sismogénica.

5
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

PAÍS MÉXICO MÉXICO PERÚ

Unidad de Instrumentación Sísmica del Instituto
Referencia de Ingeniería de la UNAM CISMID
Fecha 19-09-1985 21-09-1985 23-06-2001
Latitud 18.081 18.021 -16.08
Longitud 102.942 101.479 -73.77
Profundidad (km) 15 km 15 km 33 km
Mw 8.1 7.5 6.9
Mecanismo Inverso Inverso Normal
Código estación UNIO PAPN MOQ-001
Latitud estacion 17.987 17.328 -17.186
Longitud estacion 101.811 101.040 -70.928
Tasa de Muestreo 200 100 100
Distancia Epic.
(km) 120 89 290
Tabla # 6. Características de los sismos de diseño de subducción superficial
Fuente: Unidad de Instrumentación Sísmica del Instituto de Ingeniería de la UNAM y Centro Peruano Japonés de
Investigaciones Sísmicas y Mitigación de desastres.

PAÍS PERÚ CHILE JAPÓN

Referencia CISMID MASSONE K-NET 2005
Fecha 5-01-1974 13-06-2005 26-05-2003
Latitud -12.39 -19.89 38.9
Longitud -76.29 -69.13 141.6
Profundidad (km) 91.7 117 61
Mw 6.5 7.9 7.0

Mecanismo ND Inverso Inverso

Código estación ZAR Iquique MYG005
Latitud estacion -12.02 -20.22 38.8
Longitud estacion -77.01 -70.14 140.7
Tasa de Muestreo 50 200 100

Distancia Epic. 90 119 89

6
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

(km)
PGA -138.94 218.82 171.99
Tabla # 7. Características de los sismos de diseño de subducción Profunda
Fuente: Centro Peruano Japonés de Investigaciones Sísmicas y Mitigación de desastres

PAÍS NUEVA ZELANDA JAPÓN MÉXICO

Referencia Geonet Project 2005 K-NET 2005 UNAM
Fecha 17-08-2004 20-03-2005 30-09-1999
Latitud -42 33.7 15.95
Longitud 171.9 130 97.03
Profundidad (km) 15 12 16
Mw 5.8 6.6 7.5
Mecanismo Inverso Rumbo ND
Código estación 090A FKO001 OMTP
Latitud estacion -41.9 33.8 16.689
Longitud estacion 172.0 130.5 98.398
Tasa de Muestreo 50 100 1
Distancia Epic.
(km) 24 33 28
PGA 134.53 223.47 194.98
Tabla # 8. Características de los sismos de diseño corticales
Fuente: Centro Peruano Japonés de Investigaciones Sísmicas y Mitigación de desastres

Las especificaciones anteriormente descritas caracterizan a cada uno de los sismos, así como a la
estación que se usó para tomar los registros en aceleración de cada uno de estos movimientos
telúricos; cabe resaltar que la información de estas estaciones se obtuvo de catálogos de aceleración de
diferentes bases de datos de algunos países del mundo y en cada una de estas estaciones fue
primordial el verificar que los registros de la misma hayan sido tomados en roca, esto, con el
propósito de que los datos obtenidos no hayan sido afectados por amplificaciones propias de suelo
blandos.
Los registros de aceleraciones de estos sismos se escalaron respecto a la aceleración pico efectiva de
Santiago de Cali (Aa=0.25g) con el fin de que el máximo valor en su registro acelerómetro fuera
precisamente Aa. La componente que se utilizó en los registros de aceleración para cada uno de los
sismos fue la ​Este-Oeste ​(E-W).
Con los sismos de diseños escalonados respecto a Aa y el modelo del perfil de suelo estratigráfico
ajustado es posible, por medio del programa EFESIO, obtener la respuesta en superficie para cada uno
de los sismos de diseño. En este punto se tienen los datos necesarios para realizar un espectro de

7
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

diseño en pseudo-aceleración que deberá ser suavizado de acuerdo a la ​Norma Sismo Resistente
Colombiana del 2010 (NSR-10). Para tener un punto de referencia y poder analizar el espectro de
diseño se elaborará el espectro de diseño según la NSR-10 y la MZSC para la zona de estudio.

- ESPECTRO DE DISEÑO SUAVIZADO O ARMONIZADO.

Obtenidas las respuestas en superficie de cada uno de los sismos se dispone a hacer uso de un método
numérico por medio de un algoritmo en la herramienta MatLab. El método numérico con el cual se
trabajo es ​Newmark. Aceleración Promedio tiene el propósito de calcular la respuesta de un sistema
de 1GDL ante una excitación en su base teniendo como condiciones iniciales propiedades físicas y
dinámicas del sistema de interés. Adicional a esto, y con el propósito de elaborar el espectro de
respuesta para cada sismo, se realiza una rutina para obtener la respuesta del sistema en
Pseudo-Aceleración​ para cada evento sísmico gráfica # 3

Gráfica # 3. Espectro de respuesta de los 9 (nueve) sismos en MatLab

Fuente: propia

Seguidamente se hace un promedio de los 9 (nueve) espectros de respuesta de ​Pseudo-Aceleración y​ a

este promedio se suma la desviación estándar para cada periodo, y de esta manera obtener un solo
espectro de respuesta, así:

8
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Gráfica # 4. Espectro de respuesta del promedio de los 9 (nueve) sismos más su desviación estándar en MatLab
Fuente: propia

Después de obtener nuestro espectro de respuesta promedio más desviación estándar se realiza el
espectro de diseño identificando la siguiente información del espectro de respuesta. Cabe resaltar que
la información extraída de la gráfica 4 se realizó bajo el criterio de ajuste visual del grupo.

Gráfica # 5. Espectro de respuesta del promedio más su desviación estándar en MatLab. Selección de Tc y Sa
Fuente: propia

9
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Los valores marcados en la gráfica # corresponden a Sa=1.761 g y Tc= 0.65 s. Empleando las
siguientes ecuaciones del espectro de diseño de la NSR-10

Tc = 0.48 Av F v
Aa F a ​
(2) S a = 2.5 Aa F a I ​(Parte recta del espectro Elástico de diseño NSR-10)​ ​ ​(3)
Donde I=1. A solicitud del Docente.

De la ecuación (3) se despeja Fa para obtener el valor que le corresponde y una vez obtenido este
valor en la ecuación (2) se despeja Fv y se reemplaza el Valor de Fa para encontrar el Valor de Fv.
Los resultados se Muestran a continuación.

Fa Fv

2.817 3.815
Tabla # 9. Parámetros espectrales (Fa y Fv)
Obtenidos los valores de Tc, Sa, Fa y Fv se realiza una rutina en MatLab para graficar el Espectro de
Diseño obteniendo el siguiente espectro de diseño armonizado o suavizado.

Gráfica # 6. Espectro de respuesta del promedio más su desviación estándar y Espectro de diseño suavizado en MatLab.
Fuente: propia
- ESPECTRO DE DISEÑO NSR-10

En el capítulo A de la NSR-10 se establecen los parámetros mínimos para tener en cuenta en un

diseño estructural. En Colombia se ha clasificado a las diferentes regiones en cuanto a su amenaza
sísmica. Esta amenaza sísmica está definida por dos parámetros, la aceleración pico efectiva (Aa) y la
velocidad pico efectiva (Av) teniendo para Cali los siguientes valores:

10
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Tabla # 10 . Valor de Aa y Av para capitales de departamento

Fuente: Título A. NSR-10.
A continuación se obtiene la velocidad media de la onda cortante:

n
∑ di
i=1
Vs= n
di
; V s = 210.046 m/s (1)
∑ V si
i=1

11
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Tabla # 11 . Clasificación de los perfiles de suelo

Fuente: Título A. NSR-10 (Tabla A.2.4-1
Según la ​Tabla # 11 (Tabla A.2.4-1 NSR-10) s​ e puede identificar al perfil de suelo de la zona
estudiada como tipo D. A continuación se deben encontrar los factores Fa y Fv que son los que tienen
en cuenta la alteración en la amplificación de la zona (efectos de sitio).

Tabla # 12. Valores del coeficiente Fa, para la zona de períodos cortos del espectro
Fuente: Título A. NSR-10 (Tabla A.2.4-3)

12
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Tabla # 13. Valores del coeficiente Fv, para la zona de períodos intermedios del espectro
Fuente: Título A. NSR-10 (Tabla A.2.4-4)

Realizando interpolación lineal para encontrar los valores de Fa y Fv para un tipo de perfil D se tienen
los siguientes valores:

Fa Fv
1.3 1.9

Tabla # 14. Parámetros espectrales (Fa y Fv)

Con todos estos parámetros se puede construir el espectro de diseño basándose en las ecuaciones de la
siguiente gráfica.

Gráfica # 7. Representación esquemática de espectro elástico de aceleraciones de diseño

13
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Fuente:Titulo A.NSR-10 (Figura A.2.6-1)

Gráfica # 8. Espectro de Diseño (NSR-10)

Fuente:propia

- ESPECTRO DE DISEÑO MZSC

Para la elaboración del espectro de diseño de la microzona 6 de la ciudad de Santiago de Cali hay que
tomar como punto de partida la representación esquemática para la construcción de un espectro de
diseño como se muestra a continuación.

Gráfica # 9. Representación esquemática para la construcción del espectro de diseño

Fuente:Informe 5.2.MZSC (Figura 7.12.)

14
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

​Tabla # 15. Parámetros de control para la construcción de los espectros de diseño establecidos.
Fuente: Informe 5.2.MZSC (Tabla 7.2.)

Gráfica # 10. Espectro de diseño Para la Microzonificación Sísmica de Cali

Fuente:Propia

A continuación se presenta los tres (3) espectros de diseño.

15
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

Gráfica # 11. Espectro de Respuesta y Espectros de diseño: Armonizado, MZSC, NSR-10

Fuente:Propia

ESPECTRO ARMONIZADO NSR-10 MZSC

Sa 1,76 0,81 0,68
To 0,14 0,15 0,22
Tc 0,65 0,70 1,15
T​L 9,15 4,56 2,50
​Tabla # 16. Parámetros de los espectros de diseño
Fuente: MZSC , NSR-10

CONCLUSIONES

1. La construcción del espectro de diseño recomendada por la NSR-10 establece una

esquematización que puede ser muy general a la hora de analizar los diferentes perfiles
estratigráficos de diferentes zonas del país es por esto que la esquematización recomendada
por la microzonificación sísmica debería ser más precisa puesto que cuenta con un estudio
detallado de las diferentes zonas de una ciudad en este caso Santiago de Cali.
2. Comparando los espectros de diseño, las Pseudo-aceleraciones MZSC son menores que las
Pseudo-aceleraciones de la NSR-10 y la Armonizada. Esto se debe a que cali tiene unas
Pseudo-aceleraciones máximas establecidas en la MZSC mientras que las
Pseudo-aceleraciones de la armonizada son grandes, esto se debe a que los ajustes son
visuales y no tienen un fundamento matemático propio, ya que este espectro se construye con
las ecuaciones del espectro de diseño de la NSR-10.

16
 UNIVERSIDAD DEL VALLE - FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL Y GEOMÁTICA
PROGRAMA ACADEMICO DE INGENIERIA CIVIL
INGENIERÍA SÍSMICA
ENTREGA # 2
__________________________________________________________________________________

3. Los estudios de las microzonas sísmicas deben realizarse en todas las ciudades del país, ya
que esto permite obtener información más detallada y localizada para aumentar la seguridad
en el diseño de las estructuras y reducción de costos, a su vez, protegiendo a la poblaciòn que
las frecuentan, reduciendo la vulnerabilidad y previniendo posibles desastres.

REFERENCIAS

[1] Notas de clase. Alejandro Cruz

[2] Servicio Geológico Colombiano. Consulta catálogo de Aceleraciones.
https://www2.sgc.gov.co/sismos/sismos/ultimos-sismos.html
[3] Centro Sismológico Nacional-Universidad de Chile. Tipos de sismo en Chile.
https://www.csn.uchile.cl/tipos-de-sismos-chile/
[4] Unidad de Instrumentación Sísmica del Instituto de Ingeniería de la UNAM.
https://aplicaciones.iingen.unam.mx
[5] Normas Colombianas de Diseño y Construcción Sismo-Resistente, NSR-10. Título A

17