1

UNIVERSIDAD NACIONAL DE INGENIERIA

FACULTAD DE INGENIERIA CIVIL

ANÁLISIS DEL MOVIMIENTO OSCILATORIO EN CONSTRUCCIONES ANTISÍSMICAS

DE DOS GRADOS DE LIBERTAD MEDIANTE EL PROGRAMA OCTAVE., LIMA-2020

CURSO DE FÍSICA II. SECCIÓN J

ELABORADO POR

Espinoza Anccasi, Vladimir

Guadalupe Chuquisengo Edwin

Quiroz Centurión, Frank Xavier

Rojas Baldoceda, Carlos Teodoro

Turín Granda, Leslie Katheryn

Zapata Quispe, Edson Aaron

PROFESOR Mag. Lic. Dane Bruce Cachi Eugenio

LIMA - PERÚ 2020

 2

RESUMEN

El objetivo de este proyecto es proponer una herramienta y demostrar su aplicación en el

análisis de sistemas de suelos y estructuras de dos grados de libertad, mediante una rutina de

cálculos en Octave, el cual está enfocado en el procesamiento y análisis de señales de ondas

sísmicas obtenidas en un terreno y posteriormente discutir los resultados obtenidos.

 3

ABSTRAC

The objective of this project is to propose a tool and demonstrate its application in the analysis

of soil systems and structures with two degrees of freedom, through a routine of calculations in

Octave, which is focused on the processing and analysis of seismic wave signals. obtained in a

field and later discuss the results obtained.

 4

ÍNDICE

INTRODUCCIÓN ............................................................................................................ 6

PLANTEAMIENTO DEL PROBLEMA. ...................................................................... 7

OBJETIVO ........................................................................................................................ 7

JUSTIFICACIÓN DE LA INVESTIGACIÓN .............................................................. 7

MARCO TEÓRICO ......................................................................................................... 8

LIMITACIONES E IDEALIZACIONES DE LA METODOLOGÍA ..................... 8

TEORÍA DEL CÁLCULO MATRICIAL .................................................................. 9

PROGRAMACIÓN EN MATLAB. ............................................................................ 9

¿Qué es MATLAB? .................................................................................................... 9

Elementos programables. .......................................................................................... 10

Estructuras condicionales.......................................................................................... 11

Estructuras cíclicas.................................................................................................... 11

Sentencias y comandos propios de MATLAB. ......................................................... 12

DESARROLLO DEL SOFTWARE ............................................................................. 12

ESPECTROS DE RESPUESTA Y MODELOS SISMOLÓGICOS: .................... 12

Espectros de respuesta .............................................................................................. 12

Modelos sismológicos............................................................................................... 13
 5

Ventajas del uso de espectros de respuesta: .............................................................. 13

PROGRAMA INFORMÁTICO: .............................................................................. 14

SOLUCIÓN AL PROBLEMA DE LA DURACIÓN ESTACIONARIA

EQUIVALENTE: .................................................................................................................... 14

DURACIONES UMBRAL Y DE LA ENERGÍA.................................................... 14

PROCEDIMIENTO: ................................................................................................. 15

APLICACIÓN PRÁTICA ............................................................................................ 16

MAGNITUD 5 ............................................................................................................. 17

MAGNITUD 6 ............................................................................................................. 18

MAGNITUD 7 ............................................................................................................. 19

COMPARACIÓN ENTRE MAGNITUDES ............................................................ 20

CONCLUSIONES........................................................................................................... 21
 6

INTRODUCCIÓN

El Perú está entre las regiones de más alta actividad sísmica que existe en la tierra.

Este fenómeno conlleva a varias pérdidas socioeconómicas dependiendo de la intensidad, por ello

es necesario realizar un exhaustivo estudio ante la presencia de un fenómeno sísmico, la estructura

de dos niveles oscila en distintos modos de vibración según el nivel en el que nos encontremos.

Nos preguntamos ¿Cómo podríamos analizar los modos de vibración a consecuencia de las ondas

sísmicas en una estructura?, en esta investigación optamos por la aplicación Matlab que es un

software computacional que te permite realizar tu algoritmo que se ajuste a cualquier necesidad.

Existen principalmente dos tipos de técnicas en el análisis de estructuras que son técnicas de

análisis matricial de estructuras (AME) y análisis por el método de elementos finitos (MEF).

Ambos tipos requieren saber determinados conceptos que muchas veces resulta ser un proceso

engorroso, este proyecto busca agilizar el cálculo de estructuras mediante la aplicación programada

MATLAB, este programa en actualidad es muy usado por ingenieros, investigadores y analistas,

contiene un sistema numérico computarizado con entorno de desarrollo integrado

 7

