1.

Cálculo de los modos de vibración

En el problema se tiene una estructura de 4 pisos, por lo que se tendrá una matriz de rigidez y
masa de 4x4, para el cálculo de la masa m se iterará el comando eig hasta conseguir el periodo
de 0.5s.

Matriz de masas del sistema

[ ]
1 0 0 0
0 0.9 0 0
M= ∗m
0 0 0.9 0
0 0 0 0.7

Calculo de la Matriz de Rigidez

[ ]
2 −1 0 0
−1 1.8 −0.8 0
K= . KO
0 −0.8 1.6 −0.8
0 0 −0.8 0.8

Calibrando para Ti = 0.5 s en MATLAB

[ ]
230704981.770069 −115352490.885035 0 0
−115352490.885035 207634483.593063 −92281992.7080279 0
Kr=
0 −92281992.7080279 184563985.416056 −92281992.7080279
0 0 −92281992.7080279 92281992.7080279

Luego, se planteó el siguiente sistema de ecuaciones que proviene de la vibración libre no

amortiguada del sistema de 4GDL:
2
𝑲∅𝒋 = w j 𝑴∅𝒋

Así, se planteó la búsqueda de los autovalores y auto vectores. Con lo que se obtuvo el
siguiente set de periodos de vibración y sus correspondientes formas de vibración:
 Figura 1. Modos de vibración del sistema de 4GDL

Se tiene que el modo 1 tiene 0 puntos de inflexión, el modo 2 tiene 1 punto de inflexión, el
modo 3 tiene 2 puntos de inflexión y el modo 4 tiene 3 puntos de inflexión

CÁLCULO DE LA RESPUESTA EXACTA: CONDICIÓN INICIAL Q0

El problema implica resolver el siguiente sistema de ecuaciones de movimiento:

𝑴𝒒̈(𝒕) + 𝑪𝒒̇(𝒕) + 𝑲𝒒(𝒕) = 𝟎

Además, con las siguientes condiciones iniciales:

𝒒𝟎 = [1 2 3 4] 𝑇 …. (𝑐𝑚)

𝒒̇𝟎 = [0 0 0 0] 𝑇 …. (𝑐𝑚/𝑠)

GRAFICA DE LOS DESPLAZAMIENTOS NODALES

Figura 2: Desplazamiento en el nodo 1

Figura 3: Desplazamiento en el nodo 2

Figura 4: Desplazamiento en el nodo 3

Figura 5: Desplazamiento en el nodo 4

  El periodo de cada coordenada modal corresponde al modo de vibración
correspondiente, así se tiene que el periodo de la coordenada modal 1 (0.5 seg) es mas
largo que el de la coordenada modal 2,3 y 4.

 El decaimiento depende de la cantidad de ciclos que tiene el sistema, los modos mas
altos presentan mayor cantidad de ciclos en el mismo intervalo de tiempo que los
modos inferiores, debido a ello presentan un mayor decaimiento en el mismo
intervalo.

DEFORMACIONES DE ENTREPISOS

Figura 6: Deformaciones de entrepisos

 Se observa que en el instante inicial (t=0) la deformación relativa de cada piso es igual
a 1cm, lo cual concuerda con la condición inicial dada.

Las fuerzas elásticas sobre los GDL en cada instante de tiempo se calcularon con la siguiente
expresión:

𝒇(𝒕) = 𝑲𝒒(𝒕)

Luego, la cortante basal en el tiempo se calculó como la suma de las fuerzas elásticas sobre los
GDL en cada instante de tiempo. Asimismo, el momento volcante en el tiempo se calculó como
la suma del producto de las fuerzas elásticas sobre los GDL en cada instante de tiempo y su
correspondiente altura desde la base.

FUERZA CORTANTE BASAL

 Figura 7: Corte basal

MOMENTO VOLCANTE

Figura 8: Momento Volcante

 La cortante basal como el momento volcante decaen en su amplitud debido al

amortiguamiento de 2% definido para el sistema de 4GDL.

 Asimismo, se debe notar que en el instante inicial(t=0) existe un valor de cortante

