6) Una viga simplemente apoyada se ha idealizado como el modelo de masas concentradas

mostrado en la figura 2. Se pide:

a) Identificar los grados de libertad que representen las propiedades elásticas e

inerciales y ensamblar las matrices de rigidez y masa. Desprecie las deformaciones
por fuerza axial.

b) Formular las ecuaciones del movimiento de la viga. Aplicar condensación estática

para reducir el tamaño del problema.

c) Calcular los modos y frecuencias naturales de vibración. La viga tiene una sección de
250x550mm, E=20GPa, m=0,75 tonf/m y L=[número de letras de su segundo apellido
segundo integrante] m.

DATOS:

Viga: 0.25m x 0.55m x 6.00m

L = QUISPE =6
Elasticidad del material: E=20GPa
m=0.75Ton/m

SOLUCIÓN

1. CALCULAMOS LA INERCIA

𝒃𝒉𝟑
𝑰=
𝟏𝟐

0.25 ∗ 0.553
𝐼= = 3.4661 𝑥 10−3
12
2. IDENTIFICACIÓN DE LOS GRADOS DE LIBERTAD

3. ENSAMBLAJE DE LA MATRIZ DE MASA

𝑚𝐿
0 0 0
3
𝑚𝐿
𝒎= 0 0 0
3
0 0 0 0
[ 0 0 0 0]

0.75 𝑥 6
0 0 0
3 1.50 0 0 0
0.75 𝑥 6 0 1.50 0 0
𝒎= 0 0 0 =[ 0 ]
3 0 0 0
0 0 0 0 0 0 0 0
[ 0 0 0 0 ]
4. ENSAMBLAJE DE LA MATRIZ DE RIGIDEZ

81 𝑥 𝐸𝐼 324 𝑥 𝐸𝐼 𝟒𝟎𝟓𝑬𝑰
𝑲𝟏𝟏 = + =
𝐿3 𝐿3 𝑳𝟑
−𝟑𝟐𝟒𝑬𝑰
𝑲𝟐𝟏 =
𝑳𝟑

54𝑥 𝐸𝐼 27 𝑥 𝐸𝐼 𝟐𝟕𝑬𝑰
𝑲𝟑𝟏 = − = 𝟐
𝐿2 𝐿2 𝑳

𝟓𝟒𝑬𝑰
𝑲𝟒𝟏 =
𝑳𝟐

Para obtener la segunda columna de la matriz de rigidez se impone un

desplazamiento unitario al 2do. grado de libertad.
 −𝟑𝟐𝟒𝑬𝑰
𝑲𝟏𝟐 =
𝑳𝟑
81 𝑥 𝐸𝐼 324 𝑥 𝐸𝐼 𝟒𝟎𝟓𝑬𝑰
𝑲𝟐𝟐 = + =
𝐿3 𝐿3 𝑳𝟑
−𝟓𝟒𝑬𝑰
𝑲𝟑𝟐 =
𝑳𝟐
27 𝑥 𝐸𝐼 54 𝑥 𝐸𝐼 −𝟐𝟕𝑬𝑰
𝑲𝟑𝟏 = − =
𝐿2 𝐿2 𝑳𝟐

Para obtener la tercera columna de la matriz de rigidez se impone un giro unitario

al en el 3er grado de libertad.
 54𝑥 𝐸𝐼 27 𝑥 𝐸𝐼 𝟐𝟕𝑬𝑰
𝑲𝟏𝟑 = − = 𝟐
𝐿2 𝐿2 𝑳

−𝟓𝟒𝑬𝑰
𝑲𝟐𝟑 =
𝑳𝟐
9 𝑥 𝐸𝐼 12 𝑥 𝐸𝐼 𝟐𝟏𝑬𝑰
𝑲𝟑𝟑 = + =
𝐿 𝐿 𝑳
𝟔𝑬𝑰
𝑲𝟒𝟑 =
𝑳

Para obtener la cuarta columna de la matriz de rigidez se impone un giro unitario

al en el 4to grado de libertad.
 −𝟓𝟒𝑬𝑰
