UNIVERSIDAD PRIVADA
DOMINGO SAVIO
INGENIERÍA CIVIL
PROYECTO FINAL
“TERREMOTOS QUE SACUDEN
EDIFICIOS DE VARIOS PISOS”
ESTUDIANTE: CALLE TIRADO ELIAS
PARY MAMANI PEDRO
MATERIA: ECUACIONES DIFERENCIALES
DOCENTE: ING. LUIS ALBERTO SANCHEZ S.
SEMESTRE: QUINTO SEMESTRE
TARIJA, 27 de Julio de 2019
1.- INTRODUCCIÓN
Los terremotos o sismos así denominados a los eventos telúricos terrestres de
movimiento entre el choque de las placas tectónicas, se caracterizan por la rapidez
con que se generan, el ruido que generalmente lo acompaña, a fin de que la
población los clasificó como el fenómeno natural más trágico y temible, debido
principalmente a que ocurren en una forma repentina e inesperada y por su
capacidad de destrucción. Para comprender este fenómeno es necesario estudiar su
origen, y como se los mide. Además, debemos conocer conceptos que nos ayudaran
a comprender el modelo que analizaremos para un edificio de n pisos.
2.- OBJETIVOS
2.1.- OBJETIVO GENERAL
El propósito principal del siguiente proyecto de Ecuaciones Diferenciales
“Terremotos que Sacuden Edificios de Varios Pisos” es de comprender
matemáticamente como ocurren dichos terremotos
2.2.- OBJETIVOS ESPECÍFICOS
Conocer las ecuaciones diferenciales que nos ayuden a entender dichos
terremotos
Poder entender como son sacudidos los varios pisos de los diferentes
edificios.
Determinar las ecuaciones para estos terremotos.
3.- MARCO TEORICO
3.1.- ¿QUÉ ES UN TERREMOTO?
Un terremoto o sismo es una liberación de energía manifestado por la perturbación
(vibración) del terreno, que se produce por los desplazamientos repentinos a lo
largo de los bordes de las placas tectónicas o por el movimiento de las fallas
geológicas. La zona donde se inicia la liberación de energía se llama foco y su
proyección sobre la superficie se conoce como epicentro.
1
Se considera que los sismos con intensidad mayor o igual a 7.0 son terremotos.
3.2.- ¿CÓMO SE MIDEN LOS SISMOS?
Actualmente se utilizan dos escalas de magnitud, la de Richter y la de Mercalli. La
de Richter mide la energía liberada durante un sismo, mediante el uso de
sismógrafos y la escala de intensidad, por otro lado, la de Mercalli mide los efectos
del sismo.
Para entender el modelo matemático, analizaremos los siguientes conceptos:
3.3.- LEY DE HOOKE
Establece en el m.a.s que la fuerza de restitución, F, opuesta a la dirección del
alargamiento es proporcional a la cantidad de alargamiento x. En concreto F= -kx,
donde k es una constante de proporcionalidad.
3.4.- FRECUENCIA NATURAL
Es la frecuencia a la que un sistema mecánico seguirá vibrando después que se
quita la señal de excitación
3.5.- AMORTIGUAMIENTO DE COULOMB
Constante mecánica de amortiguación en la que la energía es absorbida por la
fricción por deslizamiento. Es un valor constante sin importar el deslizamiento o la
velocidad.
3.6.- RESONANCIA
Es un fenómeno que se produce cuando un cuerpo capaz de vibrar es sometido a
la acción de una fuerza periódica, cuyo periodo de vibración coincide con el
periodo de vibración característico de dicho cuerpo. En el cual una fuerza
relativamente pequeña aplicada en forma repentina, hace que una amplitud de un
sistema oscilante se haga muy grande.
2
4.- ANÁLISIS DEL MODELO MATEMÁTICO
Por lo general los grandes terremotos tienen un efecto devastador en los edificios por
ejemplo el famoso terremoto de 1906 en San Francisco destruyo gran parte de esa
ciudad, más recientemente, el área fue golpeada una segunda vez por el terremoto
de loma prieta.
En este proyecto moldearemos el efecto de un terremoto en un edificio de varios
pisos y luego resolvemos e interpretamos las matemáticas, sea que x, represente el
desplazamiento horizontal del i-esimo piso a partir de su posición de equilibrio, aquí
la posición de equilibrio será un punto fijo en el suelo de forma que x°= 0, durante el
terremoto la tierra se mueve horizontalmente, de modo que se considera el
desplazamiento de cada piso con respecto al suelo, supongamos que el i-esimo piso
tiene una masa m, y que los pisos sucesivos están conectados por un conector
elástico cuyo efecto se maneja al de un resorte, por lo general los elementos
estructurales en edificios grandes están hechos de acero , un material altamente
elástico, de manera que cada conector proporcione una fuerza de restauración
cuando los pisos se encuentran desplazados con respecto a los demás, suponemos
que la ley de Hooke es válida, que constituye proporcionalidad K, entre el i-esimo y el
(i+1)esimo piso, es decir la fuerza restauradora entre estos dos pisos es.
F = K¹(X(i+1) – Xi)
Donde (X(i+1)-Xi) es el desplazamiento (cambio) de (i+1)esimo piso respecto al i-
esimo. También suponemos una relación similar entre el primer piso y el suelo, con
constante de proporcionalidad Kg.
mn Kn-1
mn-1 Kn-2
: :
mi+1
m2
K1
mi
m1
K2 ki-1(xi-xi-1) mi-1
piso
Pisos del edificio Fuerzas en el i-ésimo piso
3
Se puede aplicar la segunda ley de movimiento de newton F=m.a
∂2 x
m1 =−k g x 1+ k 1 ( x 2−x1 )
∂ t2
∂2 x
m2 =−k 1 ( x 2−x1 )+k 1 (x 3−x 2 )
∂t 2
∂2 x
m n 2 =−k n−1 (x n−x n−1)
∂t
Considerando un edificio de dos pisos cada piso con una masa m = 5000kg y cada
fuerza de restauración constante con un valor k = 10000 kg/s² entonces las
ecuaciones diferenciales son;
∂2 x
=−4 x 1+2 x 2
∂ t2
∂2 x 2
=2 x1−2 x 2
∂ t2
La solución por los métodos de la sección es:
x 1 ( t )=2 c 1 cos ∞ 1 t+ 2 c2 sen ∞2 t+ 2c 3 cos ∞ 3 t+2 c 4 sen ∞ 4 t
x 2 ( t )=(4−∞ ²) c 1 cos ∞1 t+(4−∞ ²) c 2 sen ∞ 2 t +(4−∞ ²)c 3 cos ∞ 3 t+(4−∞ ²)c 4 sen ∞ 4 t
Donde ∞ 1= √3+ √ 5 = 2.288 y ∞ 2= √ 3−√ 5 = 0.874 ahora se supone que se aplican las
siguientes condiciones iniciales x 1 ( 0 )=0 x ' 1 ( 0 )=0 x 2 ( 0 )=0 x ' 2 ( 0 )=0 estas
corresponden a un edificio en la posición de equilibrio con el primer piso con una
rapidez horizontal dada de 0.2 m/s. La solución del problema con valores iniciales es:
x 1 ( t )=2 c 2 sen ∞ 2 t +2 c 4 sen ∞ 4 t
4
x 2 ( t )=(4−∞ ²) c 2 sen ∞ 2 t +(4−∞ ²)c 4 sen ∞ 4 t
2 2 2
donde c 2=( 4−∞ 2 ) 0.1/ [ (∞ 1−∞ 2 )∞ 1 ] = 0.0317 = c 4 = vea en las figuras 3 y 4 para las
graficas de
x 1 ( t ) y x 2 (t) advierta que x1 inicialmente se desplaza hacia la derecha pero es frenado
por el arrastre de x 2 esta inicialmente en reposo pero se acelera debido al jalo ¿n de
x1, pasando a x1 en menos de un segundo, continua hacia la derecha jalando hasta x 1
hasta la marca de dos segundos, en ese momento el arrastre de x 1 ha ralentizado a
x2 hasta detenerlo , después x2 se mueve hasta la izquierda pasando por el punto de
equilibrio en 3.2 segundos y continua moviéndose hacia la izquierda arrastrando a x 1
junto con el movimiento hacia atrás y adelante continua no hay ningún
amortiguamiento en el sistema, por lo que el comportamiento oscilatorio continua por
5
siempre.
FIGURA 3 Grafica de x1(t)
0.4
0.3
0.2
0.1
0
0 1 2 3 4 5 6
-0.1
-0.2
-0.3
-0.4
6
FIGURA 4 Grafica x2(t)
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
0 1 2 3 4 5 6
-0.05
-0.1
-0.15
Se aplica una fuerza horizontal oscilatoria de frecuencia ω1 o ω2, se analiza y se
espera que se produzca grandes oscilaciones del edificio, posiblemente causando un
gran daño si el terremoto dura una cantidad considerable de tiempo.
Se define las siguientes matrices y vectores:
x 1 (t )
[ ]
x ( t )= x 2 ( t )
⋮
xn(t )
7
Entonces el sistema de ecuaciones se puede escribir en forma matricial
∂2 x
M =KX O MXʺ =KX
∂ t2
Advierta que la matriz M es una matriz diagonal con la masa del i-ésimo piso
elemento diagonal. La matriz M tiene una inversa dada por
m1−1 0 0 … 0
−1
M = 0
⋮ (
m2
¿
−1 −1
0 … 0 0 ¿ … ¿ mn ¿
¿ 0 ¿ )
Por lo tanto, se puede representar con la ecuación:
X ˚= ( M −1 K ) X o Xʺ=AX
Donde A=M −1 K , la matriz M se denomina la matriz de masa, y la matriz K es la
matriz de rigidez.
Los valores propios de la matriz A revelan la estabilidad del edificio durante un
terremoto. Los valores propios de A son negativos y distintos. En el primer ejemplo,
los valore propios son −3+ √5=−0.764 y −3−√5=−5.236 . Las frecuencias naturales
del edificio son las raíces cuadradas de los negativos de los valores propios. Si λ, es
el i-ésimo valor propio, entonces ω 1=√−λ, es la i-ésima frecuencia i=1,2,...., n.
Durante un terremoto se aplica una gran fuerza horizontal al primer piso. Si ésta es
de naturaleza oscilatoria, digamos, de la forma F(t)=G cos γt, entonces se pueden
desarrollar grandes desplazamientos en el edificio, especialmente si la frecuencia del
término de forzamiento es cercana a una de las frecuencias naturales del edificio.
Esto es una reminiscencia del fenómeno de resonancia.
Analicemos otro ejemplo: supongamos que tenemos un edificio de 10 pisos, donde
cada piso tiene una masa 10000kg y cada valor k 1 es de 5000kg/s2. Entonces
8
−1 0.5 0 0 0 0 0 0 0 0
0.5 −1 0.5 0 0 0 0 0 0 0
( )
0 0.5 −1 0.5 0 0 0 0 0 0
0 0 0.5 −1 0.5 0 0 0 0 0
A=M −1 K = 0 0 0 0.5 −1 0.5 0 0 0 0
0 0 0 0 0.5 −1 0.5 0 0 0
0 0 0 0 0 0.5 −1 0.5 0 0
0 0 0 0 0 0 0.5 −1 0.5 0
0 0 0 0 0 0 0 0.5 −1 0.5
0 0 0 0 0 0 0 0 0.5 −0.5
Los valores propios de A se pueden hallar fácilmente utilizando Mathematica u otro
paquete similar. Estos valores son -1.956, -1.826, -1.623, -1.365, -1.075, -0.777, -0.5,
-0.267, -0.099 y -0.011, con las frecuencias correspondientes, 1.399, 1.351, 1.274,
1.168, 1.037, 0.881, 0.707, 0.517, 0.315 y 0.105, y periodos de oscilación (2π/ω)
4.491, 4.651, 4.932, 5.379, 6.059, 7.132, 8.887,12.153, 19.947 y 59.840. Durante un
terremoto típico cuyo periodo podría estar en el rango de 2 a 3 segundos, este
edificio no parece estar en peligro de desarrollar resonancia. Sin embargo, si los
valores k fueran 10 veces más grandes, entonces, por ejemplo, el sexto periodo
podría ser de 2.253 segundos, mientras que el quinto al séptimo son todos del orden
de 2 a 3 segundos. Este edificio sería más propenso a sufrir daños en un terremoto
típico de un periodo de 2 a 3 segundos.
9
5.- CONCLUSIONES
Es importante tener conocimiento en la forma en que se emplean fórmulas de
diferentes ramas de las ciencias a ecuaciones diferenciales, ya que el uso de las
ecuaciones diferenciales a las ciencias, permiten obtener datos precisos.
En este proyecto, se hizo énfasis en el daño que puede ocasionar que la frecuencia
del sismo coincida con la frecuencia característica de un piso del edificio, pero no
solo esto puede ocasionar daños a un edificio, existen más factores, los cuales
también se pueden analizar mediante ecuaciones diferenciales.
10
Contenido
1.- INTRODUCCIÓN......................................................................................................1
2.- OBJETIVOS..............................................................................................................1
2.1.- OBJETIVO GENERAL.......................................................................................1
2.2.- OBJETIVOS ESPECÍFICOS..............................................................................1
3.- MARCO TEORICO...................................................................................................1
3.1.- ¿QUÉ ES UN TERREMOTO?...........................................................................1
3.2.- ¿CÓMO SE MIDEN LOS SISMOS?..................................................................2
3.3.- LEY DE HOOKE.................................................................................................2
3.4.- FRECUENCIA NATURAL..................................................................................2
3.5.- AMORTIGUAMIENTO DE COULOMB..............................................................2
3.6.- RESONANCIA....................................................................................................2
4.- ANÁLISIS DEL MODELO MATEMÁTICO................................................................3
5.- CONCLUSIONES.....................................................................................................9