basal y de momento volcante diferente a 0, sucede debido a las condiciones iniciales
de desplazamiento los cuales generan fuerzas elásticas en un sistema de 4GDL.

CÁLCULO DE LA RESPUESTA EXACTA: EXCITACIÓN ARMÓNICA

El problema implica resolver el siguiente sistema de ecuaciones de movimiento:

𝑴𝒒̈(𝒕) + 𝑪𝒒̇(𝒕) + 𝑲𝒒(𝒕) = 𝑷(𝒕)

Además, con las siguientes condiciones iniciales:

 𝒒𝟎 = [0 0 0 0] 𝑇 …. (𝑐𝑚)

𝒒̇𝟎 = [0 0 0 0] 𝑇 …. (𝑐𝑚/𝑠)

Se utilizó el análisis modal, es decir, se desacoplaron las ecuaciones de movimiento

aprovechando la ortogonalidad de los modos.

𝑴̂ 𝒏̈(𝒕) + 𝑪̂𝒏̇ (𝒕) + 𝑲̂𝒏(𝒕) = 𝑷̂(𝒕)

Donde el sistema de ecuaciones anteriores no necesita ser resueltos de manera simultánea,

sino que se puede resolver cada ecuación de manera independiente para cada modo:

𝑚𝑗 ∗𝑛̈𝑗(𝑡) + 𝑐𝑗 ∗𝑛̇𝑗(𝑡) + 𝑘𝑗 ∗𝑛𝑗(𝑡) = 𝑃𝑗(𝑡) 𝑗 = 1: 4

Caso ws=4π rad/s

La excitación armónica es una aceleración en la base del sistema de 4GDL de la forma:

𝑎(𝑡) =0.2𝑔 sin (4𝜋𝑡)

GRAFICA DE LOS DESPLAZAMIENTOS NODALES

Figura 9: Desplazamiento en el nodo 1

Figura 10: Desplazamiento en el nodo 2

 Figura 11: Desplazamiento en el nodo 3

Figura 12: Desplazamiento en el nodo 4

 Se observa que el periodo de cada coordenada modal no se relaciona con el de modo

de vibración correspondiente, sino que se relaciona con el periodo de la excitación
armónica.

 Se observa claramente que en los primeros instantes está presente tanto la respuesta
transitoria como la respuesta estacionaria. Después de esos primeros instantes solo
prevalece la respuesta estacionaria.

Las fuerzas elásticas sobre los GDL en cada instante de tiempo se calcularon con la siguiente
expresión:

𝒇(𝒕) = 𝑲𝒒(𝒕)

Luego, la cortante basal en el tiempo se calculó como la suma de las fuerzas elásticas sobre los
GDL en cada instante de tiempo. Asimismo, el momento volcante en el tiempo se calculó como
la suma del producto de las fuerzas elásticas sobre los GDL en cada instante de tiempo y su
correspondiente altura desde la base.

FUERZA CORTANTE BASAL

Figura 13: Corte basal

MOMENTO VOLCANTE

Figura 14: Momento Volcante

 Tanto la cortante basal como el momento volcante estabilizan su respuesta en el

instante t=10 seg aproximadamente. Esta estabilización lenta se debe a que el
amortiguamiento modal es 2%, par amortiguamientos mayores la estabilización se
encuentra mucho más rápido.

Caso ws=10π rad/s

La excitación armónica es una aceleración en la base del sistema de 4GDL de la forma:

𝑎(𝑡) =0.2𝑔 sin (10𝜋𝑡)

GRAFICA DE LOS DESPLAZAMIENTOS NODALES

Figura 15: Desplazamiento en el nodo 1

Figura 16: Desplazamiento en el nodo 2

Figura 17: Desplazamiento en el nodo 3

 Figura 18: Desplazamiento en el nodo 4

 Se observa que el periodo de cada coordenada modal no se relaciona con el de modo

de vibración correspondiente, sino que se relaciona con el periodo de la excitación
armónica.

 Se observa claramente que en los primeros instantes está presente tanto la respuesta