PLANTEAMIENTO DEL PROBLEMA.

Ante la presencia de un fenómeno sísmico la estructura oscila en distintos modos de vibración

según sea el nivel en el cual nos encontremos. ¿Cómo podríamos analizar los modos de vibración

a consecuencia de las ondas sísmicas en una estructura?

OBJETIVO

Proponer una herramienta y demostrar su aplicación en el análisis de sistemas de suelos y

estructuras de dos grados de libertad, mediante una rutina de cálculos en Octave, el cual está

enfocado en el procesamiento y análisis de señales de ondas sísmicas obtenidas en un terreno y

posteriormente discutir los resultados obtenidos.

JUSTIFICACIÓN DE LA INVESTIGACIÓN

Los cálculos necesarios para determinación de una respuesta sísmica, sin la utilización de

herramientas computacionales, implican una serie de procesos engorrosos que derivan en posibles

errores en los resultados. Existe la necesidad de un programa que ayude a comprender los

diferentes tipos procesos que involucran el cálculo para el estudio del análisis modal y dinámico

de una estructura. El programa obtiene la respuesta dinámica de una estructura modelada como un

oscilador ante la acción del movimiento del suelo. Con este sistema es posible modelar estructuras

de uno hasta 50 grados de libertad mediante un oscilador de tipo “péndulo invertido”. La respuesta

sísmica es obtenida con el método numérico de la integral de Duhamel.

 8

MARCO TEÓRICO

LIMITACIONES E IDEALIZACIONES DE LA METODOLOGÍA

Para el desarrollo de esta metodología de análisis dinámico se tuvieron en cuenta ciertas

consideraciones que serán descritas a continuación (Maldonado & Chio Cho, 2004):

❖ “La masa de la estructura está concentra al nivel de los pisos” (p.119). Es decir que sólo se

trabaja con masas traslacionales. De esta forma, el problema se vuelve uno que tiene tantos

grados de libertad como masas concentradas.

❖ “Las vigas en los pisos no tienen deformaciones axiales debido a la rigidez infinita del

diafragma, y por ello en cada piso existe un solo grado de libertad horizontal” (p.119).

❖ “Se desprecia la deformación axial de las columnas, siempre y cuando la estructura sea poco

esbelta, baja y larga” (p.119). Esto lleva a que se desprecien los grados de libertad verticales por

la cercana nulidad en las deformaciones de los elementos verticales llamados columnas. Esta

condición garantiza que las vigas se mantengan horizontales durante el movimiento de la

estructura.

❖ Otra consideración muy importante es que “Debido a la uniformidad de las propiedades

estructurales a lo largo del edificio, y a que un mismo pórtico se repite […], el análisis de un

pórtico interior puede representar respuesta para todo el edificio” (p.140). en una dirección

(COMPLEMENTAR).

❖ Esta metodología está aplicada a estructuras aporticadas ortogonales de múltiples grados de

libertad y que se encuentran sometidas a una acción sísmica en un sentido.

 9

TEORÍA DEL CÁLCULO MATRICIAL

Generalidades

El método de la rigidez (o también conocido como cálculo matricial) está basado en la

superposición de desplazamientos, donde los desplazamientos en los nudos son considerados

como incógnitas {U} y donde las fuerzas en los extremos de los elementos de la estructura

cumplen las ecuaciones de equilibrio de los respectivos nudos; consiguiendo así un sistema de

ecuaciones lineales {𝑓}=[𝐾]∗{𝑈}.

El método de análisis matricial (AME) permite la construcción y operación de las matrices de

rigidez, vectores de fuerzas externas y vectores de desplazamiento de los elementos en forma

local y global de la estructura, mediante la Ecuación 1. {𝑓}=[𝐾]∗{𝑈}

Ecuación 1

[𝐾]= Matriz de Rigidez de la estructura.

{𝑈}= Vector de Desplazamiento de la estructura.

{𝑓} = Vector de Fuerzas Externas de la estructura.

PROGRAMACIÓN EN MATLAB.

¿Qué es MATLAB?

MATLAB es un paquete de software de computación y visualización, así como un

lenguaje de programación de cuarta generación. Este software permite realizar manipulaciones

matriciales, trazado de funciones, implementación de algoritmos y muchas operaciones

matemáticas de alto nivel. MATLAB te permite crear tu propia herramienta personalizada o

algoritmo que se ajuste a cualquier necesidad.

 10

Elementos programables.

Objetos Gráficos.

Los objetos gráficos son los componentes visuales utilizados por MATLAB® para mostrar

datos gráficamente. Cada objeto tiene un identificador único llamado Manejar. Con este

identificador, se pueden manipular las propiedades de un objeto existente.

Interfaz Gráfica del Usuario.

