Capitulo 2
Capitulo 2
Capitulo 2
DINÁMICA
INGENIERIA CIVIL
2018
INGENIERIA CIVIL 1
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Contenido
....................................................................................................................................................... 3
CAPITULO I: ................................................................................................................................... 4
1. PREGUNTAS DE REPASO .................................................................................................... 4
2. PROBLEMAS........................................................................................................................... 6
INGENIERIA CIVIL 2
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
CAPITULO 1:
FUNDAMENTOS DE
VIBRACIÓN
INGENIERIA CIVIL 3
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
CAPITULO I:
1. PREGUNTAS DE REPASO
INGENIERIA CIVIL 4
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
3. La frecuencia con la cual vibra un sistema inicialmente perturbado por
sí mismo se conoce como la frecuencia natural.
4. Cualquier función periódica puede expandirse a una serie de Fourier.
5. Un movimiento armónico es un movimiento periódico.
6. La masa equivalente de varias masas en diferentes lugares se encuentra
utilizando la equivalencia de energía cinética.
7. Las coordenadas generalizadas no necesariamente son coordenadas
cartesianas.
8. Los sistemas discretos son los mismos que los sistemas de parámetro
concentrado.
9. Considere la suma de movimientos armónicos,
10. Considere la suma de movimientos armónicos
SOLUCION:
INGENIERIA CIVIL 5
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
2. PROBLEMAS
EJERCICIO N° 1.6: Desarrolle un modelo matemático para el tractor y el arado de
la figura 1.66 considerando la masa, elasticidad y amortiguamiento de las llantas,
los amortiguadores y el arado (las hojas).
SOLUCIÓN:
INGENIERIA CIVIL 6
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
SOLUCIÓN:
Desarrollamos:
1 1 1 1
= + +
𝑘𝑒𝑞 2𝑘𝑒𝑞 𝑘2 3𝑘3
1 2𝑘1 𝑘2 𝑘3
=
𝑘𝑒𝑞 𝑘2 𝑘3 + 2𝑘1 𝑘3 + 𝑘1 𝑘2
2𝑘1 𝑘2 𝑘3
( + 𝑘4 )(𝑘5 )
𝑘2 𝑘3 + 2𝑘1 𝑘3 + 𝑘1 𝑘2
𝑘𝑒𝑞𝑡𝑜𝑡𝑎𝑙 =
2𝑘1 𝑘2 𝑘3
+ 𝑘4 + 𝑘5
𝑘2 𝑘3 + 2𝑘1 𝑘3 + 𝑘1 𝑘2
2𝑘1 𝑘2 𝑘3 + 𝑘2 𝑘3 𝑘4 + 2𝑘1 𝑘3 𝑘4 + 𝑘1 𝑘2 𝑘4
( )𝑘5
𝑘2 𝑘3 + 2𝑘1 𝑘3 + 𝑘1 𝑘2
𝑘𝑒𝑞𝑡𝑜𝑡𝑎𝑙 =
2𝑘1 𝑘2 𝑘3 + 𝑘2 𝑘3 𝑘4 + 2𝑘1 𝑘3 𝑘4 + 𝑘1 𝑘2 𝑘4 + 𝑘2 𝑘3 𝑘5 + 2𝑘1 𝑘3 𝑘5 + 𝑘1 𝑘2 𝑘5
𝑘2 𝑘3 + 2𝑘1 𝑘3 + 𝑘1 𝑘2
2𝑘1 𝑘2 𝑘3 + 𝑘2 𝑘3 𝑘4 + 2𝑘1 𝑘3 𝑘4 + 𝑘1 𝑘2 𝑘4
𝑘𝑒𝑞𝑡𝑜𝑡𝑎𝑙 =
2𝑘1 𝑘2 𝑘3 + 𝑘2 𝑘3 𝑘4 + 2𝑘1 𝑘3 𝑘4 + 𝑘1 𝑘2 𝑘4 + 𝑘2 𝑘3 𝑘5 + 2𝑘1 𝑘3 𝑘5 + 𝑘1 𝑘2 𝑘5
INGENIERIA CIVIL 7
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
SOLUICIÓN:
𝑥1 +𝑥2
∑ 𝐹𝑦 = 0 → 𝑥1 𝑘1 + 𝑥2 𝑘2 = + 𝑘𝑒𝑞
2
∑𝑀 = 0 → 𝑥1 𝑙1 𝑘1 = 𝑥2 𝑙2 𝑘2
𝑥1 𝑙2 𝑘2
=
𝑥2 𝑙1 𝑘1
2𝑥1 𝑘1 + 2𝑥2 𝑘2
= 𝑘𝑒𝑞 𝑙𝑢𝑒𝑔𝑜 𝑑𝑖𝑣𝑖𝑑𝑖𝑚𝑜𝑠 𝑒𝑛𝑡𝑟𝑒 𝑥1 𝑦 𝑥2
𝑥1 + 𝑥2
𝑥 𝑙1 𝑘2
2 𝑥1 + 2𝑘2 2( ) (𝑘1 ) + 2𝑘2
2 𝑙1 𝑘1
𝑘𝑒𝑞 = 𝑥1 =
𝑙2 𝑘2
𝑥2 + 1 𝑙1 𝑘1
+1
𝑙1
2 𝑘 + 2𝑘2
𝑙2 2
𝑘𝑒𝑞 =
𝑙2 𝑘2
𝑙1 𝑘1
INGENIERIA CIVIL 8
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
2𝑘1 𝑘2 (𝑙1 + 𝑙2 )
𝑘𝑒𝑞 =
𝑙1 𝑘1 + 𝑙2 𝑘2
solucion:
Elemento tensorial:
𝑀 𝐺𝐽
𝑀 = 𝑘𝑡 (𝜃𝐴 + 𝜃𝐵 ) 𝑘𝑡 = 𝑘𝑡 =
|𝜃𝐴 − 𝜃𝐵 | 𝐿
De donde:
𝑀: 𝑚𝑜𝑚𝑒𝑛𝑡𝑜 𝐺: 𝑀𝑜𝑑𝑢𝑙𝑜 𝑑𝑒 𝑟𝑖𝑔𝑖𝑑𝑒𝑧 𝑐𝑜𝑟𝑡𝑎𝑛𝑡𝑒
𝜃𝐴 , 𝜃𝐵 : 𝐴𝑛𝑔𝑢𝑙𝑜𝑠 𝑡𝑜𝑟𝑠𝑖𝑜𝑛𝑎𝑙𝑒𝑠 L: longitud del elemento
𝑘𝑡 : 𝑅𝑖𝑔𝑖𝑑𝑒𝑧 𝑡𝑜𝑟𝑠𝑜𝑟𝑖𝑎𝑙 J: Momento polar de inercia
Por lo tanto resolviendo:
𝑀1 = 𝐾𝑡1 𝜃
𝑀1 = 𝐾𝑡2 𝜃
De la figura tenemos:
𝑘𝑒𝑞 𝜃 = 𝑀1 + 𝑀2 + 𝑀3
𝑘𝑒𝑞 𝜃 = 𝜃𝐾𝑡1 + 𝜃𝐾𝑡2 + (𝜃𝑙1 𝑘1 + 𝜃𝑙2 𝑘3 − 𝜃𝑙1 𝑘2 )
Factorizando y luego simplificando 𝜃 tenemos:
𝑘𝑒𝑞 = 𝐾𝑡1 + 𝐾𝑡2 + (𝑙1 𝑘1 + 𝑙2 𝑘3 − 𝑙1 𝑘2 )
INGENIERIA CIVIL 9
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Solucion:
𝐺𝑖 𝐼𝑖
𝑘𝑡𝑖 = 𝑘𝑖 = 𝑘1 + 𝑘2 + 𝑘3 …
𝐿𝑖
Simplificando 𝜃 tenemos:
𝑘𝑒𝑞 = 𝑘1 + 𝑘2 + 𝑘3 + 𝑘4 + 𝑘5 − 𝑘6
EJERCICIO N° 1.11: Una máquina de masa m= 5 500 kg está montada en una viga
de acero sólo apoyada de longitud l=2m que tiene una sección transversal (de
profundidad=0.1 y ancho=1.2m) y módulo de Young 𝐸 = 2.06𝑥1011 𝑁/𝑚2 . para
reducir la deflexión vertical de la viga, se fija un resorte de rigidez k a la mitad de
su claro, como se muestra en la figura 1.71. Determine el valor de k necesario
para reducir la deflexión de la viga en:
a. 25 por ciento de su valor original.
b. 50 por ciento de su valor original.
c. 75 por ciento de su valor original.
Suponga que la masa de la viga es insignificante.
INGENIERIA CIVIL 10
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
SOLUCION:
Reduciendo el esquema tenemos:
Datos:
𝑏ℎ3
𝐼=
3
3𝐸𝐼(𝑎 + 𝑏)
𝑘𝑣𝑖𝑔𝑎 = 𝑘1 =
𝑎2 𝑏 2
3𝐸𝐼𝐿 48𝐸𝐼
𝑘1 = = 3
𝑙4 𝑙
16
48(2.06𝑥1011 )
𝑘1 = = 4944𝑥105
8
INGENIERIA CIVIL 11
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Si:
𝑘𝑒𝑞 = 𝑘1 − 𝑘
Para producir 25% tenemos:
75
(𝑘1 − 𝑘)𝑥 = 𝑘 𝑥
100 1
75
𝑘1 − 𝑘 =𝑘
100 1
𝑘 = 4944𝑥103
Para producir descuento de 50%:
50 5
𝑘= 𝑘1 = 𝑥4944𝑥105
100 10
𝑘 = 2472𝑥103
Para producir descuento de 75%:
75 25
𝑘= 𝑘1 = 𝑥4944𝑥105
100 100
𝑘 = 1225𝑥102
EJERCICIO N° 1.6 :
EJERCICIO N° 1.6:
EJERCICIO N° 1.6 :
EJERCICIO N° 1.6 :
EJERCICIO N° 1.6 :
EJERCICIO N° 1.6 :
EJERCICIO N° 1.6 :
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
EJERCICIO N° 1.6
INGENIERIA CIVIL 12
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
CAPITULO VI:
2. PREGUNTAS DE REPASO
INGENIERIA CIVIL 13
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
6.2 Indique si cada uno de los siguientes enunciados es verdadero o falso:
1. Para un sistema de varios grados de libertad, se puede escribir una ecuación
de movimiento por cada grado de libertad.
2. La ecuación de Lagrange no se puede utilizar para derivar las ecuaciones de
movimiento de un sistemas de varios grados de libertad.
3. Las matrices de masa, rigidez y amortiguamiento de un sistema de varios
grados de libertad siempre son simétricas.
4. El producto de las matrices de rigidez y flexibilidad de un sistema siempre es
una matriz de identidad.
5. El análisis modal de un sistema de n grados de libertad se puede realizar por
medio de r modos con 𝑟 < 𝑛
6. Para un sistema amortiguado de varios grados de libertad, todos los valores
eigen pueden ser complejos.
7. La relación de amortiguamiento modal indica amortiguamiento en un modo
normal particular.
8. Un sistema de varios grados de libertad puede tener seis de las frecuencias
naturales iguales a cero.
9. Las coordenadas generalizadas siempre tendrán la longitud unitaria.
10. Las coordenadas generalizadas son independientes de las condiciones de
restricción del sistema.
11. La matriz de masa generalizada de un sistema de varios grados de libertad
siempre es diagonal.
12. Las energías potencial y cinética de un sistema de varios grados de libertad
siempre son cantidades cuadráticas.
13. La matriz de masa de un sistema siempre es simétrica y definida positiva.
14. La matriz de rigidez de un sistema siempre es simétrica y definida positiva.
15. El modo de cuerpo rígido también se llama modo cero.
16. Un sistema no restringido también se conoce como sistema semidefinido.
17. La segunda ley del movimiento de Newton siempre se puede utilizar para
derivar las ecuaciones de movimiento de un sistema vibratorio.
1.3 Escriba en cada uno de los siguientes espacios en blanco la palabra correcta.
1. La constante de resorte indica la ___________ necesaria para producir un
alargamiento unitario
2. El coeficiente de influencia de flexibilidad aij indica la deflexión en el punto
__________ debido a una carga unitaria aplicada en el punto ___________.
3. La fuerza en el punto i producida por un desplazamiento unitario en el
punto j, cuando todos los puntos además del punto j están fijos, se conoce
como coeficiente de influencia de ____________.
4. Las formas de modo de un sistema de varios grados de libertad son
__________.
5. Las ecuaciones de movimiento de un sistema de varios grados de libertad se
pueden expresar en función de los coeficientes de __________.
6. Las ecuaciones de Lagrange se expresan en función de coordenadas
____________.
INGENIERIA CIVIL 14
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
7. El valor de la delta de Kronecker (𝛿𝑖𝑗 ) es 1 para todas las i = j y _________para
todas las i ≠ j.
8. La matriz de rigidez de un sistema semidefinido es ___________.
9. Un sistema de varios grados de libertad puede tener a lo sumo __________de
cuerpo rígido.
10. Cuando el vector de solución se indica como una combinación lineal de los
modos normales como 𝑥⃗(𝑡) = ∑𝑛𝑖=1 𝑞𝑖 (𝑡)𝑥⃗ 𝑖 las coordenadas generalizadas
𝑞𝑖 (𝑡) también se conocen como coeficientes de participación de ___________.
11. Cualquier conjunto de n vectores linealmente independientes en un espacio
de n dimensiones se llama __________.
12. La representación de un vector de n dimensiones arbitrario como una
combinación lineal de n vectores linealmente independientes se conoce
como teorema de ___________.
13. El análisis ___________ está basado en el teorema de expansión.
14. El análisis modal básicamente ___________ las ecuaciones de movimiento.
15. Los valores eigen de un sistema de varios grados de libertad forman un
___________ en el espacio de n dimensiones.
16. La aplicación de las ecuaciones de Lagrange requiere la disponibilidad de
___________ expresiones.
17. La ecuación en forma de determinante|[𝑘] − 𝑤 2 [𝑚]| = 0, se conoce como
ecuación ___________.
18. La simetría de las matrices de rigidez y flexibilidad se debe al teorema de
reciprocidad ___________.
19. El teorema de reciprocidad de Maxwell establece que los coeficientes de
influencia son __________.
20. La matriz de rigidez es positiva definida sólo si el sistema es ________.
21. Durante la vibración libre de un sistema no amortiguado, todas las
coordenadas tendrán movimiento ___________.
22. En amortiguamiento proporcional, se supone que la matriz de
amortiguamiento es una combinación lineal de las matrices de _________ y
___________.
INGENIERIA CIVIL 15
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
𝑡
𝑥⃗ (1)
(2)𝑡
b) [𝑥] = 𝑥⃗
⋮
[𝑥⃗ (𝑛)𝑡 ]
c) [𝑥] = [𝑘]−1 [𝑚]
5. La función de disipación de Rayleigh se utiliza para generar una
a) matriz de rigidez
b) matriz de amortiguamiento
c) matriz de masa
6. La ecuación característica de un sistema de n grados de libertad es
a) una ecuación trascendental
b) un polinomio de grado n
c) una ecuación diferencial de orden n
7. La frecuencia natural fundamental de un sistema es
a) el valor máximo
b) el valor mínimo
c) cualquier valor
8. El amortiguamiento negativo conduce a
a) Inestabilidad
b) convergencia rápida
c) oscilaciones
9. El criterio de Routh-Hurwitz se puede utilizar para investigar
a) la convergencia de un sistema
b) las oscilaciones de un sistema
c) la estabilidad de un sistema
10. Las matrices de rigidez y flexibilidad están relacionadas como
a) [k]=[a]
b) [k]=[𝑎]−1
c) [k]=[𝑎]𝑡
11. Un sistema para el cual [k] es positiva y [m] es positiva definida se llama
a) sistema semidefinido
b) sistema positivo definido
c) sistema indefinido
12. [m]-ortogonalidad de vectores modales implica
𝑡
a) 𝑥⃗ (𝑖) [𝑚]𝑥⃗ (𝑖) = 0
𝑡
b) 𝑥⃗ (𝑖) [m]𝑥⃗ (𝑗) = 0
c) [𝑥]𝑡 [𝑚][𝑥] = [𝜔𝑖2 ]
13. Se puede utilizar el análisis modal de manera conveniente para determinar
la respuesta de un sistema de varios grados de libertad.
a) sometido a condiciones forzadas arbitrarias
b) sometido a condición de vibración libre
c) implicando varios modos
INGENIERIA CIVIL 16
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
6.4 𝑐𝑜𝑟𝑟𝑒𝑙𝑎𝑐𝑖𝑜𝑛𝑒 𝑙𝑜𝑠 𝑒𝑙𝑒𝑚𝑒𝑛𝑡𝑜𝑠 𝑑𝑒 𝑙𝑎𝑠 𝑑𝑜𝑠 𝑐𝑜𝑙𝑢𝑚𝑛𝑎𝑠 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒𝑠:
1
1. 𝑥⃗̇ 𝑡 [𝑚]𝑥⃗̇ a. igual a cero produce los valores característicos.
2
1
2. 𝑥⃗ 𝑡 [𝑚]𝑥⃗ b. igual a [𝜔𝑖2 ] cuando los modos son normalizados.
2
𝑡
3. 𝑥⃗ (𝑖) [m]𝑥⃗ (𝑗) c. energía cinética del sistema.
𝑡
4. 𝑥⃗ (𝑖) [𝑚]𝑥⃗ (𝑖) d. igual a cero cuando los modos son ortogonales.
5. [𝑥]𝑡 [𝑘][𝑥] e. igual a la matriz dinámica [D]
6. [𝑚]𝑥⃗̈ + [𝑘]𝑥⃗ f. energía de deformación del sistema
7. |[𝑘] − 𝜔2 [𝑚]| g. igual al vector de fuerza aplicada 𝐹⃗
8. [𝑘]−1 [m] h. igual a uno cuando los modos son ortogonales
Problemas
Solución
INGENIERIA CIVIL 17
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
𝑚3 : − 5𝑘(𝑥3 − 𝑥1 ) − 𝑘(𝑥3 − 𝑥2 ) − 𝑘𝑥3 + 𝐹3 (𝑡) = 𝑚2 𝑥̈ 1
𝑚1 𝑥̈ 1 + 7𝑘𝑥1 − 𝑘𝑥2 − 5𝑘𝑥3 = 𝑓1 (𝑡)
𝑚2 𝑥̈ 2 − 𝑘𝑥1 − 2𝑥2 − 𝑘𝑥3 = 𝐹2 (𝑡)
𝑚3 𝑥̈ 3 − 5𝑘𝑥1 − 1𝑘𝑥2 + 7𝑘𝑥3 = 𝐹3 (𝑡)
𝑚1 0 0 𝑥̈ 1 7 −1 −5 𝑥1 𝐹1 (𝑡)
[ 0 𝑚2 0 ] {𝑥̈ 2 } + 𝑘 [−1 2 −1] {𝑥2 } = {𝐹2 (𝑡)}
0 0 𝑚3 𝑥̈ 3 −5 −1 7 𝑥3 𝐹3 (𝑡)
Ejercicio 6.12
Derive la matriz de rigidez de cada uno de los sistemas que se muestran en las
figuras 6.18 a 6.23 utilizando las coordenadas indicadas.
7𝐾 −𝐾 −5𝐾 7 −1 −5
∴ [K]=[ −𝐾 2𝐾 −𝐾 ]=K[−1 2 −1]
−5𝐾 −𝐾 7𝐾 −5 −1 7
INGENIERIA CIVIL 18
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Ejercicio 6.21
Derive la matriz de rigidez del sistema que se muestra en la figura 5.42.
Solución:
X=11, 𝜃=0
∑ 𝑀 = 0 ⇒ 𝑘𝜃𝑥 − 𝑘𝑙 = 0 ∑ 𝐹 = 0 ⇒ 𝑘𝑥𝑥 = 2𝑘
𝑘𝜃𝑥 = 𝑘𝑙
∑ M=0 ∑ M=0
−𝑘𝑙 2 − 𝑘𝑙 2 + 𝑘𝜃𝜃 = 0 𝑘𝑥𝜃 = 𝑘𝑙
𝑘𝜃𝜃 = +2𝑘𝑙 2
INGENIERIA CIVIL 19
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
2𝑘 𝑘𝑙
∴ [K]=[ ]
𝑘𝑙 𝑘𝑙 2
Ejercicio 6.22
Derive la matriz de masa del sistema que se muestra en la figura 5.42
Solución:
Calculamos la energía cinética del sistema en 6.21.
1 1 1
T= (𝑚1 )𝑣1 2 + (𝑚2 )𝑣2 2 + (𝑚8 )𝑣8 2
2 2 2
1 1 1
T=2 𝑚1 (𝜃̇𝑙)2 + 2 𝑚2 (𝜃̇𝑙)2 + 2 𝑚𝑥̇ 2
1 1 1
T= 2𝑚(𝜃̇𝑙)2 + 2𝑚(𝜃̇𝑙)2 + 𝑚𝑥̇ 2
2 2 2
1 1 1
T=𝑚𝜃 𝑙 + 𝑚𝜃̇ 2 𝑙 2 + 2 𝑚𝑥̇ 2 = 2 (4𝑚)𝜃̇ 2 𝑙 2 + 2 𝑚𝑥̇ 2
−2 2
INGENIERIA CIVIL 20
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Para: 𝑥1 = 1; 𝑥2 = 0; 𝑥3 = 0; 𝑥4 = 0
𝑘11 = 𝑘1 + 𝑘2
𝑘12 = −𝑘2
𝑘13 = 0
Para: 𝑥1 = 0; 𝑥2 = 1; 𝑥3 = ; 𝑥4 = 0
𝑘21 = −𝑘2
𝑘22 = 𝑘2 + 𝑘3
𝑘23 = −𝑘3
Para: 𝑥1 = 0; 𝑥2 = 0; 𝑥3 = 0 ; 𝑥4 = 0
𝑘31 = 0
𝑘32 = −𝑘3
𝑘33 = 𝑘3 + 𝑘4
Matriz de rigidez.
𝑚1 0 𝑜 𝑥̈ 1 𝑘1 + 𝑘2 −𝑘2 𝑜 𝑥1 0
[ 0 𝑚2 𝑜 ] {𝑥̈ 2 } + [ −𝑘2 𝑘2 + 𝑘3 −𝑘3 ] {𝑥2 } = {0}
0 0 𝑚3 𝑥̈ 3 0 −𝑘3 𝑘3 + 𝑘4 𝑥3 0
Coeficiente de influencia de flexibilidad.
[𝑐]=[𝑘]−1
INGENIERIA CIVIL 21
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Ejercicio 6.26.
Determine la matriz de flexibilidad de la barra uniforme que se muestra en la figura 6.31. Ignore la masa
de la viga comparada con las masas concentradas colocadas sobre la viga y suponga que todas las 𝑙𝑖 = l.
Deflexión de viga.
𝐹𝑏 2 𝑥2
𝑦(𝑥)𝐴𝐵 = [−𝑥(3𝑎 + 𝑏) + 3𝑎𝐿]………………………..1
𝐺𝐸𝐼𝐿3
𝐹2 𝐴2 (1−𝑥)2
𝑦(𝐴)𝐵𝐶 = [−(𝑐 − 𝑥)(3𝑎 + 𝑏) + 3𝑏𝐿]………..2
𝐺𝐸𝐼𝐿3
1 𝑙3
𝑎21 =
6 𝐼𝐸
1 𝑙3
𝑎22 =
3 𝐼𝐸
1 𝑙3
𝑎23 =
6 𝐼𝐸
13 𝑙 3
𝑎31 =
192 𝐼𝐸
1 𝑙3
𝑎32 =
6 𝐸𝐼
9 𝑙3
𝑎33 =
64 𝐸𝐼
INGENIERIA CIVIL 22
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Matriz de rigidez:
I. x1 = 0; x2 = 0; x3 = 0
𝑘11 = 𝑘1 + 𝑘2
𝑘12 = −𝑘2
𝑘31 = 0
II. x1 = 0; x2 = 0; x3 = 0
𝑘21 = −𝑘2
𝑘22 = 𝑘2 + 𝑘3
𝑘23 = −𝑘3
INGENIERIA CIVIL 23
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
III. 𝑥1 = 0; 𝑥2 = 0; 𝑥3 = 1
𝑘1 + 𝑘2 −𝑘2 0
[𝑘] = [ −𝑘2 𝑘2 + 𝑘3 −𝑘3 ]
0 −𝑘3 𝑘3
Matriz de flexibilidad:
1 1 1
𝑘1 𝑘1 𝑘1
1 1 1 1 1
[𝑎] = [𝑘]−1 = + +
𝑘1 𝑘1 𝑘2 𝑘1 𝑘2
1 1 1 1 1 1
+ + +
[𝑘1 𝑘1 𝑘2 𝑘1 𝑘2 𝑘3 ]
Ejercicio 6.41
Derive las ecuaciones de movimiento del sistema que se muestra en la figura 6.29
utilizando las ecuaciones de Lagrange.
1 1 1
Energía cinética: T=2 𝑚1 𝑥̇ 1 2 + 2 𝑚2 𝑥̇ 2 2 + 2 𝑚3 𝑥̇ 3 2
1 1 1 1
Energía potencial: V=2 𝑘1 𝑥1 2 + 2 𝑘2 (𝑥2 − 𝑥1 )2 + 2 𝑘3 (𝑥3 − 𝑥2 )2 + 2 𝑘4 𝑥3 2
𝜕𝑇 𝑑 𝜕𝑇 𝑑
= 0; ( ) = 𝑑𝑡 (𝑚1 𝑥̇ 1 ) = 𝑚1 𝑥̈ 1
𝜕𝑥1 𝑑𝑡 𝜕𝑥̇ 1
𝜕𝑇 𝑑 𝜕𝑇 𝑑
= 0; ( ) = 𝑑𝑡 (𝑚2 𝑥̇ 2 ) = 𝑚2 𝑥̈ 2
𝜕𝑥2 𝑑𝑡 𝜕𝑥̇ 2
𝜕𝑇 𝑑 𝜕𝑇 𝑑
= 0; ( ) = 𝑑𝑡 (𝑚3 𝑥̇ 3 ) = 𝑚3 𝑥̈ 3
𝜕𝑥3 𝑑𝑡 𝜕𝑥̇ 3
𝜕𝑣
=𝑘1 𝑥1 − 𝑘2 (𝑥2 − 𝑥1 )
𝜕𝑥1
𝜕𝑣
=𝑘2 (𝑥2 − 𝑥1 )-𝑘3 (𝑥3 − 𝑥2 )
𝜕𝑥2
𝜕𝑣
=𝑘3 (𝑥3 − 𝑥2 )+𝑘4 𝑥3
𝜕𝑥3
INGENIERIA CIVIL 24
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Ecuación de lagrange:
𝑑 𝜕𝑇 𝜕𝑇 𝜕𝑉
( ) − 𝜕𝑥 + 𝜕𝑥 =0
𝑑𝑡 𝜕𝑥̇ 𝑗 𝑗 𝑗
Remplacemos cada uno en la ecuación de lagrange.
𝑚1 𝑥̈ 1 + (𝑘1 + 𝑘2 )𝑥1 − 𝑘2 𝑥2 = 0
{𝑚2 𝑥̈ 2 − 𝑘2 𝑥1 + 𝑥2 (𝑘2 + 𝑘3 ) − 𝑘3 𝑥3 = 0
𝑚3 𝑥̈ 3 − 𝑘3 𝑥2 + 𝑥3 (𝑘3 + 𝑘4 ) = 0
𝜕𝑣
=𝑘1 𝑥1 − 𝑘(𝑥2 − 𝑥1 ) − 5𝑘(𝑥3 − 𝑥1 ) = 7𝑘𝑥1 − 𝑘𝑥2 − 5𝑘𝑥3
𝜕𝑥1
𝜕𝑣
=𝑘(𝑥2 − 𝑥1 ) - 𝑘(𝑥3 − 𝑥2 ) = −𝑘𝑥1 + 2𝑘𝑥2 − 𝑘𝑥3
𝜕𝑥2
𝜕𝑣
=𝑘(𝑥3 − 𝑥2 )+𝑘𝑥3 +5k(𝑥3 − 𝑥1 )=−5𝑘𝑥1 − 𝑘𝑥2 +7𝑘𝑥1
𝜕𝑥3
Ecuación de lagrange:
𝑑 𝜕𝑇 𝜕𝑇 𝜕𝑉
( ) − 𝜕𝑥 + 𝜕𝑥 =0
𝑑𝑡 𝜕𝑥̇ 𝑗 𝑗 𝑗
[𝑚]{𝑥̈ } + [𝑘]{𝑥} = 0
Ejercicio 6.55
La ecuación de frecuencia de un sistema de tres grados de libertad está dada por
𝜆 − 5 −3 −2
| −3 𝜆 − 6 −4 |=0
−1 −2 𝜆 − 6
Encuentre las raíces de esta ecuación.
Solucion:
−3 −4
(𝜆 − 5) |𝜆 − 6 −4
| + 3| |−2|
−3 𝜆 − 6
|=0
−2 𝜆−6 −1 𝜆 − 6 −1 −2
INGENIERIA CIVIL 25
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
𝜆3 − 17𝜆2 + 77𝜆 − 98 = 0
De donde tenemos:
𝜆1 =2.21398
𝜆2 =4.16929
𝜆3 = 10.6168
Ejercicio 6.56
Determine los valores eigen y vectores eigen del sistema que se muestra en la figura 6.29, con
𝑘1 = 𝑘2 = 𝑘3 = 𝑘4 = 𝑘 𝑦 𝑚1 = 𝑚2 = 𝑚3 = 𝑚
Solución:
Para el problema 6.29 tenemos.
𝑚1 0 0 𝑘1 + 𝑘2 −𝑘2 0
[m] = [ 0 𝑚2 0 ] ; [𝑘] = [ −𝑘2 𝑘2 + 𝑘3 −𝑘3 ] y [𝑚]{𝑥̈ } + [𝑘]{𝑥} = 0
0 0 𝑚3 0 −𝑘3 𝑘3 + 𝑘4
Para 𝑚1 = 𝑚2 = 𝑚3 = 𝑚 y 𝑘1 = 𝑘2 = 𝑘3 = 𝑘4 = 𝑘 𝑦
1 0 0 2 −1 0
[𝑚] = 𝑚 [0 1 0] [𝑘] = 𝑘 [−1 2 −1]
0 0 1 0 −1 2
[𝑚]{𝑥̈ } + [𝑥]{𝑥} = 0
Para movimiento armónico 𝑥𝑖 (𝑡) = 𝑥𝑖 cos(𝜔𝑡 + 𝜙)
−[𝑚]𝜔2 {𝑥} + [𝑘]{𝑥} = 0
([𝑘] − 𝜔2 [𝑚]){𝑥} = 0
[𝑚]−1 [𝑘] − 𝜔2 𝐼 = 0 Problemas de auto valores y autovectores.
Det(A-𝜆𝐼) = 0
A=[𝑚]−1 [𝑘]
𝜔2 𝑚
𝜆=
𝑘
7𝑘
𝜔=√
𝑚
1 0 0
[𝑚] = [0 1 0]
0 0 1
INGENIERIA CIVIL 26
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
2 −1 0
[𝑘] = [−1 2 −1]
0 −1 2
3.2469796 0 0 𝜆1 = 3.2469796
Autovalores. [ 0 1.5549581 0 ]⇒ 2 = 1.5549581
𝜆
0 0 0.1980623 𝜆3 = 0.1980623
1 1 1
Autovectores [−1,2469796 0.4450419 1.8019377]
0.5549581 −0.8019377 2.2469796
𝑘
𝜔1 = 1.8019377√
𝑚
𝑘
𝜔1 = 1.2469796√
𝑚
𝑘
𝜔3 = 0.4450419√
𝑚
Ejercicio 6.60
Demuestre que las frecuencias naturales del sistema que se muestra en la figura 6.6(a) con
𝑘1 = 3𝑘, 𝑘2 = 𝑘3 = 𝑘, 𝑚1 = 4𝑚, 𝑚2 = 2𝑚 𝑦 𝑚3 = 𝑚, 𝑠𝑜𝑛 𝜔1 =
0.46√𝑘⁄𝑚; 𝜔2 = √𝑘⁄𝑚 y 𝜔3 = 1.34√𝑘⁄𝑚 .
Encuentre los vectores eigen del sistema.
𝑚1 0 0 𝑘1 +𝑘2 −𝑘2 0
[𝑚] = [ 0 𝑚2 0 ] 𝑦 [𝑘] [ −𝑘2 𝑘2 +𝑘3 −𝑘3 ]
0 0 𝑚3 0 −𝑘3 𝑘3
. 𝑚1 = 4𝑚 . 𝑘1 = 3𝑘 4 0 0 4 −1 0
. 𝑚2 = 2𝑚 . 𝑘2 = 𝑘 ⇒[𝑚] = 𝑚 [0 2 0] [𝑘] [−1 2 −1]
. 𝑚3 = 𝑚 . 𝑘3 = 𝑘 0 0 1 0 −1 1
La ecuación es: ([𝑘] − 𝜔2 [𝑚]){𝑥} = 0, problemas de autovalores y autovectores:
det(𝐴 − 𝜆𝐼)=0
A=[𝑚]−1 [𝑘]
𝜔2𝑚 𝜆𝑚
𝜆= 𝜔=√
𝑘 𝑘
INGENIERIA CIVIL 27
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Aplicamos programa 1 (proyecto de diseño)
Resultado:
1 0 0 𝜆1 = 1
Autovalores =[0 1.790569 0 𝜆
]⇒ 2 = 1.790569
0 0 0.2094306 𝜆3 = 0.2094306
𝜆1 𝑚1 𝑚
𝜔1 = √ = √
𝑘1 𝑘
𝜆2 𝑚 𝑚
𝜔2 = √ = 1.3381√
𝑘 𝑘
𝜆3 𝑚 𝑚
𝜔3 = √ = 0.45763√
𝑘 𝑘
Ejercicio 6.63
La matriz de masa [𝑚] y la matriz de rigidez [k] de una barra uniforme son:
1 0 0 1 −1 0
𝜌𝐴𝑙 2𝐴𝐸
[𝑚] = [ 0 2 0] y [𝑘] = [ −1 2 −1]
4 𝑙
0 0 1 0 −1 1
Donde 𝜌 𝑒𝑠 𝑙𝑎 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑. A es el área de sección transversal, E es el módulo de Young
y l es la longitud de la barra. Encuentre las frecuencias naturales del sistema al
buscar las raíces de la ecuación característica. Incluso encuentre los modos
principales.
Solución:
1 0 0 1 0 0
𝜌𝐴𝑙 𝜌𝐴𝑙
[𝑚] = [ 0 2 0] = 𝐴1 [ 0 2 0]; Donde A=
4 4
0 0 1 0 0 1
2𝐴𝐸 1 −1 0 1 −1 0 2𝐴𝐸
[𝑘] = [−1 2 −1] = 𝐵 [−1 2 −1] ; 𝐷𝑜𝑛𝑑𝑒 𝐵 =
𝑙 𝑙
0 −1 1 0 −1 1
𝜌𝐴𝑙
𝐴 4 𝐴 𝜌𝑙2
= 2𝐴𝐸 ⇒𝐵=
𝐵 8𝐸
𝑙
La ecuación de movimiento:
([𝑘] − 𝑤 2 [𝑚]){𝑥} = 0; Problema de autovalores y autovectores.
𝐴
det(𝑚−1 [𝑘] − 𝜔2 𝐼) = 0 𝐴 = [𝑚]−1 [𝑘] 𝜆 = 𝜔2 𝐵
det(𝐴 − 𝜆𝐼) = 0
1 0 0 1 −1 0
[𝑚] = [0 2 0] y [𝑘] = [−1 2 −1]
0 0 1 0 −1 1
INGENIERIA CIVIL 28
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Resultado:
𝐵 8𝐸
𝜆1 = 2 ⇾ 𝜔1 = √𝜆1 √𝐴 = 1.41√𝜌𝑙2
2 0 0
Autovalores = [0 1 0] ⇒ 𝐵
0 0 0 𝜆2 = 1 ⇾ 𝜔2 = √𝜆2 √𝐴
𝜆3 = 0 ⇾ 𝜔 3 = 0
1 1 1
Autovectores= [−1 0 1]
1 −1 1
Modos.
1 1 1
(1) (2) (3)
𝑥⃗ (1) = {−1} 𝑥1 ; 𝑥⃗ (2) = { 0 } 𝑥1 𝑦 𝑥⃗ (3) = {1} 𝑥1
1 −1 1
Autovectores.
1 0 0 1 1 0
[𝑚] = [0 2 0] 𝑥1 = {−1} 𝑎1 𝑥2 = {1} 𝑎2 𝑥3 = {1} 𝑎3
0 0 1 1 1 2
Para ortonormalizar se debe: 𝑥𝑖𝑇 [𝑚]𝑥1 = 1; i=1,2,3
1 00 1
𝑥1𝑇 [𝑚]𝑥1 =[1 ]
−1 1 [0 20] [−1] 𝑎12 = 1 ⇒4𝑎12 = 1
0 01 1
1
𝑎1 = 2
1 0 0 1
𝑥2𝑇 [𝑚]𝑥2 =[1 1 1] [0 2 0] [1] 𝑎22 = 1 ⇒4𝑎22 = 1
0 0 1 1
1
𝑎2 = 2
1 0 0 0
3 =[0 1 2] [0 2 0] [1] 𝑎3 = 1 ⇒6𝑎3 = 1
𝑇 [𝑚]𝑥 2 2
𝑥3
0 0 1 2
1
𝑎3 =
√6
La matriz de ortonormal a [𝑚] es:
1 1
0
2 2
2 1 1
[𝑥] = −
2 2 √6
1 1 2
[ 2 2 √6]
INGENIERIA CIVIL 29
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
Ejercicio 6.83
Encuentre la respuesta de vibración libre de un sistema de tres grados de
libertad regido por la ecuación.
Suponga las condiciones iniciales como 𝑥𝑖 = 0.1 𝑦 𝑥𝑖 = 𝑜; 𝑖 = 1,2,3
1 0 0 2 −1 0 𝑥1 = 0,1 𝑥̇ 1 = 0
10[0 1 0] 𝑥⃗̈ (𝑡) =100[−1 2 −1] 𝑥⃗(𝑡) = 0 𝑥2 = 0,1 𝑥̇ 2 = 0
0 0 1 0 −1 1 𝑥3 = 0,1 𝑥̇ 3 = 0
10 0 0 200 −100 0
[ 0 10 0 ] 𝑥⃗̈ (𝑡)+[−100 200 −1 00] 𝑥⃗(𝑡) = 0
0 0 10 0 −100 100
Resolvemos mediante análisis modal.
i. Resolvemos el problema de eigen.
10 0 0 200 −100 0
𝜔2 [𝑚]𝑥⃗ = 0[𝑘]𝑥⃗ Donde [𝑚]=[ 0 10 0 ] [𝑘] = [−100 200 −1 00]
0 0 10 0 −100 100
Aplicamos programa 1 (proyecto de diseño)
32.469796 0 0
Resultado: autovalores=[ 0 15.5495813 0 ]
0 0 1.9806226
𝜔1 = 5.698 𝑟𝑎𝑑⁄𝑠
𝜆1 = 32.469796 𝜆2 = 15.5495 𝜆3 = 1.9806226 ⇾ 𝜔2 = 3.9432 𝑟𝑎𝑑⁄𝑠
𝜔3 = 1.4073 𝑟𝑎𝑑⁄𝑠
Matriz moda. [𝑥⃗]
0.1868935 0.2330523 0.1037181
Autovectores. [𝑥] = [ −𝑜. 23305 0.1037181 0.1868935]
0.1037181 −01868935 0.2330523
0.0575592 0
𝑞⃗𝑜 = [𝑥]𝑇 [𝑚]𝑥⃗𝑜 =[0.1498769] 𝑞⃗̇𝑜 = [𝑥]𝑇 [𝑚]𝑥⃗̇𝑜 =[0]
0.523669 0
La solución para 𝑞⃗𝑖 (𝑡) es:
0.0575592 cos(5,698219𝑡)
𝑞⃗𝑖 (𝑡) = [ 0.1498769 cos(3.943295𝑡) ]
0.5236639 cos 1.4073459𝑡)
INGENIERIA CIVIL 30
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
⃗⃗(𝒕) = [𝒙]𝒒
𝒙 ⃗⃗(𝒕)
↓
𝑥1 (𝑡) = 0.01075 cos(5.698𝑡) + 0.0349 cos(3.943𝑡) + 0.0543134 cos(1.4 073𝑡)
𝑥2 (𝑡) =-0.0134cos(5.698𝑡)+0.0155449cos(3.943𝑡)+0.09786cos(1.4073𝑡)
𝑥3 (𝑡) =0.0059cos(5.698𝑡) – 0.028011cos(3.943𝑡)+0.1220411cos(1.4073𝑡)
Ejercicio 6.80
Resolución.
𝑚1 = 𝑚 𝑘1 = 𝑘 𝑥1 (0) = 0 𝑥̇ 1 = 𝑥̇ 10
𝑚2 = 2𝑚 𝑘2 = 2𝑘 𝑥2 (0) = 0 𝑥̇ 2 = 0
𝑚3 = 3𝑚 𝑘3 = 3𝑘 𝑥3 (0) = 0 𝑥̇ 3 = 0
1 0 0 3 −2 0
Matriz masa: [𝑚] = 𝑚 [0 2 0] [𝑘] = 𝑘 [−2 5 −3]
0 0 3 0 −3 3
𝑘 𝑘
𝜆1 = 4.3860 ⇾ 𝜔1 = √𝜆1 √ = 2.09427√
𝑚 𝑚
𝑘 𝑘
𝜆2 = 0.11399 ⇾ 𝜔2 = √𝜆2 √ = 0.33763√
𝑚 𝑚
𝑘 𝑘
𝜆3 = 2 ⇾ 𝜔3 = √𝜆3 √ = 1.414213√
𝑚 𝑚
INGENIERIA CIVIL 31
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
1
Primer modo: 𝑥 (1)
={−0.443005}
1.628667
1
Segundo modo: 𝑥 (2)
={−0.693005}
0.2046664
1
Tercer modo: 𝑥 (3)
={ 0.5 }
−0.5
[𝑥]𝑇 [𝑚][𝑥] = 𝐼
0.6923 0. .2760 0.6667
Matriz modal. [𝑥] = [−0.4798 0.3983 0.3333 ]
0.1417 0.4496 −0.3333
La solución está dada por: 𝑥⃗(𝑡) = [𝑥]𝑞⃗(𝑡), donde los 𝑞⃗(𝑡) son soluciones de:
INGENIERIA CIVIL 32
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
𝑥⃗(𝑡) = [𝑥]𝑞⃗(𝑡)
𝑚 𝑚 𝑚
𝑥1 (𝑡)=0.2288𝑥̇10 √ sin(2.09427𝑡)+0.3142𝑥̇10 √ sin(1.4142𝑡)+0.2257𝑥̇10 √ sin(0.333763𝑡)
𝑘 𝑘 𝑘
𝑚 𝑚 𝑚
𝑥2 (𝑡)=0.1586𝑥̇10 √ sin(2.09427𝑡)+0.15713𝑥̇10 √ sin(1.4142𝑡)+0.3257𝑥̇10 √ sin(0.333763𝑡)
𝑘 𝑘 𝑘
𝑚 𝑚 𝑚
𝑥2 (𝑡)=0.046895𝑥̇10 √ sin(2.09427𝑡)–0.15713𝑥̇10 √ sin(1.4142𝑡)+6.3676023𝑥̇10 √ sin(0.333763𝑡)
𝑘 𝑘 𝑘
Ejercicio 6.12
Ejercicio 6.12
Ejercicio 6.12
Ejercicio 6.12
Ejercicio 6.12
Ejercicio 6.12
Ejercicio 6.12
Ejercicio 6.12
INGENIERIA CIVIL 33
CAPITULO 1: FUNDAMENTOS DE VIBRACIÓN
INGENIERIA CIVIL 34