PREGUNTA 1 :
Una viga infinitamente rígida de longitud total L=12m está apoyada en cuatro
columnas de concreto armado. Las columnas tienen las siguientes características:
b12 ≔ 30 cm b34 ≔ 30 cm
h12 ≔ 60 cm h34 ≔ 45 cm
b12 ⋅ h12 3 b34 ⋅ h34 3
I12 ≔ ―――= ⎛⎝5.4 ⋅ 10 5 ⎞⎠ cm 4 I34 ≔ ―――= ⎛⎝2.278 ⋅ 10 5 ⎞⎠ cm 4
12 12
H1 ≔ 3.5 m
H2 ≔ 3 m
H3 ≔ 3 m
H4 ≔ 3 m
L ≔ 6.5 m + 5.5 m = 12 m
kg
f'c ≔ 210 ――
cm 2
‾‾‾‾‾‾‾
kg kg
E ≔ 15000 ⋅ f'c ⋅ ―― 2
= ⎛⎝2.174 ⋅ 10 5 ⎞⎠ ――
cm cm 2
kg kg
W1 ≔ 5000 ― W4 ≔ 4000 ―
m m
kg kg
W1 ≔ 5000 ― W4 ≔ 4000 ―
m m
kg kg
W2 ≔ 5000 ― W3 ≔ 4000 ―
m m
Masa
W1 ⋅ L kg ⋅ s 2 W4 ⋅ L kg ⋅ s 2
m1 ≔ ――― = 61.183 ――― m4 ≔ ――― = 48.946 ―――
g cm g cm
W2 ⋅ L kg ⋅ s 2 W3 ⋅ L kg ⋅ s 2
m2 ≔ ――― = 61.183 ――― m3 ≔ ――― = 48.946 ―――
g cm g cm
Matriz de masas
⎡ m1 0 0 0 ⎤ ⎡ 6.118 ⋅ 10 3 0 0 0 ⎤
⎢ ⎥ ⎢ ⎥
0 m2 0 0
⎥=⎢0 6.118 ⋅ 10 3 0 0 kg ⋅ s 2
m≔⎢ ⎥ ―――
⎢ 0 0 m3 0 ⎥ ⎢0 0 4.895 ⋅ 10 3 0 ⎥ m
⎢⎣ 0 0 0 m4 ⎥⎦ ⎢⎣ 0 0 0 4.895 ⋅ 10 3 ⎥⎦
Rigidez
12 ⋅ E ⋅ I12 kg
k1 ≔ 3 ⋅ ―――― 3
= ⎛⎝9.856 ⋅ 10 4 ⎞⎠ ――
H1 cm
12 ⋅ E ⋅ I12 kg
k2 ≔ 3 ⋅ ―――― 3
= ⎛⎝1.565 ⋅ 10 5 ⎞⎠ ――
H2 cm
12 ⋅ E ⋅ I34 kg
k3 ≔ 3 ⋅ ―――― 3
= ⎛⎝6.603 ⋅ 10 4 ⎞⎠ ――
H3 cm
12 ⋅ E ⋅ I34 kg
k4 ≔ 3 ⋅ ―――― 3
= ⎛⎝6.603 ⋅ 10 4 ⎞⎠ ――
H4 cm
⎡ k1 + k2 -k2 0 0 ⎤ ⎡ 2.551 ⋅ 10 5 -1.565 ⋅ 10 5 0 0 ⎤
⎢ ⎥ ⎢ ⎥
-k2 k2 + k3 -k3 0 5 5 4
kg
k≔⎢ ⎥ = ⎢ -1.565 ⋅ 10 2.225 ⋅ 10 -6.603 ⋅ 10 0 ⎥ ――
⎢ 0 -k3 k3 + k4 -k4 ⎥ ⎢ 0 -6.603 ⋅ 10 4 1.321 ⋅ 10 5 -6.603 ⋅ 10 4 ⎥ cm
⎢ 0 0 -k4 k4 ⎥⎦ ⎢⎣ 0 0 -6.603 ⋅ 10 4 6.603 ⋅ 10 4 ⎥⎦
⎣
Análisis modal
Frecuencia angular
⎡ 4.169 ⋅ 10 3 -2.558 ⋅ 10 3 0 0 ⎤
⎢ 3 3 3 ⎥
-2.558 ⋅ 10 3.637 ⋅ 10 -1.079 ⋅ 10 0 1
A ≔ m -1 ⋅ k = ⎢ ⎥―
⎢ 0 -1.349 ⋅ 10 3 2.698 ⋅ 10 3 -1.349 ⋅ 10 3 ⎥ s2
⎢⎣ 0 0 -1.349 ⋅ 10 3 1.349 ⋅ 10 3 ⎥⎦
⎡ 6.663 ⋅ 10 3 ⎤
⎢ 3 ⎥ 1
B ≔ eigenvals ((A)) = ⎢ 3.621 ⋅ 10 3 ⎥―
⎢ 1.337 ⋅ 10 ⎥ s2
⎢⎣ 232.515 ⎥⎦
⎡ 6.663 ⋅ 10 3 ⎤
⎢ 3 ⎥ 1
ω2 = B ≔ eigenvals ((A)) = ⎢ 3.621 ⋅ 10 3 ⎥―
⎢ 1.337 ⋅ 10 ⎥ s2
⎢⎣ 232.515 ⎥⎦
⎡ 81.626 ⎤
⎢ 60.174 ⎥ 1
ω ≔ ‾‾
B =⎢ ⎥―
⎢ 36.561 ⎥ s
⎣ 15.248 ⎦
Periodos
rad 2⋅π
ω1 ≔ ω ((3)) = 15.248 ―― T1 ≔ ―― = 0.412 s
s ω1
rad 2⋅π
ω2 ≔ ω ((2)) = 36.561 ―― T2 ≔ ―― = 0.172 s
s ω2
rad 2⋅π
ω3 ≔ ω ((1)) = 60.174 ―― T3 ≔ ―― = 0.104 s
s ω3
rad 2⋅π
ω4 ≔ ω ((0)) = 81.626 ―― T4 ≔ ―― = 0.077 s
s ω4
Modos de la estructura
⎡ 0.232 ⎤
⎢ 0.357 ⎥
Forma del primer modo ϕ1 ≔ eigenvec ((A , B ((3)))) = ⎢ ⎥
⎢ 0.577 ⎥
⎣ 0.697 ⎦
⎡ -0.536 ⎤
⎢ -0.594 ⎥
Forma del segundo modo ϕ2 ≔ eigenvec ((A , B ((2)))) = ⎢ ⎥
⎢ 0.005 ⎥
⎣ 0.599 ⎦
⎡ -0.341 ⎤
⎢ -0.073 ⎥
Forma del tercer modo ϕ3 ≔ eigenvec ((A , B ((1)))) = ⎢ ⎥
⎢ 0.806 ⎥
⎣ -0.479 ⎦
⎡ -0.692 ⎤
⎢ 0.674 ⎥
Forma del cuarto modo ϕ4 ≔ eigenvec ((A , B ((0)))) = ⎢ ⎥
⎢ -0.251 ⎥
⎣ 0.064 ⎦
⎡ 0.232 -0.536 -0.341 -0.692 ⎤
⎢ 0.357 -0.594 -0.073 0.674 ⎥
Modos sin ϕ ≔ augment ⎛⎝ϕ1 , ϕ2 , ϕ3 , ϕ4⎞⎠ = ⎢ ⎥
normalizar ⎢ 0.577 0.005 0.806 -0.251 ⎥
⎣ 0.697 0.599 -0.479 0.064 ⎦
Factor de participación de masas
Normalización de modos con respecto a la masa
⎡ 0.032 ⎤ 1
―
Primer modo ϕ1 ⎢ 0.05 ⎥ cm 2
normalizado a la masa ϕ'1 ≔ ――――― =⎢ ⎥ ―――
ϕ1 T ⋅ m ⋅ ϕ1 ⎢ 0.081 ⎥
‾‾‾‾‾‾‾‾‾ 1
―
2
⎣ 0.097 ⎦ kg ⋅ s
⎡ -0.071 ⎤ 1
―
ϕ2 ⎢ -0.079 ⎥ cm
2
Segundo modo ϕ'2 ≔ ――――― =⎢ -4 ⎥ ―――
ϕ2 T ⋅ m ⋅ ϕ2 ⎢ 7.228 ⋅ 10 ⎥
‾‾‾‾‾‾‾‾‾ 1
normalizado a la masa ―
2
⎢⎣ 0.08 ⎥⎦ kg ⋅ s
⎡ -0.048 ⎤ 1
―
ϕ3 ⎢ -0.01 ⎥ cm 2
Tercer modo ϕ'3 ≔ ――――― =⎢ ⎥ ―――
ϕ3 T ⋅ m ⋅ ϕ3 ⎢ 0.114 ⎥
‾‾‾‾‾‾‾‾‾ 1
normalizado a la masa ―
2
⎣ -0.067 ⎦ kg ⋅ s
⎡ -0.089 ⎤ 1
―
ϕ4 ⎢ 0.087 ⎥ cm 2
Cuarto modo ϕ'4 ≔ ――――― =⎢ ⎥ ―――
‾‾‾‾‾‾‾‾‾
T -0.032 1
normalizado a la masa ϕ4 ⋅ m ⋅ ϕ4 ⎢ ⎥ ―
2
⎣ 0.008 ⎦ kg ⋅ s
Modos normalizados a la masa
⎡ 0.032 -0.071 -0.048 -0.089 ⎤ 1
―
⎢ 0.05 -0.079 -0.01 ⎥
0.087 cm
2
⎛ ⎞
ϕ ≔ augment ⎝ϕ'1 , ϕ'2 , ϕ'3 , ϕ'4⎠ = ⎢ -4 ⎥ ―――
⎢ 0.081 7.228 ⋅ 10 0.114 -0.032 ⎥ 1
―
2
⎢⎣ 0.097 0.08 -0.067 0.008 ⎥⎦ kg ⋅ s
d.2 Factor de participación modal
d.2 Factor de participación modal
⎡ 1 ⎤ ⎡ 13.757 ⎤ ―
1
⎢ 1 ⎥ ⎢ -5.249 ⎥ kg 2 ⋅ s
FPM ≔ ϕ T ⋅ m ⋅ ⎢ ⎥ = ⎢ ⎥ ―――
⎢ 1 ⎥ ⎢ -1.305 ⎥
1
―
2
⎣ 1 ⎦ ⎣ -1.317 ⎦ cm
d.3 Factor de participación de masas
2
FPM ((0))
Modo 1 fm1 ≔ ―――――――――――――――
2 2 2 2
= 85.93%
FPM ((0)) + FPM ((1)) + FPM ((2)) + FPM ((3))
2
FPM ((1))
Modo 2 fm1 ≔ ―――――――――――――――
2 2 2 2
= 12.51%
FPM ((0)) + FPM ((1)) + FPM ((2)) + FPM ((3))
2
FPM ((2))
Modo 3 fm1 ≔ ―――――――――――――――
2 2 2 2
= 0.773%
FPM ((0)) + FPM ((1)) + FPM ((2)) + FPM ((3))
2
FPM ((3))
Modo 4 fm1 ≔ ―――――――――――――――
2 2 2 2
= 0.787%
FPM ((0)) + FPM ((1)) + FPM ((2)) + FPM ((3))
d. Factor de participación de masas (Método alterno de teoría)
d.1 Factor de participación modal
⎡1⎤
⎢1⎥
ϕ1 T ⋅ m ⋅ ⎢ ⎥
Modo 1 ⎢1⎥
⎣1⎦
r1 ≔ ――――T
= 1.923
ϕ1 ⋅ m ⋅ ϕ1
⎡1⎤
⎢1⎥
ϕ2 T ⋅ m ⋅ ⎢ ⎥
Modo 2 ⎢1⎥
⎣1⎦
r2 ≔ ――――T
= -0.697
ϕ2 ⋅ m ⋅ ϕ2
⎡1⎤
⎢1⎥
ϕ3 T ⋅ m ⋅ ⎢ ⎥
Modo 3 ⎢1⎥
⎣1⎦
r3 ≔ ――――T
= -0.184
ϕ3 ⋅ m ⋅ ϕ3
⎡1⎤
⎢1⎥
ϕ4 T ⋅ m ⋅ ⎢ ⎥
Modo 4 ⎢1⎥
⎣1⎦
r4 ≔ ―――― = -0.169
ϕ4 T ⋅ m ⋅ ϕ4
d.2 Porcentaje de masa participativa
⎛ ⎡ 1 ⎤⎞ 2
⎜ T ⎢ 1 ⎥⎟
⎜ϕ1 ⋅ m ⋅ ⎢ ⎥⎟
Modo 1 ⎜ ⎢ 1 ⎥⎟
⎝ ⎣ 1 ⎦⎠ kg ⋅ s 2
M1 ≔ ――――― = 189.268 ―――
ϕ1 T ⋅ m ⋅ ϕ1 cm
⎛ ⎡ 1 ⎤⎞ 2
⎜ T ⎢ 1 ⎥⎟
⎜ϕ2 ⋅ m ⋅ ⎢ ⎥⎟
Modo 2 ⎜ ⎢ 1 ⎥⎟
⎝ ⎣ 1 ⎦⎠ kg ⋅ s 2
M2 ≔ ――――― = 27.554 ―――
ϕ2 T ⋅ m ⋅ ϕ2 cm
⎛ ⎡ 1 ⎤⎞ 2
⎜ T ⎢ 1 ⎥⎟
⎜ϕ3 ⋅ m ⋅ ⎢ ⎥⎟
Modo 3 ⎜ ⎢ 1 ⎥⎟
⎝ ⎣ 1 ⎦⎠ kg ⋅ s 2
M3 ≔ ――――― = 1.703 ―――
ϕ3 T ⋅ m ⋅ ϕ3 cm
⎛ ⎡ 1 ⎤⎞ 2
⎜ T ⎢ 1 ⎥⎟
⎜ϕ4 ⋅ m ⋅ ⎢ ⎥⎟
Modo 4 ⎜ ⎢ 1 ⎥⎟
⎝ ⎣ 1 ⎦⎠ kg ⋅ s 2
M4 ≔ ――――― = 1.734 ―――
ϕ4 T ⋅ m ⋅ ϕ4 cm
Suma de las masas (Debe dar la suma de las masas de toda la edificación)
Suma de las masas (Debe dar la suma de las masas de toda la edificación)
kg ⋅ s 2 kg ⋅ s 2
Mt ≔ m1 + m2 + m3 + m4 = 220.259 ――― M1234 ≔ M1 + M2 + M3 + M4 = 220.259 ―――
cm cm
d.3 Factor de participación de masas
M1
Modo 1 fm1 ≔ ―― = 85.93%
Mt
M2
Modo 2 fm2 ≔ ―― = 12.51%
Mt
M3
Modo 3 fm3 ≔ ―― = 0.773%
Mt
M4
Modo 4 fm4 ≔ ―― = 0.787%
Mt
* Son los mismos resultados por ambos métodos.
e. Pseudo aceleración y pseudo desplazamiento
Leemos la pseudo aceleración del espectro de respuesta correspondiente para el
periodo de cada modo.
Entradas
Salidas
Sa1 ≔ 0.12578125 ⋅ g Sa2 ≔ 0.12578125 ⋅ g Sa3 ≔ 0.12578125 ⋅ g Sa4 ≔ 0.12578125 ⋅ g
Como los 4 periodos resultan menores al Tp entonces son iguales
Modo Periodo Pseudo-aceleración Pseudo-desplazamiento
cm Sa1
Modo 1 T1 = 0.412 s Sa1 = 123.349 ―― Sd1 ≔ ―― = 0.53 cm
s2 ω1 2
cm Sa2
Modo 2 T2 = 0.172 s Sa2 = 123.349 ―― Sd2 ≔ ―― = 0.092 cm
s2 ω2 2
cm Sa3
Modo 3 T3 = 0.104 s Sa3 = 123.349 ―― Sd3 ≔ ―― = 0.034 cm
s2 ω3 2
cm Sa4
Modo 4 T4 = 0.077 s Sa4 = 123.349 ―― Sd4 ≔ ―― = 0.019 cm
s2 ω4 2
f. Aceleración y desplazamiento por modo
f. Aceleración y desplazamiento por modo
Recordando
⎡ 137.575 ⎤ ―
1
⎢ -52.492 ⎥ kg 2 ⋅ s
FPM = ⎢ ⎥ ―――
⎢ -13.052 ⎥
1
―
2
⎣ -13.167 ⎦ m
⎡ 55.077 ⎤ ⎡ 0.237 ⎤
⎢ 84.755 ⎥ cm a1 ⎢ 0.365 ⎥
Modo 1 a1 ≔ FPM ((0)) ⋅ Sa1 ⋅ ϕ'1 = ⎢ ⎥ ―― D1 ≔ ―― = ⎢ ⎥ cm
⎢ 136.841 ⎥ s ω1 2 ⎢ 0.589 ⎥
2
⎣ 165.341 ⎦ ⎣ 0.711 ⎦
⎡ 46.095 ⎤ ⎡ 0.034 ⎤
⎢ 51.036 ⎥ cm a2 ⎢ 0.038 ⎥
Modo 2 a2 ≔ FPM ((1)) ⋅ Sa2 ⋅ ϕ'2 = ⎢ ⎥ ―― D2 ≔ ―― = ⎢ -4 ⎥ cm
⎢ -0.468 ⎥ s 2
ω 2
2
⎢ -3.501 ⋅ 10 ⎥
⎣ -51.508 ⎦ ⎢⎣ -0.039 ⎥⎦
⎡ 7.719 ⎤ ⎡ 0.002 ⎤
⎢ 1.654 ⎥ cm a3 ⎢ 4.567 ⋅ 10 -4 ⎥
Modo 3 a3 ≔ FPM ((2)) ⋅ Sa3 ⋅ ϕ'3 = ⎢ ⎥ ―― D3 ≔ ―― =⎢ ⎥ cm
⎢ -18.273 ⎥ s ω3 2 ⎢ -0.005
2
⎥
⎣ 10.849 ⎦ ⎢⎣ 0.003 ⎥⎦
⎡ 14.458 ⎤ ⎡ 0.002 ⎤
⎢ -14.095 ⎥ cm a4 ⎢ -0.002 ⎥
Modo 4 a4 ≔ FPM ((3)) ⋅ Sa4 ⋅ ϕ'4 = ⎢ ⎥ ―― D4 ≔ ―― = ⎢ -4 ⎥ cm
⎢ 5.249 ⎥ s ω4 2 ⎢ 7.878 ⋅ 10 ⎥
2
⎣ -1.332 ⎦ ⎢⎣ -2 ⋅ 10 -4 ⎥⎦
Método 1
⎡ 3.37 ⎤
⎢ 5.186 ⎥
Fuerza inercial para el modo 1 f1 ≔ m ⋅ a1 = ⎢ ⎥ tonne
⎢ 6.698 ⎥
⎣ 8.093 ⎦
Cortante basal del modo 1 Vb1 ≔ f1 ((0)) + f1 ((1)) + f1 ((2)) + f1 ((3)) = 23.346 tonne
⎡ 2.82 ⎤
⎢ 3.123 ⎥
Fuerza inercial para el modo 2 f2 ≔ m ⋅ a2 = ⎢ ⎥ tonne
⎢ -0.023 ⎥
⎣ -2.521 ⎦
Cortante basal del modo 2 Vb2 ≔ f2 ((0)) + f2 ((1)) + f2 ((2)) + f2 ((3)) = 3.399 tonne
⎡ 0.472 ⎤
⎢ 0.101 ⎥
Fuerza inercial para el modo 3 f3 ≔ m ⋅ a3 = ⎢ ⎥ tonne
⎢ -0.894 ⎥
⎣ 0.531 ⎦
Cortante basal del modo 3 Vb3 ≔ f3 ((0)) + f3 ((1)) + f3 ((2)) + f3 ((3)) = 0.21 tonne
⎡ 0.885 ⎤
⎢ -0.862 ⎥
Fuerza inercial para el modo 3 f4 ≔ m ⋅ a4 = ⎢ ⎥ tonne
⎢ 0.257 ⎥
⎣ -0.065 ⎦
Cortante basal del modo 3 Vb4 ≔ f4 ((0)) + f4 ((1)) + f4 ((2)) + f4 ((3)) = 0.214 tonne
Método 2
f.1) Fuerzas inerciales
Opción a) Obtener las fuerzas mediante la rigidez y el desplazamiento
Desplazamientos de los 4 modos
⎡ 0.237 0.034 0.002 0.002 ⎤
⎢ 0.365 0.038 -4
4.567 ⋅ 10 -0.002 ⎥
D ≔ augment ⎛⎝D1 , D2 , D3 , D4⎞⎠ = ⎢ -4
⎥ cm
⎢ 0.589 -3.501 ⋅ 10 -0.005 7.878 ⋅ 10 -4 ⎥
⎢⎣ 0.711 -0.039 0.003 -2 ⋅ 10 -4 ⎥⎦
⎡ 3.37 2.82 0.472 0.885 ⎤
⎢ 5.186 3.123 0.101 -0.862 ⎥
Fuerzas inerciales de los 4 f≔k⋅D=⎢ ⎥ tonne
modos ⎢ 6.698 -0.023 -0.894 0.257 ⎥
⎣ 8.093 -2.521 0.531 -0.065 ⎦
Opción b) Obtener las fuerzas mediante la masa y la aceleración
⎡ 55.077 46.095 7.719 14.458 ⎤
⎢ 84.755 51.036 1.654 -14.095 ⎥ cm
Aceleraciones de los 3 a ≔ augment ⎛⎝a1 , a2 , a3 , a4⎞⎠ = ⎢ ⎥ ――
⎢ 136.841 -0.468 -18.273 5.249 ⎥ s 2
modos ⎣ 165.341 -51.508 10.849 -1.332 ⎦
⎡ 3.37 2.82 0.472 0.885 ⎤
⎢ 5.186 3.123 0.101 -0.862 ⎥
Fuerzas inerciales de los 4 f≔m⋅a=⎢ ⎥ tonne
modos ⎢ 6.698 -0.023 -0.894 0.257 ⎥
⎣ 8.093 -2.521 0.531 -0.065 ⎦
f.2) Cortante basal
f.2) Cortante basal
Matriz diagonal de participación modal
⎡ FPM ((0)) 0 0 0 ⎤ ⎡ 13.757 0 0 0 ⎤ ―
1
⎢ 0 FPM ((1)) 0 0 ⎥ ⎢ 0 -5.249 0 0 ⎥ kg 2 ⋅ s
γ≔⎢ =
⎥ ⎢ ⎥ ―――
⎢ 0 0 FPM ((2)) 0 ⎥ ⎢ 0 0 -1.305 0 ⎥
1
―
⎣ 0 0 0 FPM ((3)) ⎦ ⎣ 0 0 0 -1.317 ⎦ cm 2
Matriz diagonal de aceleraciones espectrales
⎡ Sa1 0 0 0 ⎤ ⎡ 123.349 0 0 0 ⎤
⎢ ⎥ ⎢ ⎥ cm
0 Sa2 0 0
Sa ≔ ⎢ ⎥=⎢ 0 123.349 0 0
⎥ ――
⎢ 0 0 Sa3 0 ⎥ ⎢ 0 0 123.349 0 ⎥ s
2
⎢⎣ 0 0 0 ⎥
Sa3 ⎦ ⎣ 0 0 0 123.349 ⎦
Cortante basal
⎡ 23.346 0 0 0 ⎤
⎢ 0 3.399 0 0 ⎥
V ≔ γ 2 ⋅ Sa = ⎢ ⎥ tonne
⎢ 0 0 0.21 0 ⎥
⎣ 0 0 0 0.214 ⎦
Cortantes de entrepiso
(Desplazamiento del piso - Desplazamiento del piso inferior)/Altura del piso
⎡ f1 ((0)) + f1 ((1)) + f1 ((2)) + f1 ((3)) ⎤ ⎡ 23.346 ⎤
⎢ ⎥ ⎢ ⎥
f1 ((1)) + f1 ((2)) + f1 ((3))
Fuerza cortante de los Vent1 ≔ ⎢ ⎥ = ⎢ 19.976 ⎥ tonne
⎢ f1 ((2)) + f1 ((3)) ⎥ ⎢ 14.791 ⎥
entrepisos modo 1
⎢⎣ f1 ((3)) ⎥⎦ ⎣ 8.093 ⎦
⎡ f2 ((0)) + f2 ((1)) + f2 ((2)) + f2 ((3)) ⎤ ⎡ 3.399 ⎤
⎢ ⎥ ⎢ ⎥
f2 ((1)) + f2 ((2)) + f2 ((3))
Fuerza cortante de los Vent2 ≔ ⎢ ⎥ = ⎢ 0.578 ⎥ tonne
⎢ f2 ((2)) + f2 ((3)) ⎥ ⎢ -2.544 ⎥
entrepisos modo 2
⎢⎣ f2 ((3)) ⎥⎦ ⎣ -2.521 ⎦
⎡ f3 ((0)) + f3 ((1)) + f3 ((2)) + f3 ((3)) ⎤ ⎡ 0.21 ⎤
⎢ ⎥ ⎢ ⎥
f3 ((1)) + f3 ((2)) + f3 ((3))
Fuerza cortante de los Vent3 ≔ ⎢ ⎥ = ⎢ -0.262 ⎥ tonne
⎢ f3 ((2)) + f3 ((3)) ⎥ ⎢ -0.363 ⎥
entrepisos modo 3
⎢⎣ (
f3 (3)) ⎥⎦ ⎣ 0.531 ⎦
⎡ f4 ((0)) + f4 ((1)) + f4 ((2)) + f4 ((3)) ⎤ ⎡ 0.214 ⎤
⎢ ⎥ ⎢ ⎥
f4 ((1)) + f4 ((2)) + f4 ((3))
Fuerza cortante de los Vent4 ≔ ⎢ ⎥ = ⎢ -0.671 ⎥ tonne
⎢ f4 ((2)) + f4 ((3)) ⎥ ⎢ 0.192 ⎥
entrepisos modo 4
⎢⎣ f4 ((3)) ⎥⎦ ⎣ -0.065 ⎦
Derivas (Desplazamiento del piso - Desplazamiento del piso inferior)/Altura del piso
Derivas (Desplazamiento del piso - Desplazamiento del piso inferior)/Altura del piso
⎡ D1 ((0)) ⎤
⎢ ――― ⎥
⎢ H1 ⎥
⎢ D1 ((1)) - D1 ((0)) ⎥ ⎡ 6.768 ⋅ 10 -4 ⎤
⎢ ――――― ⎥ ⎢ -4 ⎥
⎢ H2 ⎥ = ⎢ 4.255 ⋅ 10 ⎥
Deriva de los pisos modo 1 Δ1 ≔
⎢ D1 ((2)) - D1 ((1)) ⎥ ⎢ 7.467 ⋅ 10 -4 ⎥
⎢ ――――― ⎥ ⎢ -4 ⎥
⎢ H3 ⎥ ⎣ 4.086 ⋅ 10 ⎦
⎢ D1 ((3)) - D1 ((2)) ⎥
⎢ ――――― ⎥
⎣ H4 ⎦
⎡ D2 ((0)) ⎤
⎢ ――― ⎥
⎢ H1 ⎥
⎢ D2 ((1)) - D2 ((0)) ⎥ ⎡ 9.853 ⋅ 10 -5 ⎤
⎢ ――――― ⎥ ⎢ ⎥
⎢ H2 ⎥ 1.232 ⋅ 10 -5 ⎥
Deriva de los pisos modo 2 Δ2 ≔ = ⎢
⎢ D2 ((2)) - D2 ((1)) ⎥ ⎢ -1.284 ⋅ 10 -4 ⎥
⎢ ――――― ⎥ ⎢ -4 ⎥
⎢ H3 ⎥ ⎣ -1.273 ⋅ 10 ⎦
⎢ D2 ((3)) - D2 ((2)) ⎥
⎢ ――――― ⎥
⎣ H4 ⎦
⎡ D3 ((0)) ⎤
⎢ ――― ⎥
⎢ H1 ⎥
⎢ D3 ((1)) - D3 ((0)) ⎥ ⎡ 6.091 ⋅ 10 -6 ⎤
⎢ ――――― ⎥ ⎢ -6 ⎥
⎢ H2 ⎥ = ⎢ -5.584 ⋅ 10 ⎥
Deriva de los pisos modo 3 Δ3 ≔
⎢ D3 ((2)) - D3 ((1)) ⎥ ⎢ -1.834 ⋅ 10 -5 ⎥
⎢ ――――― ⎥ ⎢ -5 ⎥
⎢ H3 ⎥ ⎣ 2.681 ⋅ 10 ⎦
⎢ D3 ((3)) - D3 ((2)) ⎥
⎢ ――――― ⎥
⎣ H4 ⎦
⎡ D4 ((0)) ⎤
⎢ ――― ⎥
⎢ H1 ⎥
⎢ D4 ((1)) - D4 ((0)) ⎥ ⎡ 6.2 ⋅ 10 -6 ⎤
⎢ ――――― ⎥ ⎢ -5 ⎥
H2 ⎥ = ⎢ -1.428 ⋅ 10 ⎥
Deriva de los pisos modo 4 Δ4 ≔ ⎢
⎢ D4 ((2)) - D4 ((1)) ⎥ ⎢ 9.678 ⋅ 10 -6 ⎥
⎢ ――――― ⎥ ⎢ -6 ⎥
⎢ H3 ⎥ ⎣ -3.293 ⋅ 10 ⎦
⎢ D4 ((3)) - D4 ((2)) ⎥
⎢ ――――― ⎥
⎣ H4 ⎦
cm
a1 ((3)) = 165.341 ――
s2
cm
f1 ((3)) = 8.093 tonne a1 ((2)) = 136.841 ――
f1 ((2)) = 6.698 tonne s2
cm
f1 ((1)) = 5.186 tonne a1 ((1)) = 84.755 ――
f1 ((0)) = 3.37 tonne s2
cm
a1 ((0)) = 55.077 ――
Vb1 = 23.346 tonne s2
cm
a2 ((3)) = -51.508 ――
s2
cm
f2 ((3)) = -2.521 tonne a2 ((2)) = -0.468 ――
f2 ((2)) = -0.023 tonne s2
cm
f2 ((1)) = 3.123 tonne a2 ((1)) = 51.036 ――
f2 ((0)) = 2.82 tonne s2
cm
a2 ((0)) = 46.095 ――
Vb2 = 3.399 tonne s2
cm
a3 ((3)) = 10.849 ――
s2
cm
f3 ((3)) = 0.531 tonne a3 ((2)) = -18.273 ――
f3 ((2)) = -0.894 tonne s2
cm
f3 ((1)) = 0.101 tonne a3 ((1)) = 1.654 ――
f3 ((0)) = 0.472 tonne s2
cm
a3 ((0)) = 7.719 ――
Vb3 = 0.21 tonne s2
cm
a4 ((3)) = -1.332 ――
s2
cm
f4 ((3)) = -0.065 tonne a4 ((2)) = 5.249 ――
f4 ((2)) = 0.257 tonne s2
cm
f4 ((1)) = -0.862 tonne a4 ((1)) = -14.095 ――
f4 ((0)) = 0.885 tonne s2
cm
a4 ((0)) = 14.458 ――
Vb4 = 0.214 tonne s2
Criterios de combinación
Criterios de combinación
Método SRSS (Square root of the sum of the squares = raíz cuadrada de la suma de
cuadrados)
⎡ 0.239 ⎤
⎢ 0.367 ⎥
Desplazamientos SRSS DSRSS ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
D1 2 + D2 2 + D3 2 + D4 2 = ⎢ ⎥ cm
⎢ 0.589 ⎥
⎣ 0.712 ⎦
Aceleración SRSS
⎡ 73.667 ⎤ ⎡ 0.075 ⎤
⎢ ⎥
99.947 cm ⎢ 0.102 ⎥
aSRSS ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
a1 2 + a2 2 + a3 2 + a4 2 = ⎢ ⎥ ―― a SRSS = ⎢ ⎥g
⎢ 138.156 ⎥ s ⎢ 0.141 ⎥
2
⎣ 173.523 ⎦ ⎣ 0.177 ⎦
Cortante basal SRSS VbSRSS ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
Vb1 2 + Vb2 2 + Vb3 2 + Vb4 2 = 23.594 tonne
⎡ 4.507 ⎤
⎢ 6.115 ⎥
Fuerzas inerciales SRSS fSRSS ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
f1 2 + f2 2 + f3 2 + f4 2 = ⎢ ⎥ tonne
⎢ 6.762 ⎥
⎣ 8.493 ⎦
Cortantes de entrepiso SRSS
⎡ 23.594 ⎤
⎢ 19.998 ⎥
VentSRSS ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
Vent1 2 + Vent2 2 + Vent3 2 + Vent4 2 = ⎢ ⎥ tonne
⎢ 15.014 ⎥
⎣ 8.493 ⎦
Derivas por SRSS
⎡ 6.84 ⋅ 10 -4 ⎤
⎢ ⎥
4.259 ⋅ 10 -4 ⎥
ΔSRSS ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
Δ1 2 + Δ2 2 + Δ3 2 + Δ4 2 = ⎢
⎢ 7.58 ⋅ 10 -4 ⎥
⎢⎣ 4.288 ⋅ 10 -4 ⎥⎦
Metodo de combinacion absoluta
Desplazamientos son:
⎡ 0.237 ⎤ ⎡ 0.034 ⎤ ⎡ 0.002 ⎤ ⎡ 0.002 ⎤
⎢ 0.365 ⎥ ⎢ 0.038 ⎥ ⎢ -4 ⎥ ⎢ -0.002 ⎥
4.567 ⋅ 10 ⎢
D1 = ⎢ ⎥ cm D2 = ⎢ -4 ⎥ cm D 3 = ⎢ ⎥ cm D 4 = -4 ⎥ cm
⎢ 0.589 ⎥ ⎢ -3.501 ⋅ 10 ⎥ ⎢ -0.005 ⎥ ⎢ 7.878 ⋅ 10 ⎥
⎣ 0.711 ⎦ ⎢⎣ -0.039 ⎥⎦ ⎢⎣ 0.003 ⎥⎦ ⎢⎣ -2 ⋅ 10 -4 ⎥⎦
convertir en valores absolutos
convertir en valores absolutos
⎡ 0.237 ⎤ ⎡ 0.034 ⎤
⎢ 0.365 ⎥ ⎢ 0.038 ⎥
D1abs ≔ ‾‾‾‾
D1 2 = ⎢ ⎥ cm D2abs ≔ ‾‾‾‾
D2 2 = ⎢ -4 ⎥ cm
⎢ 0.589 ⎥ ⎢ 3.501 ⋅ 10 ⎥
⎣ 0.711 ⎦ ⎢⎣ 0.039 ⎥⎦
⎡ 0.002 ⎤ ⎡ 0.002 ⎤
⎢ 4.567 ⋅ 10 -4 ⎥ ⎢ 0.002 ⎥
D3abs ≔ ‾‾‾‾
D3 2 = ⎢ ⎥ cm D4abs ≔ ‾‾‾‾
D4 2 = ⎢ -4 ⎥ cm
⎢ 0.005 ⎥ ⎢ 7.878 ⋅ 10 ⎥
⎢⎣ 0.003 ⎥⎦ ⎢⎣ 2 ⋅ 10 -4 ⎥⎦
⎡ 0.276 ⎤
⎢ 0.405 ⎥
DABS ≔ D1abs + D2abs + D3abs + D4abs = ⎢ ⎥ cm
⎢ 0.595 ⎥
⎣ 0.753 ⎦
Respuesta máxima alternativa de la norma E 030
⎡ 0.248 ⎤
⎢ 0.376 ⎥
DE_030 ≔ 0.25 ⋅ DABS + 0.75 ⋅ DSRSS = ⎢ ⎥ cm
⎢ 0.59 ⎥
⎣ 0.722 ⎦
b) Proponga un cambio de dimensiones para reducir las deflexiones al 50% de los
encontrados en el punto a)
Con mayor rigidez hay menor desplazamiento y para aumentar la rigidez se puede
ampliar las dimensiones de las columnas
Con las siguientes dimensiones:
Obtenemos:
Obtenemos:
Entonces habiendo reducido los desplazamientos alrededor del 50% las dimensiones
propuestas son de columnas de 40 cm x 70 cm para los pisos 1-2 y de 40 cm x 55
cm para los pisos 3-4.
EJERCICIO 2:
1) Obtención de la matriz de rigidez
⎡ 125 -41.66666667 1250 1250 -1250 -1250 ⎤
⎢ -41.666667 41.66666667 1250 1250 1250 1250 ⎥
⎢ ⎥
1250 1250 400000 50000 50000 0 ⎥
K≔⎢
⎢ 1250 1250 50000 400000 0 50000 ⎥
⎢ 1250 1250 50000 0 150000 25000 ⎥
⎢⎣ 1250 1250 0 50000 25000 150000 ⎥⎦
Condensación estática
⎡ 125 -41.66666667 ⎤
ktt ≔ ⎢
⎣ -41.66666667 41.66666667 ⎥⎦
⎡ 1250 1250 ⎤
⎢ 1250 1250 ⎥
kto ((T)) ≔ ⎢ ⎥
⎢ -1250 1250 ⎥
⎣ -1250 1250 ⎦
⎡ 1250 1250 -1250 -1250 ⎤
ktto ≔ ⎢
⎣ 1250 1250 1250 1250 ⎥⎦
⎡ 400000 50000 50000 0 ⎤
⎢ 50000 400000 0 50000 ⎥
koo ≔ ⎢ ⎥
⎢ 50000 0 150000 25000 ⎥
⎣ 0 50000 25000 150000 ⎦
⎡ 2.66269E-06 -3.67611E-07 -9.33929E-07 2.78192E-07 ⎤
⎢ -3.67611E-07 2.66269E-06 2.78192E-07 -9.33929E-07 ⎥
koo - 1 ⎢ ⎥
⎢ -9.33929E-07 2.78192E-07 7.19324E-06 -1.2916E-06 ⎥
⎣ 2.78192E-07 -9.33929E-07 -1.2916E-06 7.19324E-06 ⎦
⎡ 0.003688525 0.00204918 ⎤
⎢ 0.003688525 0.00204918 ⎥
koo - 1 ⋅ kto ⎢ ⎥
⎢ -0.008196721 0.006557377 ⎥
⎣ -0.008196721 0.006557377 ⎦
⎡ 29.71311475 -11.2704918 ⎤
kto ⋅ koo - 1 ⋅ kto ⎢
⎣ -11.2704918 21.51639344 ⎥⎦
⎡ 95.28688525 -30.39617486 ⎤
k≔⎢
⎣ -30.39617486 20.15027322 ⎥⎦
2) Matriz de masa
⎡ 1.03626943 0 ⎤ ⎡ 0.965 0 ⎤
M≔⎢ Minversa ≔ ⎢
⎣ 0 0.518134715 ⎥⎦ ⎣ 0 1.93 ⎥⎦
3) Modos de vibración
⎡ 91.95184426 -29.33230874 ⎤
Minversa ⋅ ks ⎢
⎣ -58.66461749 38.89002732 ⎥⎦
Resolviendo con HP se obtiene:
⎛ rad ⎞ 2 ⎛ rad ⎞ 2
w12 ≔ 16.18 ⎜―― ⎟ w1 ≔ 4.022 ⎜―― ⎟
⎝ s ⎠ ⎝ s ⎠
⎛ rad ⎞ 2 ⎛ rad ⎞ 2
w22 ≔ 114.662 ⎜―― ⎟ w2 ≔ 10.70 ⎜―― ⎟
⎝ s ⎠ ⎝ s ⎠
⎡ 0.3871 1 ⎤
Matriz φ: ⎢⎣ 1 -0.7742 ⎥⎦
4) Calculamos las aceleraciones espectrales para cada modo de vibración
5) Calculamos la contribución de cada modo y la matriz de fuerzas laterales
Masa generalizada para cada modo (Mn):
⎡ 0.3871 ⎤ ⎡ 1 ⎤
ϕ1 ≔ ⎢ ϕ2 ≔ ⎢
⎣ 1 ⎥⎦ ⎣ -0.7742 ⎥⎦
M1 ≔ 0.673415969 M2 ≔ 1.346831938
Calculamos Ln para cada modo: Calculamos el aporte de cada modo:
r_1 ≔ 1.365091791
Ln_1 ≔ 0.919274611
Ln_2 ≔ 0.635129534 r_2 ≔ 0.471572968
Obtenemos las matrices S para cada modo:
⎡ 0.54759278 ⎤ ⎡ 0.48867665 ⎤
S1 ≔ ⎢ S2 ≔ ⎢
⎣ 0.707301446 ⎥⎦ ⎣ -0.189166731 ⎥⎦
En cada piso la sumatoria de los "s" de cada
modo debe dar la masa de cada nivel
1.03626943 OK!!
0.518134715 OK!!
Si sumamos todos los valores de las dos matrices "S",
obtendremos la masa total de la estructura
S total 1.55440415
Para obtener la masa participativa de cada modo, obtendremos la
sumatoria de cada matriz "s" y la dividimos entre la masa total
%S1 = 80.7315285%
% S2 = 19.2684715%
Calculamos las fuerzas laterales efectivas para cada modo de vibración
⎡ 18.17661841 ⎤ ⎡ 42.18678769 ⎤
f1 ≔ ⎢ f2 ≔ ⎢
⎣ 23.47793646 ⎥⎦ ⎣ -16.33050551 ⎥⎦
6) Calculamos los desplazamientos laterales y distorsiones de entrepiso
⎡ 18.17661841 ⎤ ⎡ 42.18678769 ⎤
f1 ≔ ⎢ f2 ≔ ⎢
⎣ 23.47793646 ⎥⎦ ⎣ -16.33050551 ⎥⎦
⎡ 1.084099485 ⎤ ⎡ 0.355064737 ⎤
U1 ≔ ⎢ U2 ≔ ⎢
⎣ 2.800478849 ⎥⎦ ⎣ -0.274829806 ⎥⎦
Los desplazamientos con la combinación SRSS son:
⎡ 1.140764069 ⎤ ⎡ 0.009506367 ⎤
UT ≔ ⎢ drift ≔ ⎢
⎣ 1.828312134 ⎥⎦ ⎣ 0.015235934 ⎥⎦
