UNIVERSIDAD DE CONCEPCIÓN
FACULTAD DE INGENIERÍA, DEPARTAMENTO DE INGENIERÍA CIVIL
MAGISTER EN CIENCIAS DE LA INGENIERÍA, MENCIÓN INGENIERÍA
CIVIL
TAREA N°1
Cálculo de propiedades dinámicas edificio de corte
Nicool Astroza Méndez - Nicole Ávila Vidal
Dinámica de Estructura
Concepción, 27 de Abril 2020
Universidad de Concepción
Facultad de Ingeniería, Departamento de Ingeniería Civil
Programa Magíster en Ciencias de la Ingeniería, mención Ingeniería Civil
Desarrollo Modelo 2
La Figura 1 corresponde a un edificio comercial de cuatro pisos en donde cada piso se puede representar por
masas concentradas (figura 2). El edificio se encuentra estructurado en base a marcos de hormigón armado
además de tener una altura entrepiso de 3 m, y una separación entre columnas de 6 m. La calidad del
hormigón corresponde a un H30 con una resistencia característica de 𝑓′𝑐 = 25 𝑀𝑃𝑎, donde su módulo de
elasticidad se determinó según lo dispuesto en el ACI318-14.
Las columnas son inextensibles, esto se refiere a elementos indeformables desde el punto de vista axial.
También a nivel de cada piso existe un diafragma horizontal rígido, no existe deformación, solamente se
traslada al ser sometido a un desplazamiento relativo horizontal.
En el caso de la estructura (figura 1), es necesario construir la matriz de rigidez horizontal del sistema y la
matriz de masa, para obtener ambas es fundamental iniciar con un desplazamiento unitario en el grado de
libertad correspondiente y dejando fijo los grados de libertad restantes de la estructura, como se ve en la
Figura 3.
25 T 25 T 25 T 2525TT 25 T
m4=75 T m4=75 T
3m
3m
50 X 50 50 X 50
50 X 50 50 X50
50X 50 50 X 50 50 X 50 50 X 50
50 T 50 T50 T 5050
TT 50 T m3=150 T m3=150 T
3m
50 X 50
3m
50 X 50 50 X 50
50 X 50
50 T 50 X 50
50 T m2= 50 T
50 X 50 m2= 50 T
50 X 50 50 X 50
3m
50 X 50 50 X 50
50 X 50
3m
50 X 50
50 T 50 T 50 T m1=150 T
50 T 50 T 50 T m1=150 T
50 X 50
3m
50 X 50 50 X 50 50 X 50
50 X 50
3m
50 X 50 50 X 50 50 X 50
EMP EMP EMP EMP
EMP EMP EMP EMP
Figura 1: Estructura edificio comercial, modelo 2 Figura 2: Diagrama de cuerpo libre
Figura 3: (a) Fijo piso 1; (b) Fijo piso 2; (c) Fijo piso 3; (d) Fijo piso 4
La matriz de masas concentradas y la matriz de rigidez horizontal del sistema obtenida de la Figura 2 y Figura
3 son las siguientes:
150 0 0 0
𝑀=[ 0 50 0 0 ] ( 𝑇)
0 0 150 0
0 0 0 75
6𝑘 + 2𝑘1 −2𝑘 −2𝑘1 0 339990 −108800 0 0
𝐾 = [ −2𝑘 4𝑘 −2𝑘 0 ] (𝑘𝑁 ) = [−108800 217590 −108800 0 𝑘𝑁
]( )
−2𝑘1 −2𝑘 6𝑘 + 2𝑘1 −4𝑘 𝑚 0 −108800 339990 −217590 𝑚
0 0 −4𝑘 4𝑘 0 0 −217590 2177590
Universidad de Concepción
Facultad de Ingeniería, Departamento de Ingeniería Civil
Programa Magíster en Ciencias de la Ingeniería, mención Ingeniería Civil
1. Diagrama de flujo para la programación en MATLAB de los scripts de cálculo de la matriz
modal-normal.
Iniciar análisis
Determinar [K], [M]
Determinar valores y vectores
propios
Reordenamiento de las
frecuencias (𝜔𝑖 )
Normalización de matriz
modal, con respecto al piso 1
Determinar matriz de masa
generalizada [M*]
Cálculo de la matriz modal
normal
Verificación propiedades de
ortogonalidad de [M] y [K]
2. Frecuencias y modos de vibrar.
- Frecuencia modal:
Para la determinación de la frecuencia de los diferentes modos de vibrar de la estructura se resolvió el
problema característico de la estructura mediante la siguiente ecuación:
det([𝐾 ] − 𝜔2 [𝑀]) = 0
Como resultado se obtuvieron los siguientes autovalores, los que corresponden a la frecuencia de cada modo
de vibrar.
𝜔1 14.08
𝜔2
Ω = (𝜔 ) = (42.20) (𝑟𝑎𝑑⁄𝑠)
3 65.58
𝜔4 75.08
Universidad de Concepción
Facultad de Ingeniería, Departamento de Ingeniería Civil
Programa Magíster en Ciencias de la Ingeniería, mención Ingeniería Civil
- Modos de Vibrar:
En base a los autovalores obtenidos en el paso anterior, se obtuvieron los autovectores asociados a cada modo
resolviendo el siguiente sistema:
([𝐾 ] − 𝜔2 [𝑀]){∅} = 0
Como resultado se obtuvieron los siguientes autovectores de la estructura. Estos autovectores se encuentran
normalizados con respecto al piso 1.
𝜙11 𝜙12 𝜙13 𝜙14 1 1 1 1
𝜙 𝜙22 𝜙23 𝜙24 2.40 0.69 −2.47 −4.88
𝜙 = ( 21 )=( )
𝜙31 𝜙32 𝜙33 𝜙34 3.59 −0.18 −1.21 1.88
𝜙41 𝜙42 𝜙43 𝜙44 3.85 −0.47 2.76 −1.20
Donde la primera columna de la matriz está asociada al modo de vibrar 1, por ende a la frecuencia 1, la
segunda columna está asociada al modo de vibrar 2, la tercera al modo de vibrar 3 y la cuarta al modo de
vibrar 4.
3. Dibujo a escala de los respectivos modos de vibrar.
En la siguiente figura se presentan de forma esquemática los diferentes modos de vibrar que presenta la
estructura en estudio.
Figura 4
4. Matriz modal-normal.
Para la determinación de la matriz modal normal es necesario en primera instancia determinar la matriz de
masa generalizada:
3480.9 0 0 0
0 195.5 0 0
𝜙𝑇 ∗ 𝑀 ∗ 𝜙 = ( )
0 0 1244.7 0
0 0 0 2173.1
Luego se resuelve lo siguiente y se obtiene la matriz modal normal:
1 1 1 1
√3480.9 √195.5 √1244.7 √2173.1
2.40 0.69 −2.47 −4.88 0.017 0.072 0.028 0.022
√3480.9 √195.5 √1244.7 √2173.1 0.041 0.050 −0.070 −0.105
𝜙= 3.59 −0.18 −1.21 1.88 =( )
0.061 −0.013 −0.034 0.040
√3480.9 √195.5 √1244.7 √2173.1 0.065 −0.034 0.078 −0.043
3.85 −0.47 2.76 −1.20
(√3480.9 √195.5 √1244.7 √2173.1)
Universidad de Concepción
Facultad de Ingeniería, Departamento de Ingeniería Civil
Programa Magíster en Ciencias de la Ingeniería, mención Ingeniería Civil
5. Verificación de las propiedades de ortogonalidad de la matriz de masa y rigidez.
- Verificación propiedades de ortogonalidad de la matriz de masa:
1 0 0 0
0 1 0 0
𝜙𝑇 ∗ 𝑀 ∗ 𝜙 = ( )
0 0 1 0
0 0 0 1
- Verificación propiedades de ortogonalidad de la matriz de rigidez:
𝜙𝑇 ∗ 𝐾 ∗ 𝜙 = Ω2
198.3 0 0 0 198.3 0 0 0
0 1780.6 0 0 0 1780.6 0 0
( )=( )
0 0 4170.4 0 0 0 4170.4 0
0 0 0 5637.0 0 0 0 5637.0
Dado lo anterior, se determina que se cumplen con las propiedades de ortogonalidad para la matriz de masa y
la matriz de rigidez. Siendo ɸ la matriz de masa generalizada.
6. Conclusiones
Se analiza la estructura en estudio con la finalidad de obtener sus propiedades dinámicas determinando
rigidez, frecuencia y los modos de vibrar asociados a las frecuencias.
En base al análisis de la matriz de rigidez se observa que el piso 1 y 3 presentan la mayor rigidez, esto está
relacionado a que dichos pisos están conectados a una mayor cantidad de columnas oponiéndose a su
movimiento. Se observa que dicha matriz es positiva en su diagonal y presenta una forma bandeada y
simétrica. En el caso de la matriz de masa presenta una forma simétrica y es definida positiva.
Por otro lado, se determinó la matriz espectral y la matriz modal del sistema. En base a esta última se generó
una representación gráfica de los diferentes modos de vibrar que presenta la estructura. El modo 1 está
asociado a la frecuencia más baja, mientras que el modo 4 está asociado a la frecuencia más alta. Se observa
que a medida que aumenta la frecuencia el comportamiento de la estructura en base a los desplazamientos se
vuelve más errático. También se observa que existe una correlación entre el número del modo de vibrar con la
cantidad de nudos que éste presenta. En específico, se tiene la cantidad de nudos que presenta la estructura en
un determinado modo es igual al número del modo menos 1. Por ejemplo para el modo 2 se presenta 1 nudo,
mientras que para el modo 3 se observan 2 nudos.
Finalmente, se obtuvo la matriz modal normal de la estructura, en base a la cual se verificaron las propiedades
de ortogonalidad tanto para la matriz de masa, como para la matriz de rigidez del sistema, verificando que
ambas matrices cumplen con las propiedades de ortogonalidad.
Universidad de Concepción
Facultad de Ingeniería, Departamento de Ingeniería Civil
Programa Magíster en Ciencias de la Ingeniería, mención Ingeniería Civil
7. ANEXO 1: Scripts cálculo de la matriz modal-normal
clear,clc;
n=4;%Grados de libertad del sistema
fc=25; %Resistencia hormigón MPa
H1=3; %Altura de entrepiso en m
H2=6; %Altura columna larga en m
a=0.50; %Largo y ancho columna en m
I=(a^4)/12; %Inercia cuadrado en m4
E=4700*sqrt(fc)*1000; %Elasticidad del hormigón en kN/m
k1=(12*E*I)/(H1^3); %Determinación rigidez EMP-EMP columna corta
k2=(12*E*I)/(H2^3); %Determinación rigidez EMP-EMP columna larga
K=[6*k1+2*k2 -2*k1 -2*k2 0;-2*k1 4*k1 -2*k1 0;-2*k2 -2*k1 6*k1+2*k2 -
4*k1;0 0 -4*k1 4*k1]; %Matriz de rigidez
M=[150 0 0 0;0 50 0 0;0 0 150 0;0 0 0 75]; %Matriz de masa
%Determinación de valores y vectores propios
[MOD1,ESP2]=eig(K,M);
%Determinación de las frecuencias
[ESP,h]=sort(sqrt(diag(ESP2)));
disp('La matriz espectral es : '),disp(ESP)
%Determinación de los periodos
T=zeros(n,1);
for i=1:n
T(i)=2*pi./ESP(i);
end
%Reordenamiento de la matriz modal según la matriz de frecuencias
MOD=zeros(n,n);
for i=1:n
for j=1:n
MOD(j,i)=MOD1(j,h(i));
end
end
%Normalizacion de la Matriz Modal respecto al piso 1
for i=1:n
for j=2:n
MOD(j,i)=MOD(j,i)./MOD(1,i);
end
MOD(1,i)=1;
end
disp('La matriz modal es : '),disp(MOD)
%Determinación de la matriz de masa generalizada
MG=(MOD')*M*MOD;
%Cálculo de la matriz modal normal
MN=zeros(n,n);
for i=1:n
for j=1:n
MN(j,i)=MOD(j,i)./sqrt(MG(i,i));
end
end
disp('La matriz modal normal es : '),disp(MN)
%Propiedades de ortogonalidad respecto a la matriz modal normal
ortM1=(MN')*M*MN
ortK1=(MN')*K*MN