𝑲𝟐𝟑 =
𝑳𝟐
27𝑥 𝐸𝐼 54 𝑥 𝐸𝐼 −𝟐𝟕𝑬𝑰
𝑲𝟏𝟑 = − =
𝐿2 𝐿2 𝑳𝟐
𝟔𝑬𝑰
𝑲𝟒𝟑 =
𝑳
9 𝑥 𝐸𝐼 12 𝑥 𝐸𝐼 𝟐𝟏𝑬𝑰
𝑲𝟑𝟑 = + =
𝐿 𝐿 𝑳
Finalmente, la matriz de rigidez queda de la siguiente manera:

𝑲𝟏𝟏 𝑲𝟏𝟐 𝑲𝟏𝟑 𝑲𝟏𝟒

𝑲 𝑲𝟐𝟐 𝑲𝟐𝟑 𝑲𝟐𝟒
𝑲 = [ 𝟐𝟏 ]
𝑲𝟑𝟏 𝑲𝟑𝟐 𝑲𝟑𝟑 𝑲𝟑𝟒
𝑲𝟒𝟏 𝑲𝟒𝟐 𝑲𝟒𝟑 𝑲𝟒𝟒
 405𝐸𝐼 −324𝐸𝐼 27𝐸𝐼 54𝐸𝐼
𝐿3 𝐿3 𝐿2 𝐿2
−324𝐸𝐼 405𝐸𝐼 −54𝐸𝐼 −27𝐸𝐼
𝑲= 𝐿3 𝐿3 𝐿2 𝐿2
27𝐸𝐼 −54𝐸𝐼 21𝐸𝐼 6𝐸𝐼
𝐿2 𝐿2 𝐿 𝐿
54𝐸𝐼 −27𝐸𝐼 6𝐸𝐼 21𝐸𝐼
[ 𝐿2 𝐿2 𝐿 𝐿 ]

135 −108 9𝐿 18𝐿

3𝐸𝐼 −108 135 −18𝐿 −9𝐿
𝐾= 3 =[ ]
𝐿 9𝐿 −18𝐿 7𝐿2 2𝐿2
18𝐿 −9𝐿 2𝐿2 7𝐿2

135 −108 9(6) 18(6)

6 −3
3(20 𝑥 10 )(3.466 𝑥 10 ) −108 135 −18(6) −9(6)
𝐾= =
(6) 3 9(6) −18(6) 7(6)2 2(6)2
[18(6) −9(6) 2(6)2 7(6)2 ]

𝟏𝟑𝟓 −𝟏𝟎𝟖 𝟓𝟒 𝟏𝟎𝟖

𝟖𝟔𝟔𝟓 −𝟏𝟎𝟖 𝟏𝟑𝟓 −𝟏𝟎𝟖 −𝟓𝟒
𝑲= =[ ]
𝟗 𝟓𝟒 −𝟏𝟎𝟖 𝟐𝟓𝟐 𝟕𝟐
𝟏𝟎𝟖 −𝟓𝟒 𝟕𝟐 𝟐𝟓𝟐

ENSAMBLAJE DEL VECTOR DE CARGAS EXTERNAS

𝑃1 (𝑡) 0 0 0
0 𝑃2 (𝑡) 0 0
𝑷(𝒕) = [ ] KN
0 0 0 0
0 0 0 0

ECUACIÓN MATRICIAL DEL MOVIMIENTO

̈
𝒎 𝒙 𝒖̈ (𝒕) + 𝒌 𝒙 𝒖(𝒕) = 𝑷(𝒕)
1.10 0 0 0 𝑢̈ 1 (𝑡) 135 −108 54 108 𝑢1 (𝑡) 𝑃1 (𝑡)
0 1.10 0 0 𝑢̈ (𝑡) 8665 −108 135 −108 −54 𝑢 (𝑡) 𝑃 (𝑡)
[ ][ 2 ]+ [ ][ 2 ] = [ 2 ]
0 0 0 0 𝑢̈ 3 (𝑡) 9 54 −108 252 72 𝑢3 (𝑡) 0
0 0 0 0 𝑢̈ 4 (𝑡) 108 −54 72 252 𝑢4 (𝑡) 0

CONDENSACIÓN ESTÁTICA

Se tienen grados de libertad rotacionales sin masa ni carga aplicada, por lo tanto, se
puede dividir las matrices de masa y rigidez considerando solo grados de libertad
traslacionales de la siguiente manera.

𝟏. 𝟏𝟎 𝟎 𝒌 𝒌𝒕𝒐
𝒎=[ ] 𝒌 = [ 𝒕𝒕 ]
𝟎 𝟏. 𝟏𝟎 𝒌𝒐𝒕 𝒌∞
Donde la nueva matriz de rigidez condensa los grados de libertad rotacionales