La interfaz gráfica de usuario, también conocida como GUI o (por las siglas Grafical User

Interface), es la forma en la que el usuario puede acceder e interactuar con las aplicaciones y

programas diseñados en el software MATLAB.

MATLAB dispone de una utilidad para el desarrollo interactivo de interfaces de usuario

(GUI) llamado GUIDE. Las Guide se componen de diferentes elementos tales como botones,

campos editables de texto, menús, sliders (deslizadores), gráficos, etc. El editor permite construir

interfaces arrastrando y soltando componentes en el área de diseño de la GUI.

Variables.

Las variables son elementos que se denotan con palabras claves, y que tienen la capacidad

de generar un espacio de memoria para que los llamamientos hagan uso de dicha reserva para

guardar datos.

Para crear una nueva variable en MATLAB, se introduce el nombre de la variable en la

ventana editor, seguida de un signo igual (=) y el valor que desea asignar a la variable.
 11

FIGURA 1 1. Ejemplo de variables

Estructuras condicionales.

En informática las declaraciones condicionadas son características de un lenguaje de

programación, que permiten realizan diferentes cálculos o acciones dependiendo de si una

condición booleana especificada por el programador, evaluándola como verdadera o falsa.

Se tienen en este grupo la estructura if, la estructura if-else, los operadores lógicos, la

estructura select case y la estructura with.

Estructuras cíclicas

For y while.

Permiten repetir las mismas operaciones o análogas sobre datos distintos.

Ciclo while.

Repite un conjunto de sentencias mientras haya elementos distintos de cero en condición,

es decir, mientras haya algún o algunos elementos “true”. El bucle se termina cuando todos los

elementos de condición son “false” (es decir, cero).

Ciclo for.

Repite un conjunto de sentencias un número predeterminado de veces.

 12

Sentencias y comandos propios de MATLAB.

Los comandos son palabras claves dentro de los lenguajes de programación, que tienen

como función realizar ciertas acciones que son propias del lenguaje usado y que están intrínsecas

en la estructura del mismo. Algunos comandos son:

Comandos find.

Busca índices correspondientes a elementos de vectores que cumplen una determinada

condición. El resultado es un vector con los índices de los elementos que cumplen la condición.

Comandos sum.

Este comando permite sumar de los elementos de un vector o matriz, según los

parámetros de rango y los paramentos de selección definidos para la función sum, es decir, qué

elementos quieren que se sumen.

Comandos text.

Esta función permite agregar descripciones de texto a puntos de datos que están siendo

graficados mediante el comando plot, lo que permite darles etiquetas a los puntos graficados.

DESARROLLO DEL SOFTWARE

ESPECTROS DE RESPUESTA Y MODELOS SISMOLÓGICOS:

Espectros de respuesta

Es la respuesta en aceleraciones de un sistema de un grado de libertad de periodo

propio T a los sismos de cierta región.

 13

FIGURA 1 2

Modelos sismológicos

Están basados en los fenómenos físicos que producen un terremoto y nos dan el

espectro de amplitudes de Fourier de un sismo en función de parámetros sísmicos como

magnitud, distancia epicentral, aceleración.

Ventajas del uso de espectros de respuesta:

❖ Permite un cálculo en el tiempo

❖ Gracias a esta herramienta, se pueden obtener sismogramas sintéticos
❖ Mayor capacidad de cálculo
❖ Relaciona la respuesta de la estructura, con parámetros sismológicos2.
Programa Informático
 14

PROGRAMA INFORMÁTICO:

FIGURA 3

SOLUCIÓN AL PROBLEMA DE LA DURACIÓN ESTACIONARIA

EQUIVALENTE:

DURACIONES UMBRAL Y DE LA ENERGÍA

❖ Duración umbral “Threshold duration” Tg,t es definida como el intervalo de tiempo en el

que el acelerograma pasa por primera vez y por última por un cierto nivel umbral

Duración de la energía de Husid

𝑡
𝐻𝑎𝑔 (𝑡) = ∫0 𝑎𝑔2 (𝑢)𝑑𝑢

En la bibliografía podemos encontrar expresiones que relacionan la duración umbral y la

duración de la energía de Husid, con parámetros sísmicos

 15

❖ Duración estacionaria equivalente

La duración estacionaria equivalente es la más apropiada para nuestro programa puesto que

facilita la obtención de espectros de varianza, que son necesarios para el cálculo de espectro

de respuesta

𝑡𝑔𝑡
2 (𝑡)𝑑𝑡
𝑇𝑔𝑠 = ∫ 𝐼𝑎𝑔,𝑢
0

PROCEDIMIENTO:

❖ Se va a seleccionar un número de sismos para analizar su duración

