UNIVERSIDAD CATÓLICA BOLIVIANA “SAN PABLO”

UNIDAD ACADÉMICA SANTA CRUZ

DEPARTAMENTO DE POSTGRADO

TRABAJO FINAL

Módulo: Principios de Dinámica Estructural

Docente: MSc. Ing. Adhemar Flores Valdivia

Participante: Ing. Tomás Wilson Alemán Ramírez

Santa Cruz – Bolivia

 DATOS DE ENTRADA
Número de Pisos: 4
Altura de Pisos: 3.5 m.
Amortiguamiento: 5% del Amortiguamiento Critico.
Columnas HoAo de 30x40cm.
Vigas HoAo: 20x60cm.
Muro Pantalla de 15cm. de espesor.
Losa de HoAo: Espesor Equivalente de 18 cm.
Carga Adicional sobre Losa: 407.70 Kg/m2
Densidad de HoAo: 2,400 Kg/m3
Matriz de Rigidez [N/m] Plano ZY:

Matriz de Rigidez [N/m] Plano ZX:

DETERMINAR:
- Eigenvalores, eigenvectores en ambas Direcciones Principales.
- Periodos de vibración para cada modo en ambas Direcciones Principales.
- Formas de los Modos de Vibración en ambas Direcciones principales (Normalizadas).
A la dirección del Plano ZY, aplicar el Espectro de Respuesta de Aceleración del Sismo de Kobe-
Japón (NGA1111_KobeJapan_NishiAkashi_00).
- Determinar las ecuaciones desacopladas de cada modo de vibración.
- Determinar Factores de Participación Modal.
- Determinar Masa Modal Efectiva (Establecer cuantos modos serían suficientes para tener una
participación de masa efectiva de 90%)
- Determinar los desplazamientos y la Deriva en ejes Globales para cada nivel.
- Determinar las Fuerzas Resultantes en cada Nivel.
OPCIONAL:
Determinar las ecuaciones de movimiento en coordenadas Globales de cada nivel si se tiene una
fuerza aplicada en el piso 4: 𝐹4(𝑡) = 100 ∗ 103 ∗ sin⁡(90 ∗ 𝑡) [N].
 SOLUCIÓN

1ERO: EIGENVALORES, EIGENVECTORES EN AMBAS DIRECCIONES PRINCIPALES.

a) Modelo Físico

u1 u2 u3 u4
K1 K2 K3 K4
m1 m2 m3 m4

b) Ecuaciones

Diagrama de cuerpo libre para m1

∑ 𝐹 = 𝑚1 ∙ 𝑢̈ 1 → ⨁
K1∙u1 K2(u2- u1)
−𝐾1 ∙ 𝑢1 + 𝐾2 (𝑢2 − 𝑢1 ) = 𝑚1 ∙ 𝑢̈ 1
m1
𝑚1 ∙ 𝑢̈ 1 + 𝐾1 ∙ 𝑢1 − 𝐾2 ⋅ 𝑢2 + 𝐾2 ⋅ 𝑢1 = 0

𝑚1 ∙ 𝑢̈ 1 + (𝐾1 + 𝐾2 )𝑢1 − 𝐾2 ⋅ 𝑢2 = 0

Diagrama de cuerpo libre para m2

∑ 𝐹 = 𝑚2 ∙ 𝑢̈ 2 → ⨁

K2(u2- u1) K3(u3- u2) −𝐾2 (𝑢2 − 𝑢1 ) + 𝐾3 (𝑢3 − 𝑢2 ) = 𝑚2 ∙ 𝑢̈ 2

m2
𝑚2 ∙ 𝑢̈ 2 + 𝐾2 (𝑢2 − 𝑢1 ) − 𝐾3 (𝑢3 − 𝑢2 ) = 0

𝑚2 ∙ 𝑢̈ 2 − 𝐾2 ⋅ 𝑢1 + (𝐾2 + 𝐾3 ) ⋅ 𝑢2 − 𝐾3 ⋅ 𝑢3 = 0

Diagrama de cuerpo libre para m3

∑ 𝐹 = 𝑚3 ∙ 𝑢̈ 3 → ⨁
K3(u3- u2) K4(u4- u3)
m3 −𝐾3 (𝑢3 − 𝑢2 ) + 𝐾4 (𝑢4 − 𝑢3 ) = 𝑚3 ∙ 𝑢̈ 3

𝑚3 ∙ 𝑢̈ 3 + 𝐾3 (𝑢3 − 𝑢2 ) − 𝐾4 (𝑢4 − 𝑢3 ) = 0

𝑚3 ∙ 𝑢̈ 3 − 𝐾3 ⋅ 𝑢2 + (𝐾3 + 𝐾4 ) ⋅ 𝑢3 − 𝐾4 ⋅ 𝑢4 = 0
Diagrama de cuerpo libre para m4
∑ 𝐹 = 𝑚4 ∙ 𝑢̈ 4 → ⨁

K4(u4- u3) −𝐾4 (𝑢4 − 𝑢3 ) = 𝑚4 ∙ 𝑢̈ 4

m4
𝑚4 ∙ 𝑢̈ 4 + 𝐾4 ⋅ 𝑢4 − 𝐾4 ⋅ 𝑢3 = 0

𝑚4 ∙ 𝑢̈ 4 − 𝐾4 ⋅ 𝑢3 + 𝐾4 ⋅ 𝑢4 = 0

Resumen de las ecuaciones obtenidas:

𝑚1 ∙ 𝑢̈ 1 + (𝐾1 + 𝐾2 )𝑢1 − 𝐾2 ⋅ 𝑢2 = 0

𝑚2 ∙ 𝑢̈ 2 − 𝐾2 ⋅ 𝑢1 + (𝐾2 + 𝐾3 ) ⋅ 𝑢2 − 𝐾3 ⋅ 𝑢3 = 0

𝑚3 ∙ 𝑢̈ 3 − 𝐾3 ⋅ 𝑢2 + (𝐾3 + 𝐾4 ) ⋅ 𝑢3 − 𝐾4 ⋅ 𝑢4 = 0

𝑚4 ∙ 𝑢̈ 4 − 𝐾4 ⋅ 𝑢3 + 𝐾4 ⋅ 𝑢4 = 0

Aplicando un formato matricial obtenemos:

𝑚1 0 0 0 𝑢̈ 1 𝐾1 + 𝐾2 −𝐾2 0 0 𝑢1
0 𝑚2 0 0 𝑢̈ −𝐾2 (𝐾2 + 𝐾3 ) −𝐾3 0 𝑢2
[ ]⋅{ 2}+[ ] ⋅ {𝑢 } = 0
0 0 𝑚3 0 𝑢̈ 3 0 −𝐾3 𝐾3 + 𝐾4 −𝐾4 3
0 0 0 𝑚4 𝑢̈ 4 0 0 −𝐾4 𝐾4 𝑢4

[𝑀] ⋅ {𝑢̈ } + [𝐾] ⋅ {𝑢} = 0

OBSERVACIÓN:

Se observa que en la matriz K hay términos nulos que no coinciden con la matriz K que fueron dadas
como dato de entrada.

c) Analizamos la masa en cada piso

407.70 kg/m2
m5

407.70 kg/m2 m4

407.70 kg/m2 m2

407.70 kg/m2 m1
Primero analizamos la masa por cada piso y luego le sumamos la columna y el muro pantalla que tiene
por encima.

Losa