135 −108 9𝐿 18𝐿

𝑘 𝑘𝑡𝑜 3𝐸𝐼 −108 135 −18𝐿 −9𝐿
𝒌 = [ 𝑡𝑡 ]= 3 [ ]
𝑘𝑜𝑡 𝑘∞ 𝐿 9𝐿 −18𝐿 7𝐿2 2𝐿2
18𝐿 −9𝐿 2𝐿2 7𝐿2

𝑲 = 𝒌𝒕𝒕 − 𝒌𝒕𝒐 𝒙𝒌∞ −𝟏 𝒙𝒌𝒐𝒕

3𝐸𝐼 2 −1
𝒌= ([
135 −108
]−[
9𝐿 18𝐿
] 𝑥 ([7𝐿2 2𝐿2 ] 𝑥[
9𝐿 −18𝐿
]))
𝐿3 −108 135 −18𝐿 −9𝐿 2𝐿 7𝐿2 18𝐿 −9𝐿

3 12
3𝐸𝐼 −
135 −108 9𝐿 18𝐿
𝒌 = 3 ([ ]−[ ] 𝑥 [5𝐿 5𝐿])
𝐿 −108 135 −18𝐿 −9𝐿 12 3
−
5𝐿 5𝐿
432 378
3𝐸𝐼 −
𝒌 = 3 ([ 5 5 ])
𝐿 378 432
−
5 5
162𝐸𝐼 8 −7
𝒌= 𝒙[ ]
5𝐿3 −7 8
162(2 𝑥 106 )(3.466 𝑥 10−3 ) 8 −7
𝒌= 𝒙[ ]
5(6)3 −7 8
5199 8 −7
𝒌= 𝒙[ ]
5 −7 8
Finalmente, la ecuación del movimiento sin los grados de libertad rotacionales es:
̈
𝒎 𝒙 𝒖̈ (𝒕) + 𝒌 𝒙 𝒖(𝒕) = 𝑷(𝒕)
1.10 0 𝑢̈ (𝑡) 5199 8 −7 𝑢1 (𝑡) 𝑃 (𝑡)
[ ][ 1 ]+ [ ][ ]=[ 1 ]
0 1.10 𝑢̈ 2 (𝑡) 5 −7 8 𝑢2 (𝑡) 𝑃2 (𝑡)
1.10 0 𝑢̈ 1 (𝑡) 8318.4 −7275.8 𝑢1 (𝑡) 𝑃 (𝑡)
[ ][ ]+[ ][ ]=[ 1 ]
0 1.10 𝑢̈ 2 (𝑡) −7275.8 8318.4 𝑢2 (𝑡) 𝑃2 (𝑡)

FRECUENCIAS NATURALES DE VIBRACIÓN

|𝒌 − 𝝎𝒏 𝟐 𝒙 𝒎| = 𝟎
Sea 𝝀 = 𝝎𝒏𝟐
8318.4 −7275.8 1.50 0
|[ ]−𝝀 [ ]| = 𝟎
−7275.8 8318.4 0 1.50

8318.4 −7275.8 1.50𝜆 0

|[ ]−[ ]| = 𝟎
−7275.8 8318.4 0 1.50𝜆

𝟖𝟑𝟏𝟖. 𝟒 − 𝟏. 𝟓𝟎𝝀 −𝟕𝟐𝟕𝟓. 𝟖

|[ ]|
−𝟕𝟐𝟕𝟓. 𝟖 𝟖𝟑𝟏𝟖. 𝟒 − 𝟏. 𝟓𝟎𝝀

La determinante del sistema nos da la siguiente ecuación cuadrática

2.25𝜆2 − 24955.2𝜆 + 16258512.92 = 0
𝝀𝟏 = 𝟔𝟗𝟓. 𝟎𝟔𝟕
 𝝀𝟐 = 𝟏𝟎𝟑𝟗𝟔. 𝟏𝟑𝟑

5. CALCULAMOS LA FRECUENCIA CIRCULAR DE VIBRACIÓN ES:

𝝎𝒏𝟏 = √𝝀𝟏 = √695.067 = 𝟐𝟔. 𝟑𝟔𝟒 𝒓𝒂𝒅/𝒔𝒆𝒈

