Tema: Dinámica de Estructuras

Tópico: Sistemas de "n" grados de Libertad

Ing. Eliud Hernández

PÓRTICO IDEALIZADO DE 2 PISOS

(MODOS DE VIBRACIÓN)

Masa por Nivel:

m2 ≔ 19.62 tonne m1 ≔ 2 ⋅ m2 = 39.24 tonne

Rigidez por Nivel:

tonnef tonnef
k2 ≔ 15 ――― k1 ≔ 2 ⋅ k2 = 30 ―――
cm cm

1. Matriz de Masa

⎡ m 0 ⎤ ⎡ 39.24 0 ⎤
m≔⎢ 1 ⎥=⎢ ⎥ tonne
⎣ 0 m2 ⎦ ⎣ 0 19.62 ⎦

2. Matriz de Rigidez

⎡ k + k2 −k2 ⎤ ⎡ 45 −15 ⎤ tonnef

k≔⎢ 1 ⎥= ―――
⎣ −k2 k2 ⎦ ⎢⎣ −15 15 ⎥⎦ cm

3. Cálculo de Autovalores

Se plantea el determinante

(k − λ ⋅ m) ＝ 0

⎡ k1 + k2 −k2 ⎤ ⎡ m1 0 ⎤ ⎡ ⎛⎝k1 + k2⎞⎠ − λ ⋅ m1 −k2 ⎤

⎢ ⎥−λ⋅⎢ ⎥＝0 ⎢ ⎥＝0
⎣ −k2 k2 ⎦ ⎣ 0 m2 ⎦ ⎣ −k2 k2 − λ ⋅ m2 ⎦

⎛⎝⎛⎝k1 + k2⎞⎠ − λ ⋅ m1⎞⎠ ⋅ ⎛⎝k2 − λ ⋅ m2⎞⎠ − ⎛⎝−k2⎞⎠ ⋅ ⎛⎝−k2⎞⎠ ＝ 0

⎛⎝k1 + k2⎞⎠ ⋅ k2 − ⎛⎝k1 + k2⎞⎠ ⋅ ⎛⎝λ ⋅ m2⎞⎠ − λ ⋅ m1 ⋅ k2 + ⎛⎝λ ⋅ m1⎞⎠ ⋅ ⎛⎝λ ⋅ m2⎞⎠ − k2 2 ＝ 0

INESA adiestramiento c.a. www.inesa-adiestramiento.com

Derechos Reservados: No se permite hacer un uso comercial de este trabajo sin autorización previa del autor. La misma deberá estar firmada.
Tema: Dinámica de Estructuras
Tópico: Sistemas de "n" grados de Libertad
Ing. Eliud Hernández

λ 2 ⋅ ⎛⎝m1 ⋅ m2⎞⎠ − λ ⋅ ⎛⎝⎛⎝k1 + k2⎞⎠ ⋅ m2 + m1 ⋅ k2⎞⎠ + ⎛⎝k1 + k2⎞⎠ ⋅ k2 − k2 2 ＝ 0

λ 2 ⋅ ⎛⎝m1 ⋅ m2⎞⎠ − λ ⋅ ⎛⎝⎛⎝k1 + k2⎞⎠ ⋅ m2 + m1 ⋅ k2⎞⎠ + k1 ⋅ k2 + k2 2 − k2 2 ＝ 0

λ 2 ⋅ ⎛⎝m1 ⋅ m2⎞⎠ − λ ⋅ ⎛⎝⎛⎝k1 + k2⎞⎠ ⋅ m2 + m1 ⋅ k2⎞⎠ + k1 ⋅ k2 ＝ 0

tonne 2
a ≔ m1 ⋅ m2 = 769.889 tonne 2 c ≔ k1 ⋅ k2 = ⎛⎝4.328 ⋅ 10 8 ⎞⎠ ―――
s4