𝑚𝐿 = 2400 ⋅ 18 ⋅ 11 ⋅ 0,18 = 85536𝐾𝑔

Encadenado de Vigas

𝑚𝑉 = 2400 ⋅ 0.2 ⋅ 0.6 ⋅ (18 ⋅ 3 + 11 ⋅ 4) = 28224𝐾𝑔

Carga adicional sobre Losa

𝑚𝑎 = 407.70 ⋅ 18 ⋅ 11 = 80724.6𝐾𝑔

Total

𝑚 = 85536 + 28224 + 80724.6 = 194484.6𝐾𝑔

Masa de la Primera Planta

𝑚1 = 194484.6 + 12(2400 ⋅ 0.3 ⋅ 0.4 ⋅ 3.5) + 2400 ⋅ 6 ⋅ 0.15 ⋅ 3.5

𝑚1 = 214140.6𝐾𝑔

Masa de la Segunda Planta

𝑚2 = 214140.6𝐾𝑔

Masa de la Tercera Planta

𝑚3 = 214140.6𝐾𝑔

Masa de la Cuarta Planta

𝑚4 = 194484.6𝐾𝑔

La matriz de Masa será el mismo para el análisis de ambas direcciones

214140.6 0 0 0
[𝑀] = [ 0 214140.6 0 0 ]
0 0 214140.6 0
0 0 0 194484.6
d) Analizamos la Dirección ZY

1.638 ⋅ 109 −7.768 ⋅ 108 1.008 ⋅ 108 2.978 ⋅ 107

8
[𝐾] = [−7.768 ⋅ 108 1.215 ⋅ 109 −7.646 ⋅ 108 1.655 ⋅ 108 ]
1.008 ⋅ 10 −7.646 ⋅ 108 1.044 ⋅ 109 −4.342 ⋅ 108
2.978 ⋅ 107 1.655 ⋅ 108 −4.342 ⋅ 108 2.502 ⋅ 108
|[𝐾] − 𝜆 ∙ [𝑀]| = 0

1.638 ⋅ 109 −7.768 ⋅ 108 1.008 ⋅ 108 2.978 ⋅ 107 214140.6 0 0 0

8
|[−7.768 ⋅ 108 1.215 ⋅ 109 8 8
−7.646 ⋅ 10 1.655 ⋅ 10 ] − 𝜆 ∙ [ 0 214140.6 0 0 ]| = 0
1.008 ⋅ 10 −7.646 ⋅ 108 1.044 ⋅ 10 9
−4.342 ⋅ 10 8 0 0 214140.6 0
2.978 ⋅ 107 1.655 ⋅ 108 −4.342 ⋅ 108 2.502 ⋅ 108 0 0 0 194484.6

1.638 ⋅ 109 −7.768 ⋅ 108 1.008 ⋅ 108 2.978 ⋅ 107 214140.6 ∙ 𝜆 0 0 0

8
|[−7.768 ⋅ 108 1.215 ⋅ 109 8 8
−7.646 ⋅ 10 1.655 ⋅ 10 ] − [ 0 214140.6 ∙ 𝜆 0 0 ]| = 0
1.008 ⋅ 10 −7.646 ⋅ 108 1.044 ⋅ 10 9
−4.342 ⋅ 10 8 0 0 214140.6 ∙ 𝜆 0
2.978 ⋅ 107 1.655 ⋅ 108 −4.342 ⋅ 108 2.502 ⋅ 108 0 0 0 194484.6 ∙ 𝜆

1.638 ⋅ 109 − 214140.6 ∙ 𝜆 −7.768 ⋅ 108 1.008 ⋅ 108 2.978 ⋅ 107

| −7.768 ⋅ 108 1.215 ⋅ 109 − 214140.6 ∙ 𝜆 −7.646 ⋅ 108 1.655 ⋅ 108 |=0
1.008 ⋅ 108 −7.646 ⋅ 108 1.044 ⋅ 109 − 214140.6 ∙ 𝜆 −4.342 ⋅ 108
2.978 ⋅ 107 1.655 ⋅ 108 −4.342 ⋅ 108 2.502 ⋅ 108 − 194484.6 ∙ 𝜆

Factorizamos el valor de 107 para reducir el tamaño de los coeficientes obtenidos después de resolver el determinante

163.8 − 0.021414 ∙ 𝜆 −77.68 10.08 2.978

| −77.68 121.5 − 0.021414 ∙ 𝜆 −76.46 16.55 |=0
10.08 −76.46 104.4 − 0.021414 ∙ 𝜆 −43.42
2.978 16.55 −43.42 25.02 − 0.0194485 ∙ 𝜆

Resolviendo el determinante obtenemos la siguiente ecuación de 4to grado

 1.9098 ∙ 10−7 𝜆4 − 0.003721 ∙ 𝜆3 + 19.1799 ∙ 𝜆2 − 20822.6513 ∙ 𝜆 + 1.4601 ∙ 106 = 0

Resolviendo la ecuación de cuarto grado obtenemos los Eigenvalores

𝜆1 = 75.262
𝜆2 = 1367.124
𝜆3 = 6361.903
𝜆4 = 11679.426
Con los valores de λ, procedemos a calcular los Eigenvectores

Reemplazamos λ1 = 75.262

163,8 − 0,021414 ∙ 𝟕𝟓. 𝟐𝟔𝟐 −77,68 10,08 2,978 𝑈1

−77,68 121,5 − 0,021414 ∙ 𝟕𝟓. 𝟐𝟔𝟐 −76,46 16,55 𝑈
| | ∙ [ 2]
10,08 −76,46 104,4 − 0,021414 ∙ 𝟕𝟓. 𝟐𝟔𝟐 −43,42 𝑈3
2,978 16,55 −43,42 25,02 − 0,0194485 ∙ 𝟕𝟓. 𝟐𝟔𝟐 𝑈4

El sistema de ecuaciones obtenidos es:

162.188 ∙ 𝑈1 − 77.68 ∙ 𝑈2 + 10.08 ∙ 𝑈3 + 2.978 ∙ 𝑈4 = 0
−77.68 ∙ 𝑈1 + 119.888 ∙ 𝑈2 − 76.46 ∙ 𝑈3 + 16.55 ∙ 𝑈4 = 0
10.08 ∙ 𝑈1 − 76.46 ∙ 𝑈2 + 102.788 ∙ 𝑈3 − 43.42 ∙ 𝑈4 = 0
2.978 ∙ 𝑈1 + 16.55 ∙ 𝑈2 − 43.42 ∙ 𝑈3 + 23.556 ∙ 𝑈4 = 0
Para normalizar U4 = 1

162.188 ∙ 𝑈1 − 77.68 ∙ 𝑈2 + 10.08 ∙ 𝑈3 + 2.978 = 0

−77.68 ∙ 𝑈1 + 119.888 ∙ 𝑈2 − 76.46 ∙ 𝑈3 + 16.55 = 0
10.08 ∙ 𝑈1 − 76.46 ∙ 𝑈2 + 102.788 ∙ 𝑈3 − 43.42 = 0
2.978 ∙ 𝑈1 + 16.55 ∙ 𝑈2 − 43.42 ∙ 𝑈3 + 23.556 = 0
Resolviendo el sistema de ecuaciones obtenemos el primer Eigenvector:
U1 = 0.124 0.124
U2 = 0.388 {∅1 } = {0.388}
0.699
U3 = 0.699
1
U4 =1
Reemplazamos λ2 = 1367.124

