INFORME
INFORME
INFORME
FIC-2017
PORTIFIC
INTEGRANTES:
General
Introducción………................................................................................................04
Objetivos…………………………………………………………………………………05
Informe Técnico…………………………………………………………………………06
ANEXOS
Introducción
m1 0 0
M := 0 m2 0
0 0 m3
0.767 0 0
M= 0 0.757 0
0 0 0.738
Matriz de Rigides:
Ingresando los datos de las rigidecez equivalentes de cada piso calculadas anteriormente,
formaremos la matriz de rigides de la estructura para un modelo de MASA - CORTANTE como
se muestra a continuacion:
k1 + k2 -k2 0
K := -k2 k2 + k3 -k3
0 -k3 k3
5688.69 -2844.345 0
K = -2844.345 5688.69 -2844.345
0 -2844.345 2844.345
Utilizando el calculo por valores y vectores porpios de una matriz se procedera a calcular las
frecuencias, Periodos y las formas modales, para ello a partir de la ecuacion de movimiento
obtendremos la matriz K1 definida a continuacion:
-1
K1 := M K
7416.806 -3708.403 0
K1 = -3757.391 7514.782 -3757.391
0 -3854.126 3854.126
A la matriz K1 le calcularemos sus valores propios los cuales so equivalentes a una matriz de
las frecuencias elevadas al cuadrado, esto lo realizaremos en el Mathcad utilizando el comando
eifenvals pero el resultado se genera desordenadamente, a esta matriz la denominaremos evp.
Procedemos a ordenar los datos de menor a mayor utilizando el comando sort del Mathcad y a
esta matriz la denominaremos Freq.
Para el calculo de las frecuencias naturales de la estructura sacaremos la raiz cuadrada a los
valores de la matriz Freq y los denominaremos ωi, para ello extraeremos los datos de la matriz
definiendo primero un valor de rango i (0,1,2) para luego efectuar la operacion.
Con los valores de las frecuencias procederemos al calculo de los periodos de vibracion
denominados Ti.
i := 0 .. 2
2π
ω := Freq T :=
i i i ω
i
ω = 27.448 T = 0.229
0 0
ω = 76.503 T = 0.082
1 1
ω = 110.362 T = 0.057
2 2
El siguiente paso sera el calculo de las formas modales de vibracion, para ello calcularemos los
vectores propios de la matriz K1 utilizando el comando eigenvecs del Mathcad y denominaremos
a esta matriz como v.
v := eigenvecs( K1)
1 1 1
V = 1.797 V = 0.422 V = -1.284
0 1 2
2.233 -0.813 0.595
Una vez ordenados y normalizados las formas modales procedemos a ensamblar la matriz de
modos a la que denominaremos ϕ como se muestra a continuacion:
1 1 1
ϕ = 1.797 0.422 -1.284
2.233 -0.813 0.595
6.892 0 0
T
Me := ϕ M ϕ Me = 0 1.39 0
0 0 2.277
T 5192.519 -0 -0
Ke := ϕ K ϕ
Ke = -0 8134.991 0
-0 0 27728.276
1 3.775
I := 1 L = 0.486
T
L := ϕ M I
1 0.234
2
L( 0 , 0)
Le := Le = 2.068
0, 0 Me 0, 0
( 0 , 0)
2
L( 1 , 0)
Le := Le = 0.17
1, 0 Me 1, 0
( 1 , 1)
2
L( 2 , 0)
Le := Le = 0.024
2, 0 Me 2, 0
( 2 , 2)
2
Le
Mtotal :=
i , 0
Mtotal = 2.262
i = 0
L
( 0 , 0)
r := r = 0.548
0 Me 0
( 0 , 0)
L
( 1 , 0)
r := r = 0.35
1 Me 1
( 1 , 1)
L
( 2 , 0)
r := r = 0.103
2 Me 2
( 2 , 2)
r =1
i
i=0
P :=
(Le0 , 0)100 P = 91.429
0 Mtotal 0
P :=
(Le1 , 0)100 P = 7.512
1 Mtotal 1
( )
P :=
(Le2 , 0)100 P = 1.059
2 Mtotal 2
2
P
Suma :=
i Suma = 100
i = 0