tonne 2
b ≔ −⎛⎝⎛⎝k1 + k2⎞⎠ ⋅ m2 + m1 ⋅ k2⎞⎠ = −1.443 ⋅ 10 6 ―――
s2

λ2 ⋅ a + λ ⋅ b + c＝0

−b − ‾‾‾‾‾‾‾‾‾
b2 − 4 ⋅ a ⋅ c 1 −b + ‾‾‾‾‾‾‾‾‾
b2 − 4 ⋅ a ⋅ c ⎛ 1
λ1 ≔ ――――――= 374.872 ― λ2 ≔ ――――――= ⎝1.499 ⋅ 10 3 ⎞⎠ ―
2⋅a s2 2⋅a s2

Alternativamente, para cada modo, se puede obtener el autovalor a través de una función de
iteración en Mathcad. Para el modo 1 se fija un valor inicial de λ =1 y la función lo va variando
hasta encontrar el valor final λ que cumple la condición de que el determinante sea igual a cero.
Luego, para el modo 2, se define un valor inicial de λ mayor al obtenido para el modo 1 a fin de
realizar nuevamente la iteración.

Modo 1 Modo 2

1 1
λ≔1⋅― λ ≔ 1100 ⋅ ―
s2 s2

det (k − λ ⋅ m) ＝ 0 det (k − λ ⋅ m) ＝ 0

1 1
λ1 ≔ find (λ) = 374.872 ― λ2 ≔ find (λ) = 1499.488 ―
s2 s2

INESA adiestramiento c.a. www.inesa-adiestramiento.com

Derechos Reservados: No se permite hacer un uso comercial de este trabajo sin autorización previa del autor. La misma deberá estar firmada.
Tema: Dinámica de Estructuras
Tópico: Sistemas de "n" grados de Libertad
Ing. Eliud Hernández

4. Frecuencias y Períodos
1
ω1 ≔ ‾‾
λ1 = 19.362 ― Frecuencia del Modo 1
s

2π
T1 ≔ ―― = 0.325 s Período del Modo 1
ω1

1
ω2 ≔ ‾‾
λ2 = 38.723 ― Frecuencia del Modo 2
s

2π
T2 ≔ ―― = 0.162 s Período del Modo 2
ω2

5. Cáculo de Autovectores

ϕ se define como el Autovector.

Análisis del Modo 1 Análisis del Modo 2

1 1
para λ1 = 374.872 ― para λ2 = ⎛⎝1.499 ⋅ 10 3 ⎞⎠ ―
s2 s2

⎡ ⎛⎝k1 + k2⎞⎠ − λ1 ⋅ m1 −k2 ⎤ ⎡ ϕ11 ⎤ ⎡ ⎛⎝k1 + k2⎞⎠ − λ2 ⋅ m1 −k2 ⎤ ⎡ ϕ12 ⎤

⎢ ⎥⋅⎢ ⎥＝0 ⎢ ⎥⋅⎢ ⎥＝0
⎣ −k2 k2 − λ1 ⋅ m2 ⎦ ⎣ ϕ21 ⎦ ⎣ −k2 k2 − λ2 ⋅ m2 ⎦ ⎣ ϕ22 ⎦

Resolviendo con la primera fila se tiene: Resolviendo con la primera fila se tiene:

⎛⎝⎛⎝k1 + k2⎞⎠ − λ1 ⋅ m1⎞⎠ ⋅ ϕ11 − k2 ⋅ ϕ21 ＝ 0 ⎛⎝⎛⎝k1 + k2⎞⎠ − λ2 ⋅ m1⎞⎠ ⋅ ϕ12 − k2 ⋅ ϕ22 ＝ 0

si ϕ21 ≔ 1 si ϕ22 ≔ 1

k2 ⋅ ϕ21 k2 ⋅ ϕ22
ϕ11 ≔ ――――――= 0.5 ϕ12 ≔ ――――――= −1
⎛⎝⎛⎝k1 + k2⎞⎠ − λ1 ⋅ m1⎞⎠ ⎛⎝⎛⎝k1 + k2⎞⎠ − λ2 ⋅ m1⎞⎠