transitoria como la respuesta estacionaria. Después de esos primeros instantes solo
prevalece la respuesta estacionaria.

Las fuerzas elásticas sobre los GDL en cada instante de tiempo se calcularon con la siguiente
expresión:

𝒇(𝒕) = 𝑲𝒒(𝒕)

Luego, la cortante basal en el tiempo se calculó como la suma de las fuerzas elásticas sobre los
GDL en cada instante de tiempo. Asimismo, el momento volcante en el tiempo se calculó como
la suma del producto de las fuerzas elásticas sobre los GDL en cada instante de tiempo y su
correspondiente altura desde la base.

FUERZA CORTANTE BASAL

Figura 19: Corte basal

MOMENTO VOLCANTE

Figura 20: Momento Volcante

 Tanto la cortante basal como el momento volcante estabilizan su respuesta en el

instante t=10 seg aproximadamente. Esta estabilización lenta se debe a que el
amortiguamiento modal es 2%, par amortiguamientos mayores la estabilización se
encuentra mucho más rápido.

Problema 2
En el problema se tiene una estructura de 1 pisos, por lo que se tendrá una matriz de rigidez y
masa de 3x3, para el cálculo de la masa m se iterará el comando eig hasta conseguir el periodo
de 2s.

Matriz de masas del sistema

[ ]
1 0 0
M = 0 1 0 ∗m
0 0 36

Cálculo de la Matriz de Rigidez

[ ]
3.5 0 2
K= 0 3 4 . K
2 4 34

Calibrando para Ti = 0.5 s en MATLAB

[ ]
35677879.7145770 0 20387359.8369011
K= 0 30581039.7553517 40774719.6738022
20387359.8369011 40774719.6738022 346585117.227319

CÁLCULO DE LA RESPUESTA EXACTA: CONDICIÓN INICIAL

El problema implica resolver el siguiente sistema de ecuaciones de movimiento:

𝑴𝒒̈(𝒕) + 𝑪𝒒̇(𝒕) + 𝑲𝒒(𝒕) = 𝟎

Además, con las siguientes condiciones iniciales:

 𝒒𝟎 = [0 0 0] 𝑇 …. (𝑐𝑚 𝑦 𝑟𝑎𝑑)

𝒒̇𝟎 = [50 20 0] 𝑇 …. (𝑐𝑚/𝑠 𝑦 𝑟𝑎𝑑/𝑠)

Se utilizó el análisis modal, es decir, se desacoplaron las ecuaciones de movimiento

aprovechando la ortogonalidad de los modos.

𝑴̂ 𝒏̈(𝒕) + 𝑪̂𝒏̇(𝒕) + 𝑲̂𝒏(𝒕) = 𝟎

Donde el sistema de ecuaciones anteriores no necesita ser resueltos de manera simultánea,

sino que se puede resolver cada ecuación de manera independiente para cada modo:

𝑚𝑗 ∗𝑛̈𝑗(𝑡) + 𝑐𝑗 ∗𝑛̇𝑗(𝑡) + 𝑘𝑗 ∗𝑛𝑗(𝑡) = 0 𝑗 = 1: 3

GRAFICA DE LOS DESPLAZAMIENTOS NODALES

Figura 21: Desplazamiento en el nodo 1

Figura 22: Desplazamiento en el nodo 2

  El periodo de cada coordenada modal corresponde al modo de vibración
correspondiente, así se tiene que el periodo de la coordenada modal 1 (0.5 seg) es más
largo que el de la coordenada modal 2,3 y 4.

 El decaimiento depende de la cantidad de ciclos que tiene el sistema, los modos más
altos presentan mayor cantidad de ciclos en el mismo intervalo de tiempo que los
modos inferiores, debido a ello presentan un mayor decaimiento en el mismo
intervalo.

Las fuerzas elásticas sobre los GDL en cada instante de tiempo se calcularon con la siguiente
expresión:

𝒇(𝒕) = 𝑲𝒒(𝒕)

Luego, la cortante basal en X, Y y el momento torsor son directamente 𝒇(𝒕) en cada grado de
libertad correspondiente.