163,8 − 0,021414 ∙ 𝟏𝟑𝟔𝟕. 𝟏𝟐𝟒 −77,68 10,08 2,978 𝑈1

−77,68 121,5 − 0,021414 ∙ 𝟏𝟑𝟔𝟕. 𝟏𝟐𝟒 −76,46 16,55 𝑈
| | ∙ [ 2]
10,08 −76,46 104,4 − 0,021414 ∙ 𝟏𝟑𝟔𝟕. 𝟏𝟐𝟒 −43,42 𝑈3
2,978 16,55 −43,42 25,02 − 0,0194485 ∙ 𝟏𝟑𝟔𝟕. 𝟏𝟐𝟒 𝑈4

El sistema de ecuaciones obtenidos es:

134.524 ∙ 𝑈1 − 77.68 ∙ 𝑈2 + 10.08 ∙ 𝑈3 + 2.978 ∙ 𝑈4 = 0
−77.68 ∙ 𝑈1 + 92.224 ∙ 𝑈2 − 76.46 ∙ 𝑈3 + 16.55 ∙ 𝑈4 = 0
10.08 ∙ 𝑈1 − 76.46 ∙ 𝑈2 + 75.124 ∙ 𝑈3 − 43.42 ∙ 𝑈4 = 0
2.978 ∙ 𝑈1 + 16.55 ∙ 𝑈2 − 43.42 ∙ 𝑈3 − 1.569 ∙ 𝑈4 = 0
Para normalizar U4 = 1
134.524 ∙ 𝑈1 − 77.68 ∙ 𝑈2 + 10.08 ∙ 𝑈3 + 2.978 = 0
−77.68 ∙ 𝑈1 + 92.224 ∙ 𝑈2 − 76.46 ∙ 𝑈3 + 16.55 = 0
10.08 ∙ 𝑈1 − 76.46 ∙ 𝑈2 + 75.124 ∙ 𝑈3 − 43.42 = 0
2.978 ∙ 𝑈1 + 16.55 ∙ 𝑈2 − 43.42 ∙ 𝑈3 − 1.569 = 0
Resolviendo el sistema de ecuaciones obtenemos el segundo Eigenvector:

U1 = -0.661 −0.661
U2 = -1.175 {∅2 } = {−1.175}
−0.529
U3 = -0.529 1
U4 =1
Reemplazamos λ3 = 6361.903

163,8 − 0,021414 ∙ 𝟔𝟑𝟔𝟏. 𝟗𝟎𝟑 −77,68 10,08 2,978 𝑈1

−77,68 121,5 − 0,021414 ∙ 𝟔𝟑𝟔𝟏. 𝟗𝟎𝟑 −76,46 16,55 𝑈
| | ∙ [ 2]
10,08 −76,46 104,4 − 0,021414 ∙ 𝟔𝟑𝟔𝟏. 𝟗𝟎𝟑 −43,42 𝑈3
2,978 16,55 −43,42 25,02 − 0,0194485 ∙ 𝟔𝟑𝟔𝟏. 𝟗𝟎𝟑 𝑈4

El sistema de ecuaciones obtenidos es:

27.566 ∙ 𝑈1 − 77.68 ∙ 𝑈2 + 10.08 ∙ 𝑈3 + 2.978 ∙ 𝑈4 = 0
−77.68 ∙ 𝑈1 − 14.73379 ∙ 𝑈2 − 76.46 ∙ 𝑈3 + 16.55 ∙ 𝑈4 = 0
10.08 ∙ 𝑈1 − 76.46 ∙ 𝑈2 − 31.83379 ∙ 𝑈3 − 43.42 ∙ 𝑈4 = 0
2.978 ∙ 𝑈1 + 16.55 ∙ 𝑈2 − 43.42 ∙ 𝑈3 − 98.70947 ∙ 𝑈4 = 0
Para normalizar U4 = 1

27.566 ∙ 𝑈1 − 77.68 ∙ 𝑈2 + 10.08 ∙ 𝑈3 + 2.978 = 0

−77.68 ∙ 𝑈1 − 14.73379 ∙ 𝑈2 − 76.46 ∙ 𝑈3 + 16.55 = 0
10.08 ∙ 𝑈1 − 76.46 ∙ 𝑈2 − 31.83379 ∙ 𝑈3 − 43.42 = 0
2.978 ∙ 𝑈1 + 16.55 ∙ 𝑈2 − 43.42 ∙ 𝑈3 − 98.70947 = 0
Resolviendo el sistema de ecuaciones obtenemos el tercer Eigenvector:

U1 = 2.023 2.023
U2 = 0.505 {∅3 } = { 0.505 }
−1.936
U3 = -1.936 1
U4 =1
Reemplazamos λ4 =11679.426
163,8 − 0,021414 ∙ 𝟏𝟏𝟔𝟕𝟗. 𝟒𝟐𝟔 −77,68 10,08 2,978 𝑈1
−77,68 121,5 − 0,021414 ∙ 𝟏𝟏𝟔𝟕𝟗. 𝟒𝟐𝟔 −76,46 16,55 𝑈2
| |∙[ ]
10,08 −76,46 104,4 − 0,021414 ∙ 𝟏𝟏𝟔𝟕𝟗. 𝟒𝟐𝟔 −43,42 𝑈3
2,978 16,55 −43,42 25,02 − 0,0194485 ∙ 𝟏𝟏𝟔𝟕𝟗. 𝟒𝟐𝟔 𝑈4
El sistema de ecuaciones obtenidos es:
−86.303 ∙ 𝑈1 − 77.68 ∙ 𝑈2 + 10.08 ∙ 𝑈3 + 2.978 ∙ 𝑈4 = 0
−77.68 ∙ 𝑈1 − 128.603 ∙ 𝑈2 − 76.46 ∙ 𝑈3 + 16.55 ∙ 𝑈4 = 0
10.08 ∙ 𝑈1 − 76.46 ∙ 𝑈2 − 145.703 ∙ 𝑈3 − 43.42 ∙ 𝑈4 = 0
2.978 ∙ 𝑈1 + 16.55 ∙ 𝑈2 − 43.42 ∙ 𝑈3 − 202.127 ∙ 𝑈4 = 0
Para normalizar U4 = 1

−86.303 ∙ 𝑈1 − 77.68 ∙ 𝑈2 + 10.08 ∙ 𝑈3 + 2.978 = 0

−77.68 ∙ 𝑈1 − 128.603 ∙ 𝑈2 − 76.46 ∙ 𝑈3 + 16.55 = 0
10.08 ∙ 𝑈1 − 76.46 ∙ 𝑈2 − 145.703 ∙ 𝑈3 − 43.42 = 0
2.978 ∙ 𝑈1 + 16.55 ∙ 𝑈2 − 43.42 ∙ 𝑈3 − 202.127 = 0
Resolviendo el sistema de ecuaciones obtenemos el cuarto Eigenvector:

U1 = -4.731 −4.731
U2 = 4.881 {∅4 } = { 4.881 }
−3.186
U3 = -3.186 1
U4 =1
Sintetizamos los Eigenvectores
0.124 −0.661 2.023 −4.731
[∅]𝑍𝑌 = [0.388 −1.175 0.505 4.881 ]
0.699 −0.529 −1.936 −3.186
1 1 1 1
e) Analizamos la Dirección ZX

2.061 ⋅ 108 −8.478 ⋅ 107 7.969 ⋅ 106 −2.54 ⋅ 106

