Tarea N°1: “Cálculo de propiedades dinámica de edificios de corte”
Integrantes: Magdalena Espinoza, Carla Ferrada, Silvia Gálvez.
Introducción al problema
Uno de los modelos más comunes de análisis dinámico es suponer que la masa de la estructura está
concentrada en varios puntos (una masa concentrada por cada nivel) y que se tendrán
desplazamientos laterales únicamente en cada nivel. Este comportamiento está gobernado por la
ecuación del movimiento. Cuando el sistema estudiado no posee amortiguamiento y consiste en
una vibración libre (sin acción de fuerzas externas) se tendrá la siguiente ecuación:
[𝑀]𝑦̈ + [𝐾]𝑦 = 0
En donde:
[M]= Es la matriz de masas concentradas 𝑦̈ = es el vector de aceleración de la estructura
[K]= Es la matriz de rigidez 𝑦 = es el vector de desplazamiento
Este último caso es el que representa el problema propuesto en donde se entrega la elevación de
un edificio comercial de marcos de hormigón armado con diafragma rígido (imagen 1), para el cual
se debe realizar un informe técnico con el cálculo de las frecuencias, modos de vibrar, matriz modal,
matriz modal normal y la verificación de ortogonalidad de las matrices de masas y rigidez.
Imagen 1. Modelo N° 4.
Para plantear la ecuación del movimiento se requiere que la estructura cumpla con los siguientes
supuestos: 1. La estructura se compone de diafragmas rígidos.
2. Las columnas son inextensibles.
3. La masa traslacional se concentra a nivel de cada piso.
Con esto, es posible determinar la matriz de masas concentradas y la matriz de rigidez con ayuda
del software Mathcad.
1
La ecuación del módulo de elasticidad se obtiene del ACI, y las inercias se obtienen considerando
una sección transversal cuadrada para los dos tipos de columnas presentes.
Con estos datos es posible obtener la rigidez horizontal de cada columna, como todas se comportan
como marco rígido y hay dos opciones de dimensiones solo se tendrán dos valores posibles, a partir
de los cuales se podrá construir la matriz de rigidez del sistema.
Cálculo de las frecuencias y modos de vibrar
Para el cálculo de las frecuencias y modos de vibrar se debe desarrollar un problema de valores y
vectores propios. Se conocen las matrices de rigidez y masas, las cuales son reales, simétricas y
definidas positivas, por lo que siempre se pueden obtener valores de frecuencias reales.
Realizando los cálculos correspondientes se obtienen las frecuencias (𝜔𝑖 ), ordenando las
frecuencias de menor a mayor y agrupándolas en una matriz se conforma la matriz espectral (Ω).
16.4414 16.4414 0 0 0
47.5626 𝑟𝑎𝑑
Ω=[ 0 47.5626 0 0 𝑟𝑎𝑑
𝜔=[ ][ ] , ][ ]
70.5503 𝑠 0 0 70.5503 0 𝑠
89.7139 0 0 0 89.7139
Luego para calcular las componentes del vector de amplitudes (Φ𝑖𝑗 ) se procede a resolver la
siguiente ecuación: (𝐾 − 𝜔𝑖2 ∙ 𝑀) ∙ Φ𝑖 = 0 , 𝑖 = 1,2,3,4.
Según los resultados obtenidos con Matlab, para el sistema en estudio de 4 grados de libertad se
tiene:
1.0000 1.0000
Para i=1, 𝜔1 = 16.4414, Φ1 = [ 1.9255 ] Para i=2, 𝜔2 = 47.5626, Φ2 = [ 1.3762 ]
2.4526 0.0593
2.7558 −2.0147
2
1.0000 1.0000
Para i=3, 𝜔3 = 70.5503, Φ3 = [ 0.6275 ] Para i=4, 𝜔4 = 89.7139, Φ4 = [ −0.2194 ]
−1.6472 0.1028
0.9694 −0.0225
Agrupando los diferentes modos comenzando por el asociado a la frecuencia fundamental y
terminando con la frecuencia máxima se obtiene la matriz modal, que corresponde a:
1.0000 1.0000 1.0000 1.0000
Φ = [1.9255 1.3762 0.6275 −0.2194 ]
2.4526 0.0593 −1.6472 0.1028
2.7558 −2.0147 0.9694 −0.0225
En Matlab este cálculo se realiza utilizado el comando eig, el cual devuelve directamente los valores
y vectores propios, almacenando los valores propios en la diagonal de una matriz, correspondiente
a la matriz espectral al cuadrado y almacenando los vectores propios en la matriz modal.
Dibujo a escala de los respectivos modos de vibrar
En la imagen 2 se representan los distintos modos de vibrar asociados a las frecuencias (𝜔𝑖 ), fueron
realizados en Matlab a partir de la matriz modal obtenida.
Imagen 2. Representación de los modos de vibrar.
3
Diagrama de flujo
Inicio
M, K
Cálculo de la matriz modal (Φ) y
espectral al cuadrado (Ω 2 ).
Cálculo de frecuencias wi a partir de √Ω2
wi
Cálculo de los periodos: Ti=2π/wi
Normalización de la Matriz modal con respecto al primer piso
Matriz modal: Φ
Cálculo matriz masa generalizada: Φ𝑗 𝑇 ∗ [𝑀] ∗ Φ𝑗 = [𝑀]∗𝑗
Calcular √[𝑀]∗𝑗
Dividir elementos de la matriz modal por √[𝑀]∗𝑗
Matriz modal-normal: Φ
Verificación de ortogonalidad
2
𝐼, Ω
FIN
4
Cálculo de la matriz modal-normal
Para poder escribir de forma más compacta las ecuaciones de movimiento se recomienda
normalizar la matriz modal. Este procedimiento se realiza diviendo la matriz modal por un factor de
escala que corresponde a la raíz cuadrada de la diagonal de la matriz de masa generalizada. La masa
generalizada y los factores de escalas entregados por el sript son:
Finalmente, la matriz modal normalizada queda de la siguiente forma:
Verificación propiedades de ortogonalidad
Los vectores modales normalizados satisfacen las siguientes condiciones importantes de
ortogonalidad:
Φ𝑖 𝑇 ∗ [𝑀] ∗ Φ𝑗 = 0 𝑝𝑎𝑟𝑎 𝑖 ≠ 𝑗 Φ𝑖 𝑇 ∗ [𝐾] ∗ Φ𝑗 = 0 𝑝𝑎𝑟𝑎 𝑖 ≠ 𝑗
Φ𝑖 𝑇 ∗ [𝑀] ∗ Φ𝑗 = 1 𝑝𝑎𝑟𝑎 𝑖 = 𝑗 Φ𝑖 𝑇 ∗ [𝐾] ∗ Φ𝑗 = 𝜔𝑖2 𝑝𝑎𝑟𝑎 𝑖 = 𝑗
donde Φ𝑖 y Φ𝑗 corresponden a vectores modales normalizados.
Las relaciones anteriores son equivalentes a: Φ𝑇 ∗ [𝑀] ∗ Φ = [Ι]
Φ𝑇 ∗ [𝐾] ∗ Φ = [Ω2 ]
Estas son las ecuaciones que resolverá el sript. Remplazando las matrices por las correspondientes
al problema propuesto se tiene:
0.0269708 0.0519314 0.0661487 0.0743268 30 0 0 0 0.0269708 0.0455413 0.053701 0.16628 1 0 0 0
[0.0455413 0.0626748 0.0027023 − 0.091752 ] [ 0 110 0 0 ] [0.0519314 0.06626748 0.033699 −0.036473 ] = [0 1 0 0]
0.053701 0.033699 −0.088457 0.052056 0 0 80 0 0.0661487 0.0027023 −0.088457 0.017099 0 0 1 0
0.16628 −0.036473 0.017099 0.0037387 0 0 0 60 0.0743268 −0.091752 0.052056 0.0037387 0 0 0 1
0.0269708 0.0519314 0.0661487 0.0743268 217593 −108796 0 0 0.0269708 0.0455413 0.053701 0.16628
[0.0455413 0.0626748 0.0027023 − 0.091752 ] ∗ [−108796 314291 −163194 −14100 ] ∗ [0.0519314 0.06626748 0.033699 −0.036473 ] = [Ω2 ]
0.053701 0.033699 −0.088457 0.052056 0 −163194 271991 −108796 0.0661487 0.0027023 −0.088457 0.017099
0.16628 −0.036473 0.017099 0.0037387 0 −14100 −108796 122896 0.0743268 −0.091752 0.052056 0.0037387
Finalmente, con los resultados obtenidos en matlab se verifica que las matrices de masa y rigidez
cumplen con la propiedades de ortogonalidad con respecto a los vectores modales normales.
5
Conclusiones
Las edificaciones presentan distintos modos de vibración y a cada uno le corresponde un periodo
de vibración diferente, el cual depende de las características propias de la estructura, como altura,
masa y rigidez.
El período fundamental de la estructura analizada es de 0,38 [s]. Este resultado es uno de los más
importantes del cálculo realizado ya que corresponde al período de mayor duración, esto significa
que es el que generará mayores daños a la estructura (por mayores desplazamientos) y por ende el
que mejor se debe analizar ya que permite predecir las fuerzas a aplicar sobre la estructura con el
fin de dimensionar su sistema de resistencia sísmica. Además, podemos comprobar que a medida
que el período disminuye, los modos de vibrar se van haciendo más insignificantes, lo que explica
que para edificios con más grados de libertad se busque un método que permita despreciar los
modos de vibración que no contribuyan al modelo.
Anexos
Script en Matlab
clear all;close all;clc a(:,3)=[a13,a23,a33,a43];
%Del libro Basic Structural Dynamics. Anderson and a(:,4)=[a14,a24,a34,a44];
Naeim,2012 disp('MATRIZ MODAL');disp(a)
%Modos de vibrar en un sistema de 4 grados de %Calcula el factor de normalizacion asociado a la
libertad masa
%Ingreso Matriz de Masa Concentrada en Toneladas aMa=a'*M*a; %matriz de masa generalizada
M=[30 0 0 0; 0 110 0 0; 0 0 80 0; 0 0 0 60]; norm_1=sqrt(aMa(1,1));
%Ingreso Matriz de Rigidez Horizontal en kN/m norm_2=sqrt(aMa(2,2));
K=[217593 -108796 0 0; -108796 314291 -163194 - norm_3=sqrt(aMa(3,3));
14100; norm_4=sqrt(aMa(4,4));
0 -163194 271991 -108796; 0 -14100 -108796 122896]; %Divide para normalizar
%Calculo de la Matriz modal (a)y Matriz espectral al nor_phi(:,1)=1./norm_1.*a(:,1);
cuadrado (D) nor_phi(:,2)=1./norm_2.*a(:,2);
[a,D]=eig(K,M) nor_phi(:,3)=1./norm_3.*a(:,3);
%Calculo de las frecuencias wi nor_phi(:,4)=1./norm_4.*a(:,4);
[omegas, k]=sort(sqrt(diag(D))); %Matriz modal-normal
disp('FRECUENCIAS');disp(omegas) nor_phi;
T=2*i./omegas; disp('MATRIZ MODAL NORMAL');disp(nor_phi)
%Calculo de los periodos modo 1, 2, 3 y 4 %Verifica que la matriz de masa sea unitaria y la
T1=2*pi./omegas(1,1); %matriz rigidez igual a la espectral al cuadrado.
T2=2*pi./omegas(2,1); orth_M=nor_phi'*M*nor_phi;
T3=2*pi./omegas(3,1); ortt_k=nor_phi'*K*nor_phi;
T4=2*pi./omegas(4,1); disp('VERIFICACION ORTOGONALIDAD')
%Normalizacion de la Matriz Modal disp('igual a matriz identidad');disp(orth_M)
%ai elementos de las columnas disp('Matriz espectral al cuadrado');disp(round(D))
a1=a(:,1); disp('verificar igual a matriz espectral al
a2=a(:,2); cuadrado')
a3=a(:,3); disp(round(ortt_k))
a4=a(:,4);
%Normaliza 1°columna matriz modal (c/r al piso 1) %Gráfico de modos de vibrar
a11=1; figure (1)
a21=a1(2,1)./a1(1,1); circ(:,1)=[3 6 9 12];
a31=a1(3,1)./a1(1,1); subplot(2,2,1),plot([0;a(:,1)],[0;circ(:,1)],'b',[0
a41=a1(4,1)./a1(1,1); 0 0 0 0],[0 3 6 9 12],'ko-
%Normaliza 2°columna matriz modal (c/r al piso 1) ','MarkerFaceColor','k'),title('Modo 1')
a12=1; xlim ([-4 4])
a22=a2(2,1)./a2(1,1); ylim ([0 12])
a32=a2(3,1)./a2(1,1); subplot(2,2,2),plot([0;a(:,2)],[0;circ(:,1)],'r',[0
a42=a2(4,1)./a2(1,1); 0 0 0 0],[0 3 6 9 12],'ko-
%Normaliza 3°columna matriz modal (c/r al piso 1) ','MarkerFaceColor','k');title('Modo 2')
a13=1; xlim ([-4 4])
a23=a3(2,1)./a3(1,1); ylim ([0 12])
a33=a3(3,1)./a3(1,1); subplot(2,2,3),plot([0;a(:,3)],[0;circ(:,1)],'g',[0
a43=a3(4,1)./a3(1,1); 0 0 0 0],[0 3 6 9 12],'ko-
%Normaliza 4°columna matriz modal (c/r al piso 1) ','MarkerFaceColor','k');title('Modo 3')
a14=1; xlim ([-4 4])
a24=a4(2,1)./a4(1,1); ylim ([0 12])
a34=a4(3,1)./a4(1,1); subplot(2,2,4),plot([0;a(:,4)],[0;circ(:,1)],'c',[0
a44=a4(4,1)./a4(1,1); 0 0 0 0],[0 3 6 9 12],'ko-
%Matriz Modal ','MarkerFaceColor','k');title('Modo 4')
a=[]; xlim ([-4 4])
a(:,1)=[a11,a21,a31,a41]; ylim ([0 12])
a(:,2)=[a12,a22,a32,a42];