FUERZA CORTANTE BASAL

Figura 23: Cortante basal eje X

 Figura 24: Cortante basal eje Y

ANALISIS ARMONICO DEL SISTEMA

El problema implica resolver el siguiente sistema de ecuaciones de movimiento:

𝑴𝒒̈(𝒕) + 𝑪𝒒̇(𝒕) + 𝑲𝒒(𝒕) = 𝑷(𝒕)

Además, con las siguientes condiciones iniciales:

𝒒𝟎 = [0 0 0 0] 𝑇 …. (𝑐𝑚)

𝒒̇𝟎 = [0 0 0 0] 𝑇 …. (𝑐𝑚/𝑠)

Se utilizó el análisis modal, es decir, se desacoplaron las ecuaciones de movimiento

aprovechando la ortogonalidad de los modos.

𝑴̂ 𝒏̈(𝒕) + 𝑪̂𝒏̇ (𝒕) + 𝑲̂𝒏(𝒕) = 𝑷̂(𝒕)

Donde el sistema de ecuaciones anteriores no necesita ser resueltos de manera simultánea,

sino que se puede resolver cada ecuación de manera independiente para cada modo:

𝑚𝑗 ∗𝑛̈𝑗(𝑡) + 𝑐𝑗 ∗𝑛̇𝑗(𝑡) + 𝑘𝑗 ∗𝑛𝑗(𝑡) = 𝑃𝑗(𝑡) 𝑗 = 1: 3

Caso ws=π rad/s

La excitación armónica es una aceleración en la base del sistema de 3GDL de la forma:

𝑎(𝑡) =0.2𝑔 sin (𝜋𝑡)

GRAFICA DE LOS DESPLAZAMIENTOS NODALES

Figura 25: Desplazamiento en el nodo 1

 Figura 26: Desplazamiento en el nodo 2

Figura 27: Desplazamiento en el nodo 3

 Se observa que el periodo de cada coordenada modal no se relaciona con el de modo

de vibración correspondiente, sino que se relaciona con el periodo de la excitación
armónica.

 Se observa claramente que en los primeros instantes está presente tanto la respuesta
transitoria como la respuesta estacionaria. Después de esos primeros instantes solo
prevalece la respuesta estacionaria.

FUERZA CORTANTE BASAL

 Figura 28: Cortante Basal Eje X

Figura 29: Cortante Basal Eje Y

 Tanto la cortante basal como el momento volcante estabilizan su respuesta en el

instante t=3 seg aproximadamente.

Caso ws=1.5π rad/s

La excitación armónica es una aceleración en la base del sistema de 3GDL de la forma:

𝑎(𝑡) =0.2𝑔 sin (1.5𝜋𝑡)

Lo interesante de esta excitación es que su frecuencia circular de vibración coincide

exactamente con la frecuencia circular del primer modo de vibración.

La excitación armónica a la que se somete al sistema de 3GDL es una aceleración en su base

a(t). Esto es equivalente tener al sistema fijo en su base y sometido a un sistema de cargas
equivalentes en cada grado de libertad. Así, en el sistema de ecuaciones original P(t) tiene la
siguiente forma:

𝑷(𝒕) = −𝑴 ∙ 1t ∙ 𝒂(𝒕)

GRAFICA DE LOS DESPLAZAMIENTOS NODALES

Figura 30: Desplazamiento en el nodo 1

Figura 31: Desplazamiento en el nodo 2

 Figura 32: Desplazamiento en el nodo 3

 Se observa que el periodo de cada coordenada modal no se relaciona con el de modo

de vibración correspondiente, sino que se relaciona con el periodo de la excitación
armónica.

 Se observa claramente que en los primeros instantes está presente tanto la respuesta
transitoria como la respuesta estacionaria. Después de esos primeros instantes solo
prevalece la respuesta estacionaria.

FUERZA CORTANTE BASAL

Figura 33: Cortante Basal Eje X

 Figura 34: Cortante Basal Eje Y

 Tanto la cortante basal como el momento volcante estabilizan su respuesta en el

instante t=3 seg aproximadamente.