7
[𝐾] = [−8.478 ⋅ 106 1.864 ⋅ 108 −8.986 ⋅ 107 4.637 ⋅ 106 ]
7.969 ⋅ 10 −8.986 ⋅ 107 1.725 ⋅ 108 −7.813 ⋅ 107
−2.54 ⋅ 106 4.637 ⋅ 106 −7.813 ⋅ 107 7.537 ⋅ 107
 |[𝐾] − 𝜆 ∙ [𝑀]| = 0

2.061 ⋅ 108 −8.478 ⋅ 107 7.969 ⋅ 106 −2.54 ⋅ 106 214140.6 0 0 0

7
|[−8.478 ⋅ 106 1.864 ⋅ 108 7 6
−8.986 ⋅ 10 4.637 ⋅ 10 ] − 𝜆 ∙ [ 0 214140.6 0 0 ]| = 0
7.969 ⋅ 10 −8.986 ⋅ 107 1.725 ⋅ 10 8
−7.813 ⋅ 10 7 0 0 214140.6 0
−2.54 ⋅ 106 4.637 ⋅ 106 −7.813 ⋅ 107 7.537 ⋅ 107 0 0 0 194484.6

2.061 ⋅ 108 −8.478 ⋅ 107 7.969 ⋅ 106 −2.54 ⋅ 106 214140.6 ∙ 𝜆 0 0 0

7
|[−8.478 ⋅ 106 1.864 ⋅ 108 7 6
−8.986 ⋅ 10 4.637 ⋅ 10 ] − [ 0 214140.6 ∙ 𝜆 0 0 ]| = 0
7.969 ⋅ 10 −8.986 ⋅ 107 1.725 ⋅ 10 8
−7.813 ⋅ 10 7 0 0 214140.6 ∙ 𝜆 0
−2.54 ⋅ 106 4.637 ⋅ 106 −7.813 ⋅ 10 7
7.537 ⋅ 10 7 0 0 0 194484.6 ∙ 𝜆

2.061 ⋅ 108 − 214140.6 ∙ 𝜆 −8.478 ⋅ 107 7.969 ⋅ 106 −2.54 ⋅ 106

| −8.478 ⋅ 107 1.864 ⋅ 108 − 214140.6 ∙ 𝜆 −8.986 ⋅ 107 4.637 ⋅ 106 |=0
7.969 ⋅ 106 −8.986 ⋅ 107 1.725 ⋅ 108 − 214140.6 ∙ 𝜆 −7.813 ⋅ 107
−2.54 ⋅ 106 4.637 ⋅ 106 −7.813 ⋅ 107 7.537 ⋅ 107 − 194484.6 ∙ 𝜆

Factorizamos el valor de 106 para reducir el tamaño de los coeficientes obtenidos después de resolver el determinante

206.1 − 0.2141406 ∙ 𝜆 −84.78 7.969 −2.54

| −84.78 186.4 − 0.2141406 ∙ 𝜆 −89.86 4.637 |=0
7.969 −89.86 172.5 − 0.2141406 ∙ 𝜆 −78.13
−2.54 4.637 −78.13 75.37 − 0.1944846 ∙ 𝜆

Resolviendo el determinante obtenemos la siguiente ecuación de 4to grado

1.9098 ∙ 10−3 𝜆4 − 5.77896 ∙ 𝜆3 + 5452.9812 ∙ 𝜆2 − 1.71054 ∙ 106 ∙ 𝜆 + 1.08122 ∙ 108 = 0

Resolviendo la ecuación de cuarto grado obtenemos los Eigenvalores

𝜆1 = 83.543
𝜆2 = 464.905
𝜆3 = 961.462
𝜆4 = 1516.088
Con los valores de λ, procedemos a calcular los Eigenvectores

Reemplazamos λ1 = 83.543

206.1 − 0.2141406 ∙ 𝟖𝟑. 𝟓𝟒𝟑 −84.78 7.969 −2.54 𝑈1

−84.78 186.4 − 0.2141406 ∙ 𝟖𝟑. 𝟓𝟒𝟑 −89.86 4.637 𝑈
| | ∙ [ 2]
7.969 −89.86 172.5 − 0.2141406 ∙ 𝟖𝟑. 𝟓𝟒𝟑 −78.13 𝑈3
−2.54 4.637 −78.13 75.37 − 0.1944846 ∙ 𝟖𝟑. 𝟓𝟒𝟑 𝑈4

El sistema de ecuaciones obtenidos es:

188.21 ∙ 𝑈1 − 84.78 ∙ 𝑈2 + 7.969 ∙ 𝑈3 − 2.54 ∙ 𝑈4 = 0
−84.78 ∙ 𝑈1 + 168.51 ∙ 𝑈2 − 89.86 ∙ 𝑈3 + 4.637 ∙ 𝑈4 = 0
7.969 ∙ 𝑈1 − 89.86 ∙ 𝑈2 + 154.61 ∙ 𝑈3 − 78.13 ∙ 𝑈4 = 0
−2.54 ∙ 𝑈1 + 4.637 ∙ 𝑈2 − 78.13 ∙ 𝑈3 + 59.122 ∙ 𝑈4 = 0
Para normalizar U4 = 1
188.21 ∙ 𝑈1 − 84.78 ∙ 𝑈2 + 7.969 ∙ 𝑈3 − 2.54 = 0
−84.78 ∙ 𝑈1 + 168.51 ∙ 𝑈2 − 89.86 ∙ 𝑈3 + 4.637 = 0
7.969 ∙ 𝑈1 − 89.86 ∙ 𝑈2 + 154.61 ∙ 𝑈3 − 78.13 = 0
−2.54 ∙ 𝑈1 + 4.637 ∙ 𝑈2 − 78.13 ∙ 𝑈3 + 59.122 = 0
Resolviendo el sistema de ecuaciones obtenemos el primer Eigenvector:

U1 = 0.201 0.201
U2 = 0.489 {∅1 } = {0.489}
0.779
U3 = 0.779
1
U4 =1
Reemplazamos λ2 = 464.905

206.1 − 0.2141406 ∙ 𝟒𝟔𝟒. 𝟗𝟎𝟓 −84.78 7.969 −2.54 𝑈1

−84.78 186.4 − 0.2141406 ∙ 𝟒𝟔𝟒. 𝟗𝟎𝟓 −89.86 4.637 𝑈
| | ∙ [ 2]
7.969 −89.86 172.5 − 0.2141406 ∙ 𝟒𝟔𝟒. 𝟗𝟎𝟓 −78.13 𝑈3
−2.54 4.637 −78.13 75.37 − 0.1944846 ∙ 𝟒𝟔𝟒. 𝟗𝟎𝟓 𝑈4

El sistema de ecuaciones obtenidos es:

106.545 ∙ 𝑈1 − 84.78 ∙ 𝑈2 + 7.969 ∙ 𝑈3 − 2.54 ∙ 𝑈4 = 0
−84.78 ∙ 𝑈1 + 86.845 ∙ 𝑈2 − 89.86 ∙ 𝑈3 + 4.637 ∙ 𝑈4 = 0
7.969 ∙ 𝑈1 − 89.86 ∙ 𝑈2 + 72.945 ∙ 𝑈3 − 78.13 ∙ 𝑈4 = 0
−2.54 ∙ 𝑈1 + 4.637 ∙ 𝑈2 − 78.13 ∙ 𝑈3 − 15.047 ∙ 𝑈4 = 0

Para normalizar U4 = 1

106.545 ∙ 𝑈1 − 84.78 ∙ 𝑈2 + 7.969 ∙ 𝑈3 − 2.54 = 0