❖ Para ello, utilizamos la base de datos del European Strong Motion Database
❖ El objetivo es lograr una relación entre la duración estacionaria equivalente y la duración
umbral o la duración de la energía, para así relacionarla con parámetros sísmicos

❖ Criterios de selección de sismos:

❖ Magnitud entre 4 y 9
❖ Distancia al epicentro entre 10 y 200km
❖ Aceleración del terreno, entre 0,06 hasta más de 0,3g
❖
❖
❖ Base de datos: European Strong Motion Database
❖ Número de sismos: 217
 16

APLICACIÓN PRÁTICA

❖ Magnitud: se va a analizar a evolución en 3 casos, magnitud 5, 6, y 7

❖ Tensión efectiva: 𝜎 = 200 𝑏𝑎𝑟
❖ Distancia hipocentral: 25km
❖ Densidad: 𝜌 = 2000 𝑘𝑔/𝑚3
❖ Frecuencia: 2Hz
❖ Tipo de suelo: Capa de tierra sobre roca
❖ Para el cálculo de la función de intensidad se va a utilizar la expresión de Saragoni y Hart
(1974)

𝐼 (𝑡) = 𝑘 + 𝑡 𝑛 + 𝑒 −𝑘+𝑡
 17

MAGNITUD 5

Espectro respuesta M5
0.12

0.1

0.08
Rsa

0.06

0.04

0.02

0
0 3 6 9 12 15 18 21 24 27 30
Tn (s)

FIGURA 4 Espectro respuesta M5

Espectro respuesta unitaria M5

3.5

2.5

2
Rsa,u

1.5

0.5

0
0 3 6 9 12 15 18 21 24 27 30 33
Tn (s)

FIGURA 5 Espectro respuesta unitaria M5

 18

MAGNITUD 6

Espectro respuesta M6
0.4
0.35
0.3
0.25
Rsa

0.2
0.15
0.1
0.05
0
0 3 6 9 12 15 18 21 24 27 30
Tn (s)

FIGURA 1 6Espectro de respuesta M6

Espectro respuesta unitario M6

3.5

2.5

2
Rsa,u

1.5

0.5

0
0 3 6 9 12 15 18 21 24 27 30
Tn (s)

FIGURA 7 1 Espectro respuesta unitario M6

 19

MAGNITUD 7

Espectro respuesta M7
1.2

0.8
Rsa

0.6

0.4

0.2

0
0 3 6 9 12 15 18 21 24 27 30
Tn (s)

FIGURA 8 Espectro respuesta M7

Espectro respuesta unitario M7

3

2.5

2
Rsa,u

1.5

0.5

0
0 3 6 9 12 15 18 21 24 27 30
Tn (s)

FIGURA 9 Espectro respuesta unitario M7

 20

COMPARACIÓN ENTRE MAGNITUDES

1.2 Espectros respuesta

1
0.8
Rsa

0.6 M5
0.4 M6
0.2 M7
0
0 3 6 9 12 15 18 21 24 27 30
Tn (s)
FIGURA 10 Espectro de respuesta

Espectros respuesta unitaria

4
3
Rsa,u

2 M5
1 M6
0 M7
0 3 6 9 12 15 18 21 24 27 30 33
Tn (s)

FIGURA 11 Espectros respuesta unitaria

 21

CONCLUSIONES

❖ Se ha conseguido obtener una relación válida entre duraciones, lo que permite relacionar
la duración estacionaria equivalente con parámetros sísmicos
❖ Gracias a ello, se puede completar el diseño del programa de cálculo de espectros de
respuesta sísmica
❖ Mediante este programa, se puede mejorar la construcción de estructuras ya que el
espectro de respuesta se basa en parámetros sísmicos y no en aproximaciones a otros
países como se realiza en la NCSE-02
❖ Con la intervención de un Informático, se podría llegar a comercializar el programa
 22

BIBLIOGRAFÍA

❖ Hernández Martínez, M. D., & Mercado López, O. J. (2015). Diseño y aplicación de

códigos en Matlab para cálculo de parámetros dinámicos en análisis sísmico de suelo y
estructuras (Doctoral dissertation, Universidad Centroamericana).
❖ Hurtado Gómez, J. E. (2000). Dinámica de estructuras: aplicaciones a la ingeniería
sísmica. Centre Internacional de Mètodes Numèrics en Enginyeria (CIMNE).
❖ Coral, H. A., Sandoval, J. A., Rosero, E. E., Ramírez, J. M., & Gómez, D. Simulador
Sísmico Uniaxial Tele-Operable para Modelos Estructurales de Pequeña Escala (Doctoral
dissertation, Tesis de pregrado, Universidad del Valle, Cali, Colombia).