𝝎𝒏𝟐 = √𝝀𝟐 = √10396.133 = 𝟏𝟎𝟏. 𝟗𝟔𝟏 𝒓𝒂𝒅/𝒔𝒆𝒈

6. PERIODO NATURAL DE VIBRACIÓN

𝟐𝝅 2𝜋
𝑻𝒏𝟏 = = = 𝟎. 𝟐𝟑𝟖 𝒓𝒂𝒅/𝒔𝒆𝒈
𝝎𝒏𝟏 26.364
𝟐𝝅 2𝜋
𝑻𝒏𝟐 = = = 𝟎. 𝟎𝟔𝟐 𝒓𝒂𝒅/𝒔𝒆𝒈
𝝎𝒏𝟏 101.961

7. FRECUENCIA NATURAL DE VIBRACIÓN

𝟏 1
𝑭𝒏𝟏 = = = 𝟒. 𝟐𝟎𝟐 𝑯𝒁
𝑻𝒏𝟏 0.238
𝟏 1
𝑭𝒏𝟐 = = = 𝟏𝟔. 𝟏𝟐𝟗 𝑯𝒁
𝑻𝒏𝟐 0.062

8. CALCULAMOS PARA EL MODO 1

Para ello usaremos los siguientes datos: 𝝎𝒏𝟏 = 𝟐𝟔. 𝟑𝟔𝟒 𝒓𝒂𝒅/𝒔𝒆𝒈 y 𝝀𝟏 = 𝟔𝟗𝟓. 𝟎𝟔𝟕

|𝒌 − 𝝎𝒏 𝟐 𝒙 𝒎| 𝒙 𝝓𝟏 = 𝟎
8318.4 −7275.8 1.50𝜆 0 𝜙
|[ ] − [ 0 1 1.50𝜆 ]|)𝑥 [ 𝜙11 ]
−7275.8 8318.4 1 21
8318.4 − 1.50(695.067) −7275.8 𝜙
[ ] 𝑥 [ 𝜙11 ]
−7275.8 8318.4 − 1.50(695.067) 21
𝟕𝟐𝟕𝟓. 𝟖 −𝟕𝟐𝟕𝟓. 𝟖 𝝓𝟏𝟏
[ ]𝒙 [𝝓 ]
−𝟕𝟐𝟕𝟓. 𝟖 𝟕𝟐𝟕𝟓. 𝟖 𝟐𝟏
𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟏𝟏 −𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟐𝟏 = 𝟎
{ }
−𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟏𝟏 𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟐𝟏 = 𝟎
(A)
REMPLAZAMOS EN LA ECUACIÓN A ASUMIENDO 𝝓𝟏𝟏 = 𝟏
𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟏𝟏 − 𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟐𝟏 = 𝟎 (B)
7275.8 𝑥 (1) − 7275.8 𝑥 𝜙21 = 0
𝝓𝟐𝟏 = 𝟏

CALCULAMOS PARA EL MODO 2

Para ello usaremos los siguientes datos: 𝝎𝒏𝟐 = 𝟏𝟎𝟏. 𝟗𝟔𝟏 𝒓𝒂𝒅/𝒔𝒆𝒈 y 𝝀𝟐 = 𝟏𝟎𝟑𝟗𝟔. 𝟏𝟑𝟑
|𝒌 − 𝝎𝒏 𝟐 𝒙 𝒎| 𝒙 𝝓𝟐 = 𝟎
8318.4 −7275.8 1.50𝜆 0 𝜙
|[ ] − [ 0 2 1.50𝜆 ]|)𝑥 [ 𝜙11 ]
−7275.8 8318.4 2 21
8318.4 − 1.50(10396.133) −7275.8 𝜙
[ ] 𝑥 [ 𝜙11 ]
−7275.8 8318.4 − 1.50(10396.133) 21
 −𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟏𝟐 −𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟐𝟐 = 𝟎
{ }
−𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟏𝟐 −𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟐𝟐 = 𝟎

REMPLAZAMOS EN LA ECUACIÓN A ASUMIENDO 𝝓𝟏𝟐 = 𝟏

-𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟏𝟏 − 𝟕𝟐𝟕𝟓. 𝟖 𝒙 𝝓𝟐𝟏 = 𝟎

7275.8 𝑥 (1) − 7275.8 𝑥 𝜙21 = 0
𝝓𝟐𝟏 = −𝟏

MODO1

MODO2