−84.78 ∙ 𝑈1 + 86.845 ∙ 𝑈2 − 89.86 ∙ 𝑈3 + 4.637 = 0
7.969 ∙ 𝑈1 − 89.86 ∙ 𝑈2 + 72.945 ∙ 𝑈3 − 78.13 = 0
−2.54 ∙ 𝑈1 + 4.637 ∙ 𝑈2 − 78.13 ∙ 𝑈3 − 15.047 = 0

Resolviendo el sistema de ecuaciones obtenemos el segundo Eigenvector:

U1 = -0.861 −0.861
U2 = -1.134 {∅2 } = {−1.134}
−0.232
U3 = -0.232 1
U4 =1
Reemplazamos λ3 = 961.462

206.1 − 0.2141406 ∙ 𝟗𝟔𝟏. 𝟒𝟔𝟐 −84.78 7.969 −2.54 𝑈1

−84.78 186.4 − 0.2141406 ∙ 𝟗𝟔𝟏. 𝟒𝟔𝟐 −89.86 4.637 𝑈
| | ∙ [ 2]
7.969 −89.86 172.5 − 0.2141406 ∙ 𝟗𝟔𝟏. 𝟒𝟔𝟐 −78.13 𝑈3
−2.54 4.637 −78.13 75.37 − 0.1944846 ∙ 𝟗𝟔𝟏. 𝟒𝟔𝟐 𝑈4

El sistema de ecuaciones obtenidos es:

0.212 ∙ 𝑈1 − 84.78 ∙ 𝑈2 + 7.969 ∙ 𝑈3 − 2.54 ∙ 𝑈4 = 0
−84.78 ∙ 𝑈1 − 19.488 ∙ 𝑈2 − 89.86 ∙ 𝑈3 + 4.637 ∙ 𝑈4 = 0
7.969 ∙ 𝑈1 − 89.86 ∙ 𝑈2 − 33.388 ∙ 𝑈3 − 78.13 ∙ 𝑈4 = 0
−2.54 ∙ 𝑈1 + 4.637 ∙ 𝑈2 − 78.13 ∙ 𝑈3 − 111.619 ∙ 𝑈4 = 0
Para normalizar U4 = 1
0.212 ∙ 𝑈1 − 84.78 ∙ 𝑈2 + 7.969 ∙ 𝑈3 − 2.54 = 0
−84.78 ∙ 𝑈1 − 19.488 ∙ 𝑈2 − 89.86 ∙ 𝑈3 + 4.637 = 0
7.969 ∙ 𝑈1 − 89.86 ∙ 𝑈2 − 33.388 ∙ 𝑈3 − 78.13 = 0
−2.54 ∙ 𝑈1 + 4.637 ∙ 𝑈2 − 78.13 ∙ 𝑈3 − 111.619 = 0
Resolviendo el sistema de ecuaciones obtenemos el tercer Eigenvector:

U1 = 1.675 1.675
U2 = -0.166 {∅3 } = {−0.166}
−1.493
U3 = -1.493 1
U4 =1
Reemplazamos λ4 =1516.088

206.1 − 0.2141406 ∙ 𝟏𝟓𝟏𝟔. 𝟎𝟖𝟖 −84.78 7.969 −2.54 𝑈1

−84.78 186.4 − 0.2141406 ∙ 𝟏𝟓𝟏𝟔. 𝟎𝟖𝟖 −89.86 4.637 𝑈
| | ∙ [ 2]
7.969 −89.86 172.5 − 0.2141406 ∙ 𝟏𝟓𝟏𝟔. 𝟎𝟖𝟖 −78.13 𝑈3
−2.54 4.637 −78.13 75.37 − 0.1944846 ∙ 𝟏𝟓𝟏𝟔. 𝟎𝟖𝟖 𝑈4
El sistema de ecuaciones obtenidos es:
−118.556 ∙ 𝑈1 − 84.78 ∙ 𝑈2 + 7.969 ∙ 𝑈3 − 2.54 ∙ 𝑈4 = 0
−84.78 ∙ 𝑈1 − 138.256 ∙ 𝑈2 − 89.86 ∙ 𝑈3 + 4.637 ∙ 𝑈4 = 0
7.969 ∙ 𝑈1 − 89.86 ∙ 𝑈2 − 152.156 ∙ 𝑈3 − 78.13 ∙ 𝑈4 = 0
−2.54 ∙ 𝑈1 + 4.637 ∙ 𝑈2 − 78.13 ∙ 𝑈3 − 219.486 ∙ 𝑈4 = 0
Para normalizar U4 = 1
−118.556 ∙ 𝑈1 − 84.78 ∙ 𝑈2 + 7.969 ∙ 𝑈3 − 2.54 ∙ 𝑈4 = 0
−84.78 ∙ 𝑈1 − 138.256 ∙ 𝑈2 − 89.86 ∙ 𝑈3 + 4.637 ∙ 𝑈4 = 0
7.969 ∙ 𝑈1 − 89.86 ∙ 𝑈2 − 152.156 ∙ 𝑈3 − 78.13 ∙ 𝑈4 = 0
−2.54 ∙ 𝑈1 + 4.637 ∙ 𝑈2 − 78.13 ∙ 𝑈3 − 219.486 ∙ 𝑈4 = 0
Resolviendo el sistema de ecuaciones obtenemos el cuarto Eigenvector:

U1 = -2.485 −2.485
U2 = 3.206 {∅4 } = { 3.206 }
−2.537
U3 = -2.537 1
U4 =1
Sintetizamos los Eigenvectores
0.201 −0.861 1.675 −2.485
[∅]𝑍𝑋 = [0.489 −1.134 −0.166 3.206 ]
0.779 −0.232 −1.493 −2.537
1 1 1 1
2DO: PERIODOS DE VIBRACIÓN PARA CADA MODO EN AMBAS DIRECCIONES PRINCIPALES.
a) Analizamos la Dirección ZY

𝜆1 = 75.262
𝜆2 = 1367.124 𝑊𝑛 = √𝜆
𝜆3 = 6361.903 𝑊𝑛1 = √75.262 = 8.675 𝑟𝑎𝑑/𝑠𝑒𝑔
𝜆4 = 11679.426
𝑊𝑛2 = √1367.124 = 36.975 𝑟𝑎𝑑/𝑠𝑒𝑔

𝑊𝑛3 = √6361.903 = 79.762 𝑟𝑎𝑑/seg

𝑊𝑛4 = √11679.426 = 108.071 𝑟𝑎𝑑/𝑠𝑒𝑔

8.675
[𝑊𝑛 ] = { 36.975 } 𝑟𝑎𝑑/𝑠𝑒𝑔
79.762
108.071

2𝜋
𝑇𝑛 =
𝑊𝑛
2𝜋
𝑇𝑛1 = = 0.724 𝑠𝑒𝑔
8.675
2𝜋
𝑇𝑛2 = = 0.170 𝑠𝑒𝑔
36.975
2𝜋
𝑇𝑛3 = = 0.079 𝑠𝑒𝑔
79.762
2𝜋
𝑇𝑛4 = = 0.058 𝑠𝑒𝑔
108.071

0.724
[𝑇] = {0.170} 𝑠𝑒𝑔
0.079
0.058
b) Analizamos la Dirección ZX