INESA adiestramiento c.a. www.inesa-adiestramiento.com

Derechos Reservados: No se permite hacer un uso comercial de este trabajo sin autorización previa del autor. La misma deberá estar firmada.
Tema: Dinámica de Estructuras
Tópico: Sistemas de "n" grados de Libertad
Ing. Eliud Hernández

Alternativamente, para cada modo, se pueden obtener los autovectores a través de una función de
iteración en Mathcad. Se fija un valor inicial igual a 1.00 para ambas componentes del autovector
correspondiente y la función va variando dichos valores hasta que se satisface la ecuación.

Modo 1 Modo 2

ϕ11 ≔ 1 ϕ12 ≔ 1

ϕ21 ≔ 1 ϕ22 ≔ 1

⎡ ϕ ⎤ ⎡ ⎤ tonnef ⎡ ϕ ⎤ ⎡ ⎤ tonnef
⎛⎝k − λ1 ⋅ m⎞⎠ ⋅ ⎢ 11 ⎥ ＝ ⎢ 0 ⎥ ⋅ ――― ⎛⎝k − λ2 ⋅ m⎞⎠ ⋅ ⎢ 12 ⎥ ＝ ⎢ 0 ⎥ ⋅ ―――
⎣ 1 ⎦ ⎣0⎦ cm ⎣ 1 ⎦ ⎣0⎦ cm

⎡ 0.5 ⎤ ⎡ −1 ⎤
ϕ1 ≔ find ⎛⎝ϕ11 , ϕ21⎞⎠ = ⎢ ϕ2 ≔ find ⎛⎝ϕ12 , ϕ22⎞⎠ = ⎢
⎣ 1 ⎥⎦ ⎣ 1 ⎥⎦

6. Normalización de los modos

Modo 1 Modo 2

1 ⎡ 0.092 ⎤ 1 1 ⎡ −0.13 ⎤ 1
φ1 ≔ ――――― ⋅ ϕ1 = ⎢ ⎥ ――― φ2 ≔ ――――― ⋅ ϕ2 = ⎢ ⎥ ―――
‾‾‾‾‾‾‾‾‾
ϕ1 T ⋅ m ⋅ ϕ1 ⎣ 0.184 ⎦ 1
― ‾‾‾‾‾‾‾‾‾
ϕ2 T ⋅ m ⋅ ϕ2 ⎣ 0.13 ⎦ 1
―
2 2
tonne tonne

7. Ortogonalidad de los modos

ϕ1 T ⋅ m ⋅ ϕ2 = 0 tonne Se comprueba la ortogonalidad

tonnef
ϕ1 T ⋅ k ⋅ ϕ2 = 0 ――― Se comprueba la ortogonalidad
cm

INESA adiestramiento c.a. www.inesa-adiestramiento.com

Derechos Reservados: No se permite hacer un uso comercial de este trabajo sin autorización previa del autor. La misma deberá estar firmada.
Tema: Dinámica de Estructuras
Tópico: Sistemas de "n" grados de Libertad
Ing. Eliud Hernández

8. Formas Modales

Modo 1 Modo 2

ϕ21 ϕ22

ϕ11 ϕ12

⎡ϕ ⎤ ⎡ϕ ⎤
ϕ1 ＝ ⎢ 11 ⎥ ϕ2 ＝ ⎢ 12 ⎥
⎣ ϕ21 ⎦ ⎣ ϕ2 ⎦

⎡ 0.5 ⎤ ⎡ −1 ⎤
ϕ1 = ⎢ ϕ2 = ⎢
⎣ 1 ⎥⎦ ⎣ 1 ⎥⎦

INESA adiestramiento c.a. www.inesa-adiestramiento.com

Derechos Reservados: No se permite hacer un uso comercial de este trabajo sin autorización previa del autor. La misma deberá estar firmada.