𝜆1 = 83.543
𝜆2 = 464.905 𝑊𝑛 = √𝜆
𝜆3 = 961.462 𝑊𝑛1 = √83.543 = 9.140 𝑟𝑎𝑑/𝑠𝑒𝑔
𝜆4 = 1516.088
𝑊𝑛2 = √464.905 = 21.562 𝑟𝑎𝑑/𝑠𝑒𝑔

𝑊𝑛3 = √961.462 = 31.007 𝑟𝑎𝑑/seg

𝑊𝑛4 = √1516.088 = 38.937 𝑟𝑎𝑑/𝑠𝑒𝑔

9.140
[𝑊𝑛 ] = { 21.562 } 𝑟𝑎𝑑/𝑠𝑒𝑔
31.007
38.937

2𝜋
𝑇𝑛 =
𝑊𝑛
2𝜋
𝑇𝑛1 = = 0.687 𝑠𝑒𝑔
9.140
2𝜋
𝑇𝑛2 = = 0.291 𝑠𝑒𝑔
21.562
2𝜋
𝑇𝑛3 = = 0.203 𝑠𝑒𝑔
31.007
2𝜋
𝑇𝑛4 = = 0.161 𝑠𝑒𝑔
38.937

0.687
[𝑇] = {0.291} 𝑠𝑒𝑔
0.203
0.161
3ERO: FORMAS DE LOS MODOS DE VIBRACIÓN EN AMBAS DIRECCIONES PRINCIPALES (NORMALIZADAS).
a) Para el Plano ZY
0.124 −0.661 2.023 −4.731
[∅]𝑍𝑌 = [0.388 −1.175 0.505 4.881 ]
0.699 −0.529 −1.936 −3.186
1 1 1 1

Representamos gráficamente

1 1 1 1

0.699 0.529 1.936 3.186

1.175 0.505 4.881

0.388

2.023 4.731
0.124 0.661

Modo 1 Modo 2 Modo 3 Modo 4

 b) Para el Plano ZX
0.201 −0.861 1.675 −2.485
[∅]𝑍𝑋 = [0.489 −1.134 −0.166 3.206 ]
0.779 −0.232 −1.493 −2.537
1 1 1 1

Representamos gráficamente

1 1 1 1

0.232 2.537
0.779 1.493

0.489 1.134 0.166 3.206

1.675 2.485
0.201 0.861

Modo 1 Modo 2 Modo 3 Modo 4

A LA DIRECCIÓN DEL PLANO ZY, APLICAR EL ESPECTRO DE RESPUESTA DE ACELERACIÓN DEL SISMO DE KOBE-
JAPON (NGA1111_KOVEJAPAN_NISHIAKASHI_00)

4TO: DETERMINAR LAS ECUACIONES DESACOPLADAS DE CADA MODO DE VIBRACIÓN

[𝑀∗ ] = [∅]𝑇 ∙ [𝑀] ∙ [∅]

0.124 0.388 0.699 1 214140.6 0 0 0 0.124 −0.661 2.023 −4.731

[𝑀 ∗]
=[ −0.661 −1.175 −0.529 1] ∙ [ 0 214140.6 0 0 ] ∙ [0.388 −1.175 0.505 4.881 ]
2.023 0.505 −1.936 1 0 0 214140.6 0 0.699 −0.529 −1.936 −3.186
−4.731 4.881 −3.186 1 0 0 0 194484.6 1 1 1 1

𝟑𝟑𝟒𝟔𝟒𝟒. 𝟏 123.1 372.1 −2488.3

[𝑀∗ ] = [ 123.1 𝟔𝟒𝟑𝟔𝟐𝟎. 𝟑 380.3 −3079.4 ]
372.1 380.3 𝟏𝟗𝟐𝟖𝟎𝟗𝟏. 𝟗 −6339
−2488.3 −3079.4 −6339 𝟏𝟐𝟐𝟔𝟐𝟖𝟑𝟏. 𝟕
Cálculo del Factor de participación Modal

1
{∅𝑖 }𝑇 1
𝑃𝑖 = − ∗ ∙ [𝑀] ∙ {1}
𝑚𝑖
1
214140.6 0 0 0 1
{0.124 0.388 0.699 1} 0 214140.6
𝑃1 = − ∙[ 0 0 ] ∙ {1} = −1.356
334644.1 0 0 214140.6 0 1
0 0 0 194484.6 1

214140.6 0 0 0 1
{−0.661 −1.175 −0.529 1} 0 214140.6
𝑃2 = − ∙[ 0 0 ] ∙ {1} = 0.485
643620.3 0 0 214140.6 0 1
0 0 0 194484.6 1

214140.6 0 0 0 1
{2.023 0.505 −1.936 1} 0 214140.6
𝑃3 = − ∙[ 0 0 ] ∙ {1} = −0.167
1928091.9 0 0 214140.6 0 1
0 0 0 194484.6 1
 214140.6 0 0 0 1
{−4.731 4.881 −3.186 1} 0 214140.6
𝑃4 = − ∙[ 0 0 ] ∙ {1} = 0.037
12262831.7 0 0 214140.6 0 1
0 0 0 194484.6 1
Ecuación del Primer modo de vibración
2
𝑦̈ 1(𝑡) + 2 ∙ 𝜍 ∙ 𝑊𝑛1 ∙ 𝑦̇ 1(𝑡) + 𝑊𝑛1 ∙ 𝑦1(𝑡) = 𝑃1 ∙ 𝑢̈ 𝑔(𝑡)
Datos
ς=0.05
Wn1=8.675 rad/seg
P1=-1.356
𝑦̈ 1(𝑡) + 2 ∙ 0.05 ∙ 8.675 ∙ 𝑦̇1(𝑡) + 8.6752 ∙ 𝑦1(𝑡) = −1.356 ∙ 𝑢̈ 𝑔(𝑡)

𝑦̈ 1(𝑡) + 0.8675 ∙ 𝑦̇1(𝑡) + 75.256 ∙ 𝑦1(𝑡) = −1.356 ∙ 𝑢̈ 𝑔(𝑡)

Ecuación del Segundo modo de vibración

2
𝑦̈ 2(𝑡) + 2 ∙ 𝜍 ∙ 𝑊𝑛2 ∙ 𝑦̇ 1(𝑡) + 𝑊𝑛2 ∙ 𝑦1(𝑡) = 𝑃2 ∙ 𝑢̈ 𝑔(𝑡)
Datos
ς=0.05
Wn2=36.975 rad/seg
P2=0.485
𝑦̈ 2(𝑡) + 2 ∙ 0.05 ∙ 36.975 ∙ 𝑦̇1(𝑡) + 36.9752 ∙ 𝑦1(𝑡) = 0.485 ∙ 𝑢̈ 𝑔(𝑡)

𝑦̈ 2(𝑡) + 3.6975 ∙ 𝑦̇ 1(𝑡) + 1367.15 ∙ 𝑦1(𝑡) = 0.485 ∙ 𝑢̈ 𝑔(𝑡)

Ecuación del Tercer modo de vibración
2
𝑦̈ 3(𝑡) + 2 ∙ 𝜍 ∙ 𝑊𝑛3 ∙ 𝑦̇ 1(𝑡) + 𝑊𝑛3 ∙ 𝑦1(𝑡) = 𝑃3 ∙ 𝑢̈ 𝑔(𝑡)
Datos
ς=0.05
Wn3=79.762 rad/seg
P3=-0.167
𝑦̈ 3(𝑡) + 2 ∙ 0.05 ∙ 79.762 ∙ 𝑦̇1(𝑡) + 79.7622 ∙ 𝑦1(𝑡) = −0.167 ∙ 𝑢̈ 𝑔(𝑡)

𝑦̈ 3(𝑡) + 7.9762 ∙ 𝑦̇ 1(𝑡) + 6361.98 ∙ 𝑦1(𝑡) = −0.167 ∙ 𝑢̈ 𝑔(𝑡)

Ecuación del Cuarto modo de vibración
2
𝑦̈4(𝑡) + 2 ∙ 𝜍 ∙ 𝑊𝑛4 ∙ 𝑦̇ 1(𝑡) + 𝑊𝑛4 ∙ 𝑦1(𝑡) = 𝑃4 ∙ 𝑢̈ 𝑔(𝑡)
Datos
ς=0.05
Wn4=108.071 rad/seg
P4=0.037
𝑦̈4(𝑡) + 2 ∙ 0.05 ∙ 108.071 ∙ 𝑦̇ 1(𝑡) + 108.0712 ∙ 𝑦1(𝑡) = 0.037 ∙ 𝑢̈ 𝑔(𝑡)

𝑦̈4(𝑡) + 10.8071 ∙ 𝑦̇1(𝑡) + 11679.34 ∙ 𝑦1(𝑡) = 0.037 ∙ 𝑢̈ 𝑔(𝑡)

Resumiendo, las ecuaciones obtenidas:

𝑦̈ 1(𝑡) + 0.8675 ∙ 𝑦̇1(𝑡) + 75.256 ∙ 𝑦1(𝑡) = −1.356 ∙ 𝑢̈ 𝑔(𝑡)

𝑦̈ 2(𝑡) + 3.6975 ∙ 𝑦̇ 1(𝑡) + 1367.15 ∙ 𝑦1(𝑡) = 0.485 ∙ 𝑢̈ 𝑔(𝑡)

𝑦̈ 3(𝑡) + 7.9762 ∙ 𝑦̇ 1(𝑡) + 6361.98 ∙ 𝑦1(𝑡) = −0.167 ∙ 𝑢̈ 𝑔(𝑡)

𝑦̈4(𝑡) + 10.8071 ∙ 𝑦̇1(𝑡) + 11679.34 ∙ 𝑦1(𝑡) = 0.037 ∙ 𝑢̈ 𝑔(𝑡)

5TO: DETERMINAR FACTORES DE PARTICIPACION MODAL

1
{∅𝑖 }𝑇 1
𝑃𝑖 = − ∗ ∙ [𝑀] ∙ {1}
𝑚𝑖
1
214140.6 0 0 0 1
{0.124 0.388 0.699 1} 0 214140.6
𝑃1 = − ∙[ 0 0 ] ∙ {1} = −1.356
334644.1 0 0 214140.6 0 1
0 0 0 194484.6 1

214140.6 0 0 0 1
{−0.661 −1.175 −0.529 1} 0 214140.6
𝑃2 = − ∙[ 0 0 ] ∙ {1} = 0.485
643620.3 0 0 214140.6 0 1
0 0 0 194484.6 1
 214140.6 0 0 0 1
{2.023 0.505 −1.936 1} 0 214140.6
𝑃3 = − ∙[ 0 0 ] ∙ {1} = −0.167
1928091.9 0 0 214140.6 0 1
0 0 0 194484.6 1

214140.6 0 0 0 1
{−4.731 4.881 −3.186 1} 0 214140.6
𝑃4 = − ∙[ 0 0 ] ∙ {1} = 0.037
12262831.7 0 0 214140.6 0 1
0 0 0 194484.6 1

6TO: DETERMINAR MASA MODAL EFECTIVA (ESTABLECER CUANTOS MODOS SERIAN SUFICIENTES PARA TENER
UNA PARTICIPACIÓN DE MASA EFECTIVA DE 90%)

̅ 𝑖 = 𝑃𝑖2 ∙ 𝑚𝑖∗
𝑚

̅ 1 = 𝑃12 ∙ 𝑚1∗ = (−1.356)2 ∙ 334644.1 = 615322.15

𝑚

̅ 2 = 𝑃22 ∙ 𝑚2∗ = 0.4852 ∙ 643620.3 = 151395.59

𝑚

̅ 3 = 𝑃32 ∙ 𝑚3∗ = (−0.167)2 ∙ 1928091.9 = 53772.555

𝑚

̅ 4 = 𝑃42 ∙ 𝑚4∗ = 0.0372 ∙ 12262831.7 = 16787.817

𝑚

Masa
Modos Masa % Acumulado
Acumulada
1 615322.15 615322.15 73.5
2 151395.59 766717.74 91.6
3 53772.555 820490.29 98
4 16787.817 837278.11 100
total 837278.11

Es suficiente tener dos modos de vibración para tener una participación efectiva superior al 90%
7MO: DETERMINAR LOS DESPLAZAMIENTO Y LA DERIVA EN EJES GLOBALES DE CADA NIVEL

1.52

0.88

0.73
0.64
0.56

0.724
0.058

0.170
0.079
Entramos al Espectro de Respuesta de Aceleración Sismo de Kobe Japón con los valores Tn, para
encontrar las aceleraciones

𝑇𝑛1 = 0.724 𝑠𝑒𝑔 𝑎1 = 0.73 𝑚/𝑠 2

𝑇𝑛2 = 0.170 𝑠𝑒𝑔 𝑎2 = 0.88 𝑚/𝑠 2
𝑇𝑛3 = 0.079 𝑠𝑒𝑔 𝑎3 = 0.64 𝑚/𝑠 2
𝑇𝑛4 = 0.058 𝑠𝑒𝑔 𝑎4 = 0.56 𝑚/𝑠 2

Con las aceleraciones y las frecuencias obtenemos los desplazamientos:

𝑎𝑖
Datos 𝑑𝑖 =
𝑊𝑛𝑖 2
𝑊𝑛1 = 8.675 𝑟𝑎𝑑/𝑠𝑒𝑔
0.73
𝑊𝑛2 = 36.975 𝑟𝑎𝑑/𝑠𝑒𝑔 𝑑1 = = 0.00970 𝑚
8.6752
𝑊𝑛3 = 79.762 𝑟𝑎𝑑/seg 0.88
𝑑2 = = 0.00064 𝑚
𝑊𝑛4 = 108.071 𝑟𝑎𝑑/𝑠𝑒𝑔 36.9752
0.64
𝑑3 = = 0.00010 𝑚
79.7622
0.56
𝑑4 = = 0.00005 𝑚
108.0712
Estos desplazamientos tienen que ser afectados por el factor de Participación

𝑦𝑖 = 𝑃𝑖 ∙ 𝑑𝑖

𝑦1 = 𝑃1 ∙ 𝑑1 = 1.356 ∙ 0.00970 = 0.013154𝑚

𝑦2 = 𝑃2 ∙ 𝑑2 = 0.485 ∙ 0.00064 = 0.000312𝑚

𝑦3 = 𝑃3 ∙ 𝑑3 = 0.167 ∙ 0.00010 = 0.000017𝑚

𝑦4 = 𝑃4 ∙ 𝑑4 = 0.037 ∙ 0.00005 = 0.000002𝑚

Convertimos los desplazamientos a Coordenadas Globales

{𝑢𝑖 } = {∅𝑖 } ∙ 𝑦𝑖
Para el modo 1
0.124 0.001631
{𝑢1 } = { 0.388 } ∙ 0.013154 = {0.005104}
0.699 0.009195
1 0.013154
Para el modo 2
−0.661 −0.000206
{𝑢2 } = { −1.175 } ∙ 0.000312 = {−0.000367}
−0.529 −0.000165
1 0.000312
Para el modo 3
2.023 0.000034
{𝑢3 } = { 0.505 } ∙ 0.000017 = { 0.000009 }
−1.936 −0.000033
1 0.000017
Para el modo 4
−4.731 −0.000009
{𝑢4 } = { 4.881 } ∙ 0.000002 = { 0.000010 }
−3.186 −0.000006
1 0.000002

A nuestros resultados aplicamos SRSS

√0.0016312 + (−0.000206)2 + 0.0000342 + (−0.000009)2 = 0.00164 𝑚

√0.0051042 + (−0.000367)2 + 0.0000092 + 0.0000102 = 0.00512 𝑚

√0.0091952 + (−0.000165)2 + (−0.000033)2 + (−0.000006)2 = 0.0092 𝑚

√0.0131542 + 0.0003122 + 0.0000172 + 0.0000022 = 0.01316𝑚

0.01316m

0.0092m

0.00512m

0.00164m
Calculamos la Deriva

Para el modo 1
0.001631 0.001631
{ 0.005104 − 0.001631 } = {0.003473}
0.009195 − 0.005104 0.004091
0.013154 − 0.009195 0.003959
Para el modo 2
−0.000206 −0.000206
−0.000367 − (−0.000206)
{ } = {−0.000160}
−0.000165 − (−0.000367) 0.000202
0.000312 − (−0.000165) 0.000477

Para el modo 3
0.000034 0.000034
0.000009 − 0.000034 −0.000026
{ −0.000033 − 0.000009 } = { }
−0.000041
0.000017 − (−0.000033) 0.000050
Para el modo 4
−0.000009 −0.000009
0.000010 − (−0.000009)
{ } = { 0.000019 }
−0.000006 − 0.000010 −0.000016
0.000002 − (−0.000006) 0.000008

A nuestros resultados aplicamos SRSS

√0.0016312 + (−0.000206)2 + 0.0000342 + (−0.000009)2 = 0.00164 𝑚

√0.0034732 + (−0.000160)2 + (−0.000026)2 + 0.0000192 = 0.00348 𝑚

√0.0040912 + 0.0002022 + (−0.000041)2 + (−0.000016)2 = 0.00410 𝑚

√0.0039592 + 0.0004772 + 0.0000502 + 0.0000082 = 0.00399 𝑚

𝑢1 0.00164
𝑢2 − 𝑢1
{𝑢 − 𝑢 } = {0.00348}
3 2 0.00410
𝑢4 − 𝑢3 0.00399
8VO: DETERMINAR LAS FUERZAS RESULTANTES EN CADA NIVEL

𝐹𝑖𝑛 = 𝑎𝑛𝑜 ∙ 𝑚𝑖 ∙ 𝑃𝑖 ∙ ∅𝑛

Para el modo 1

𝑎1 = 0.73 𝑚/𝑠 2

𝐹11 214140.6 0.124 257853

𝐹21 214140.6
{ } = 0.73 ∙ 9.81 ∙ { } ∙ 1.356 ∙ { 0.388 } = { 806830 }
𝐹31 214140.6 0.699 1453543
𝐹41 194484.6 1 1888586

Para el modo 2

𝑎2 = 0.88 𝑚/𝑠 2
𝐹12 214140.6 −0.661 −592644
𝐹22
{ } = 0.88 ∙ 9.81 ∙ {214140.6} ∙ 0.485 ∙ {−1.175} = {−1053490}
𝐹32 214140.6 −0.529 −474295
𝐹42 194484.6 1 814289

Para el modo 3

𝑎3 = 0.64 𝑚/𝑠 2
𝐹13 214140.6 2.023 454214
𝐹23
{ } = 0.64 ∙ 9.81 ∙ {214140.6} ∙ 0.167 ∙ { 0.505 } = { 113385 }
𝐹33 214140.6 −1.936 −434680
𝐹43 194484.6 1 203916

Para el modo 4

𝑎4 = 0.56 𝑚/𝑠 2
𝐹14 214140.6 −4.731 −205926
𝐹24 214140.6 4.881
{ } = 0.56 ∙ 9.81 ∙ { } ∙ 0.037 ∙ { } = { 212455 }
𝐹34 214140.6 −3.186 −138677
𝐹44 194484.6 1 39532
 OPCIONAL

DETERMINAR LAS ECUACIONES DE MOVIMIENTO EN COORDENADAS GLOBALES DE

CADA NIVEL SI SE TIENE UNA FUERZA APLICADA EN EL PISO 4

𝑭𝟒(𝒕) = 𝟏𝟎𝟎 ∙ 𝟏𝟎𝟑 ∙ 𝑺𝒆𝒏(𝟗𝟎 ∙ 𝒕) [N]

∗
100103·Sen(90·t) 𝐹(𝑡) = [∅]𝑇 ∙ {𝐹(𝑡) }

Analizamos la dirección ZY

0.124 0.388 0.699 1 0

∗ −0.661 −1.175 −0.529 1] ∙ { 0
𝐹(𝑡) =[ 0 }
2.023 0.505 −1.936 1
−4.731 4.881 −3.186 1 100103 ∙ sin(90 ∙ 𝑡)

100103 ∙ sin(90 ∙ 𝑡)
∗ 100103 ∙ sin(90 ∙ 𝑡)
𝐹(𝑡) ={ }
100103 ∙ sin(90 ∙ 𝑡)
100103 ∙ sin(90 ∙ 𝑡)

100103 ∙ sin(90 ∙ 𝑡)
𝑦̈1(𝑡) + 0.8675 ∙ 𝑦̇ 1(𝑡) + 75.256 ∙ 𝑦1(𝑡) =
334644.1
100103 ∙ sin(90 ∙ 𝑡)
𝑦̈ 2(𝑡) + 3.6975 ∙ 𝑦̇ 1(𝑡) + 1367.15 ∙ 𝑦1(𝑡) =
643620.3
100103 ∙ sin(90 ∙ 𝑡)
𝑦̈ 3(𝑡) + 7.9762 ∙ 𝑦̇ 1(𝑡) + 6361.98 ∙ 𝑦1(𝑡) =
1928091.9
100103 ∙ sin(90 ∙ 𝑡)
𝑦̈4(𝑡) + 10.8071 ∙ 𝑦̇1(𝑡) + 11679.34 ∙ 𝑦1(𝑡) =
12262831.7

Haciendo operaciones obtenemos:

𝑦̈ 1(𝑡) + 0.8675 ∙ 𝑦̇ 1(𝑡) + 75.256 ∙ 𝑦1(𝑡) = 0.299 ∙ sin(90 ∙ 𝑡)

𝑦̈ 2(𝑡) + 3.6975 ∙ 𝑦̇ 1(𝑡) + 1367.15 ∙ 𝑦1(𝑡) = 0.155 ∙ sin(90 ∙ 𝑡)

𝑦̈ 3(𝑡) + 7.9762 ∙ 𝑦̇ 1(𝑡) + 6361.98 ∙ 𝑦1(𝑡) = 0.052 ∙ sin(90 ∙ 𝑡)

𝑦̈4(𝑡) + 10.8071 ∙ 𝑦̇1(𝑡) + 11679.34 ∙ 𝑦1(𝑡) = 0.0082 ∙ sin(90 ∙ 𝑡)

Resolviendo las ecuaciones obtenemos: