Anae Grupo3 TALLER3

ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO
FACULTAD DE INGENIERÍA CIVIL
ANÁLISIS ESTRUCTURAL
TALLER 3.
INTEGRANTES:
GUTIÉRREZ MARTÍNEZ CHRISTIAN
PENAGOS ZÁRATE LUISA
PORTILLO SÁNCHEZ JOSÉ
PRIETO MORENO SEBASTIÁN
RODRIGUEZ BUENO NATALIA
PRESENTADO A:
ING. HECTOR JOSÉ PEREZ BARRERA

INGENIERÍA CIVIL
BOGOTÁ, MAYO 21 DE 2021
1
1. INFORMACIÓN DE ENTRADA
La estructura de cuatro pisos mostrada abajo tiene alturas de piso típicas de 3.0m, y distancia
entre ejes de columnas de 5.0m. Las vigas son de sección 0.4m x 0.5m, y las columnas de
0.5m x 0.5m. El concreto de todos los elementos es de f’c=21 MPa, con módulo de elasticidad
Ec=21500 MPa. Las masas de cada piso se indican en la figura abajo.
2. CARACTERIZACIÓN DÍNAMICA
a. Matriz de masa del sistema
48.8 0 0 0 𝑃𝑖𝑠𝑜 1
𝑀= [ 0 48.8 0 0 ] 𝑇𝑜𝑛 𝑃𝑖𝑠𝑜 2
0 0 48.8 0 𝑃𝑖𝑠𝑜 3
0 0 0 37.5 𝑃𝑖𝑠𝑜 4
b. Matriz de rigidez lateral usando el método 2 visto en clase.
Para emplear dicho método, se plantea la aplicación de cargas laterales a la estructura para
así obtener los desplazamientos y reacciones. Dicho proceso se repetirá para cada uno de los
pisos para así obtener la matriz de rigidez de la estructura. Se sabe que la rigidez se calcula
2
como
𝐹
𝑘=
𝑑
Donde
✓ F corresponde a las cargas aplicadas o reacciones [kN]
✓ d: Desplazamiento
Piso 1
Del programa se obtuvo que el desplazamiento para el piso 1 es de 0.0003106 m.
∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = 𝐹11 → ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = 100𝑘𝑁
∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎2 = 𝐹21 + 𝐹22 + 𝐹32 + 𝐹42
→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 = −14.86𝑘𝑁 − 14.2𝑘𝑁 − 14.42𝑘𝑁 − 13.68𝑘𝑁
∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 = 𝐹31 + 𝐹32 + 𝐹33 + 𝐹34
→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 = 5.52𝑘𝑁 + 2.68𝑘𝑁 + 3.3𝑘𝑁 + 2.42𝑘𝑁
3
→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4
= −4.67𝑘𝑁 + 0.0886𝑘𝑁 − 0.758𝑘𝑁 + 3.56𝑘𝑁
REACCIONES PISO 1
ΣReacciones1 – R1 (kN) 100
ΣReacciones2 – R2 (kN) -57.16
ΣReacciones3 – R3 (kN) 13.92
ΣReacciones4 – R4 (kN) -1,7815
𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 100

𝑘11 = → 𝑘11 = → 𝑘11 = 321957.5016
𝑑 0.0003106
𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 −57.16

𝑘21 = → 𝑘21 = → 𝑘21 = −184030.9079
𝑑 0.0003106
𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 13.92

𝑘31 = → 𝑘31 = → 𝑘31 = −184030.9079
𝑑 0.0003106

𝑘41 = → 𝑘41 = → 𝑘41 = −5735.6729
𝑑 0.0003106
Piso 2
4
= −15.06𝑘𝑁 − 16.65𝑘𝑁 − 16.47𝑘𝑁 − 17.96𝑘𝑁
→ ∑ 𝑅𝑒𝑎𝑐𝑐𝑖𝑜𝑛𝑒𝑠3 = −13.66𝑘𝑁 − 16.31𝑘𝑁 − 15.67𝑘𝑁 − 16.73𝑘𝑁
∑ 𝐶𝑎𝑟𝑔𝑎𝑠 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = 𝐹41 + 𝐹42 + 𝐹43 + 𝐹44
→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = −0.72𝑘𝑁 + 3.33𝑘𝑁 + 2.48𝑘𝑁 + 7.51𝑘𝑁
REACCIONES PISO 2
ΣReacciones2 – R2 (kN) -100

𝑘12 = → 𝑘12 = → 𝑘12 = −184028.937
𝑑 0.0003594
𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 −100

𝑘22 = → 𝑘22 = → 𝑘22 = 1112966.05
𝑑 0.0003594

𝑘32 = → 𝑘32 = → 𝑘32 = 173539.232
𝑑 0.0003594

𝑘42 = → 𝑘42 = → 𝑘42 = 35058.4307
𝑑 0.0003594
Piso 3
5
= − 17.08𝑘𝑁 − 17.25𝑘𝑁 − 17.51𝑘𝑁 − 15.96𝑘𝑁
REACCIONES PISO 3
6

𝑘13 = → 𝑘13 = → 𝑘13 = 44828.469
𝑑 0.0003906

𝑘23 = → 𝑘23 = → 𝑘23 = −173579.1091
𝑑 0.0003906

𝑘33 = → 𝑘33 = → 𝑘33 = 1024065.54
𝑑 0.0003906

𝑘43 = → 𝑘43 = → 𝑘43 = −120250.8961
𝑑 0.0003906
Piso 4
Por último, en el Piso 4 en el programa se obtuvo que el desplazamiento es de 0.00111 m.
7
= −27.94𝑘𝑁 − 37.53𝑘𝑁 − 36.99𝑘𝑁 − 31.13𝑘𝑁
REACCIONES PISO 4

𝑘14 = → 𝑘14 = → 𝑘14 = −5765.76577
𝑑 0.00111

𝑘24 = → 𝑘24 = → 𝑘24 = 35099.0991
𝑑 0.00111

𝑘34 = → 𝑘34 = → 𝑘34 = −120351.351
𝑑 0.00111
𝑘44 = → 𝑘44 = → 𝑘44 = 360360.36
𝑑 0.00111
Con los valores de rigidez calculados para cada piso se arma la matriz de rigidez del sistema
𝑘11 𝑘12 𝑘13 𝑘14

𝑘21 𝑘22 𝑘23 𝑘24
𝑘= [ ]
𝑘31 𝑘32 𝑘33 𝑘34
𝑘41 𝑘42 𝑘43 𝑘44
8
321957.5 −184028.937 44828.469 −5765.7658

𝑘 = [−184030.908 278241.514 −173579.109 35099.0991 ] 𝑘𝑁
44816.4842 −173539.232 256016.39 −120351.351 𝑚
−5735.6729 35058.4307 −120250.896 90090.09
c. Matriz de rigidez lateral usando el método 3 (edificio cortante) visto en clase.
Teniendo en cuenta lo siguiente
9
𝐹11 𝐹12 𝐹13 𝐹14

𝐹 𝐹22 𝐹23 𝐹24
𝑘 = [ 21 ]
𝐹31 𝐹32 𝐹33 𝐹34
𝐹41 𝐹42 𝐹43 𝐹44
𝐾1 + 𝐾2 −𝐾2 0 0
−𝐾2 𝐾2 + 𝐾3 −𝐾3 0
𝑘= [ ]
0 −𝐾3 𝐾3 + 𝐾4 −𝐾4
0 0 −𝐾4 𝐾4
Se obtiene la matriz de rigidez
398148 −199074 0 0
𝑘= [ −199074 398148 −199074 0 ]
0 −199074 398148 −199074
0 0 −199074 199074
d. Periodos de vibración de todos los modos de vibración de la estructura, usando

la matriz de rigidez obtenida en (b).
321957.5 −184028.937 44828.469 −5765.7658
𝑘 = [−184030.908 278241.514 −173579.109 35099.0991 ] 𝑘𝑁
44816.4842 −173539.232 256016.39 −120351.351 𝑚
−5735.6729 35058.4307 −120250.896 90090.09
Mediante Microsoft Excel se obtiene las frecuencias de cada modo de vibración

mediante el comando solver, teniendo así:
➢ Modo 1
𝑤1 = 13.31627 𝑟𝑎𝑑/𝑠
➢ Modo 2
𝑤2 = 42.950524 𝑟𝑎𝑑/ 𝑠
➢ Modo 3
𝑤3 = 77.277827 𝑟𝑎𝑑/𝑠
➢ Modo 4
𝑤4 = 109,33374 𝑟𝑎𝑑/𝑠
10
Para obtener los periodos de vibración de cada modo se utiliza la siguiente ecuación
2𝜋
𝑇=
𝑤
➢ Modo 1:
2𝜋
𝑇1 = = 0.4718 𝑠
𝑤1
➢ Modo 2:
2𝜋
𝑇2 = = 0.1462 𝑠
𝑤2
➢ Modo 3:
2𝜋
𝑇3 = = 0.0813 𝑠
𝑤3
➢ Modo 4:
2𝜋
𝑇4 = = 0.05746 𝑠
𝑤4
e. Periodos de vibración de todos los modos de vibración de la estructura, usando

la matriz de rigidez obtenida en (c). Comparar los resultados, justificar y
concluir sobre las diferencias.
Mediante el aplicativo Wolfram Mathematica se obtiene las frecuencias de cada

modo de vibración.
11
➢ Modo 1
𝑤1 = 23.3601𝑟𝑎𝑑/𝑠
➢ Modo 2
𝑤2 = 66.6976𝑟𝑎𝑑/ 𝑠
➢ Modo 3
𝑤3 = 100.5845 𝑟𝑎𝑑/𝑠
➢ Modo 4
𝑤4 = 121.1339 𝑟𝑎𝑑/𝑠
Para obtener los periodos de vibración de cada modo se utiliza la siguiente ecuación
2𝜋
𝑇=
𝑤
➢ Modo 1:
2𝜋
𝑇1 = = 0.2689𝑠
𝑤1
➢ Modo 2:
2𝜋
𝑇2 = = 0.0942𝑠
𝑤2
➢ Modo 3:
2𝜋
𝑇3 = = 0.0624𝑠
𝑤3
➢ Modo 4:
12
2𝜋
𝑇4 = = 0.0518 𝑠
𝑤4
Los valores de los periodos del método 3 y el método 2 corresponden a lo esperado

ya que los periodos del método 3 son mucho menores que los del método 2 porque en
el método 3 se tiene en cuenta restricciones de tercer grado lo que hace más rígido la
estructura y se obtiene periodos de vibración menores y en el método dos se tienen
restricciones de segundo grado para hallar la rigidez lo que hace que los periodos sean
más extensos.
f. Usando un análisis de valores propios (eigen-análisis), junto con la matriz de rigidez

del literal (b), obtener frecuencias y periodos de vibración, y vectores de amplitudes
modales sin normalizar.
Para el cálculo de los eigenvalores son necesarias las matrices de masa y rigidez que se
determinaron en el primer punto. Teniendo en cuenta ambas matrices se emplea la siguiente
ecuación:
det([𝐾] − 𝜔𝑖2 [𝑀]) = 0
Se resuelve el determinante de la matriz dejando expresada la frecuencia angular ὡ y al

resolver el determinante se obtiene una ecuación que debe ser igualada a cero para obtener
los respectivos vectores de amplitudes sin normalizar para cada modo.
w1 13,3162703 rad/s T1 0,47184273

w2 42,9505245 rad/s T2 0,14628891
w3 77,2778274 rad/s T3 0,08130644
w4 109,333741 rad/s T4 0,05746794
13
Vector de amplitud
modal sin normalizar
- Modo 1
1
2,444603191
3,588469961
4,213164501
Vector de amplitud
- Modo 2
1
1,32743064
0,094415211
-1,408204306
Vector de amplitud
- Modo 3
1
-0,077516197
-0,903049542
0,748118701
14
Vector de amplitud
- Modo 4
1
-1,180929259
0,926075493
-0,442511149
g. Normalizar los modos resultantes con respecto a la matriz de masa.
Con los vectores de amplitud modal sin normalizar para cada modo hallados en el numeral
inmediatamente anterior y con la matriz de masa, se determina el factor normalizado, y con
este los vectores de amplitudes normalizados. Para lo anterior se emplea la siguiente
ecuación:
✓ Modo 1
Factor
40,42882074
Normalizado
Modo 1 Modo normalizado

1 0,02473483
2,444603191 0,060466844
3,588469961 0,088760194
4,213164501 0,104211907
✓ Modo 2
Factor
14,4771582
Normalizado

1 0,06907433
1,32743064 0,091691382
0,094415211 0,006521667
-1,408204306 -0,097270769
✓ Modo 3
15
Factor
10,48225227
Normalizado

1 0,095399345
-0,077516197 -0,007394994
-0,903049542 -0,086150335
0,748118701 0,071370034
✓ Modo 4
Factor
12,88607541
Normalizado

1 0,077603147
-1,180929259 -0,091643827
0,926075493 0,071866372
-0,442511149 -0,034340258
h. Calcular coeficientes de participación modal de cada uno de los modos.
Para calcular los coeficientes de participación modal, se trasponer la matriz de modos de cada
modo y después con el vector de masa se determinan los coeficientes α:
Modo de vibración 1:
48,8 0 0 0 1
0 48,8 0 0 1
[0,0247 0,0604 0,0887 0,10421] * [ ]∗ [ ] = 12,397
0 0 48,8 0 1
0 0 0 37,5 1
48,8 0 0 0 1
0 48,8 0 0 1
[0,0690 0,0916 0,0065 −0,09727] * [ ]∗ [ ] = 4,515
0 0 48,8 0 1
0 0 0 37,5 1
48,8 0 0 0 1
0 48,8 0 0 1
[0,0953 −0,0073 −0,0861 0,0713] * [ ]∗ [ ] =2,7668
0 0 48,8 0 1
0 0 0 37,5 1
16
48,8 0 0 0 1
0 48,8 0 0 1
[0,0776 −0,0916 0,07186 −0,0343] * [ ]∗ [ ] = 1,534
0 0 48,8 0 1
0 0 0 37,5 1
Teniendo así que los coeficientes de participación modal son:
Coef. Participación
modal
∝1 12,397
∝2 4,516
∝3 2,767
∝4 1,534
i. Calcular masas efectivas de cada modo de vibración, y comentar sobre los

resultados.
Para calcular las masas efectivas de cada modo de vibración, se toman los coeficientes de
participación modal obtenidos en el numeral inmediatamente anteriormente en donde su
resultado corresponde a elevar al cuadrado dichos coeficientes. Para determinar el porcentaje
de masa para cada uno de los modos, se determina la relación entre el coeficiente de
participación al cuadrado de cada modo entre la sumatoria de los coeficientes al cuadrado
teniendo así:
Modo α α^2(Ton) % Masa

1 12,3973 153,6927 83,485
2 4,5160 20,3940 11,078
3 2,7669 7,6555 4,158
4 1,5341 2,3536 1,278
Total 184,0957 100,000
De los resultados obtenidos se puede decir que el primer piso moviliza el 83%
aproximadamente de la masa total, mientras que el segundo piso moviliza el 11%, el tercer
piso moviliza el 4% y el ultimo nivel solo moviliza el 1% de la masa total. Se tiene entonces
que el modo 1 es el modo dominante en la estructura.
j. Verificar los resultados obtenidos (periodos de vibración, modos de vibración) con

los resultados de un programa de análisis (SAP2000, ETABS, u otro).
Sap Periodo obtenido

Modo Periodo[s] Modo Periodo[s]
17
1 0,47056 1 0,04718
2 0,1468 2 0,1462
3 0,08138 3 0,0813
4 0,05748 4 0,05746
Los periodos con el método 2 corresponden a los periodos obtenidos en sap de cada modo de
vibración.
3. RESPUESTA SÍSMICA – PARTE 1
Usando el método del análisis modal espectral, y el espectro de diseño de la

Microzonificación Sísmica de Bogotá Lacustre 200, calcular lo siguiente:
1. Fuerzas sísmicas elásticas en cada nivel de la edificación.
Para determinar fuerzas sísmicas es necesario calcular el valor de la aceleración espectral

correspondiente para cada uno de los modos y por medio de las relaciones espectrales se
calcula el desplazamiento espectral, para esto se tendrá en cuenta el valor del periodo y el
espectro de aceleración como se muestra a continuación:
ESPECTRO ELASTICO DE ACELERACIONES

ZONA BOGOTA Lacustre 200
TIPO DE SUELO 0 Media
Aceleración pico efectiva-Aa 0,15
Tabla A.2.3-2
Velocidad pico efectiva - Av 0,20
Coeficiente Fa 1,2 Tabla A.2.4-3
Coeficiente Fv 3,5 Tabla A.2.4-4
Coeficiente i 1,0 Tabla A.2.5-1
To 0,38888889
Tc 1,86666667
TL 4
T(s) Sa(g)
0 0,18
0,003 0,182082857
0,006 0,184165714
0,051 0,215408571
0,054 0,217491429
0,057 0,219574286
18
0,06 0,221657143
0,063 0,22374
0,075 0,232071429
0,078 0,234154286
0,081 0,236237143
0,084 0,23832
0,087 0,240402857
0,09 0,242485714
0,138 0,275811429
0,141 0,277894286
0,144 0,279977143
0,147 0,28206
0,15 0,284142857
0,375 0,440357143
0,378 0,44244
0,381 0,444522857
0,384 0,446605714
0,387 0,448688571
0,39 0,45
1,854 0,45
1,857 0,45
1,86 0,45
1,863 0,45
1,866 0,45
3,975 0,211320755
3,978 0,211161388
3,981 0,211002261
3,984 0,210843373
3,987 0,210684725
3,99 0,210526316
3,993 0,210368144
3,996 0,21021021
3,999 0,210052513
4,002 0,209790157
19
A continuación, se muestra en conjunto el periodo, el coeficiente de participación modal, la

frecuencia angular al cuadrado, la aceleración espectral y la relación espectral para cada piso
de la estructura. Las relaciones espectrales se calculan con la siguiente ecuación:
𝑆𝑎
𝑆𝑑 =
𝑤2
MODO T α w^2 Sa(Tl) Sa(m/s^2) Sd(Tl)

1 0,472 12,397 177,323 0,45 4,41297 0,02488661
2 0,146 4,516 1844,748 0,28 2,745848 0,00148847
3 0,081 2,767 5971,863 0,24 2,353584 0,00039411
4 0,057 1,534 11953,867 0,22 2,157452 0,00018048
Donde 1g = 9,806 m/s^2
Para calcular las fuerzas sísmicas por piso se emplea la siguiente ecuación:
➢ Modo 1
20
𝐹11 48,8 0 0 0 0,02473483

𝐹21 0 48,8 0 0 0,06046684
1 =[ 0 0 48,8 0
]𝑥 [
0,08876019
] 𝑥 12,397 𝑥 4,4145
𝐹3
1 0 0 0 37,5 0,10421191
{𝐹4 }
𝐹11 66,05974259
𝐹21 161,4898575
= [ ]
𝐹31 237,0534019
1 213,873383
{𝐹4 }
➢ Modo 2
𝐹12 48,8 0 0 0 0,06907433

𝐹22 0 48,8 0 0 0,09169138
2 =[ 0 0 48,8 0
]𝑥 [
0,00652167
] 𝑥 4,516 𝑥 2,7468
𝐹3
2 0 0 0 37,5 −0,09727077
{𝐹4 }
𝐹12 41,81331645
𝐹22 55,50427743
2 = [ 3,947813107 ]
𝐹3
2 −45,24720204
{𝐹4 }
➢ Modo 3
𝐹13 48,8 0 0 0 0,09539934

𝐹23 0 48,8 0 0 −0,00739499
3 =[ 0 0 48,8 0
]𝑥 [
−0,08615033
] 𝑥 2,767 𝑥 2,3544
𝐹3
3 0 0 0 37,5 0,07137003
{𝐹4 }
𝐹13 30,32713825
𝐹23 −2,350844428
= [ ]
𝐹33 −27,38690829
3 17,43465619
{𝐹4 }
➢ Modo 4
21
𝐹14 48,8 0 0 0 0,077603

𝐹24 0 48,8 0 0 −0,091643
=[ ]𝑥[ ] 𝑥[1,534]𝑥[2,1582]
𝐹34 0 0 48,8 0 0,0718663
4 0 0 0 37,5 −0,034340
{𝐹4 }
F14 12,538748
F24 −14,807374
=[ ]
F34 11,611827
4 −4,2637314
{ F4 }
Para obtener el valor de las fuerzas sísmicas totales para cada piso se elevan al cuadrado los
vectores obtenidos anteriormente, se suman y se aplica raíz cuadrada y se obtiene que:
𝐹1 66,0597 2 41,8133 2 30,3271 2 12,5388 2

𝐹 161,4899 55,5043 −2,3508 −14,8074
{ 2 } = √{ } +{ } +{ } +{ }
𝐹3 237,0534 3,9478 −27,3869 11,6118
𝐹4 213,8733 −45,2472 17,4347 −4,2637
𝐹1 84,7891
𝐹2 171,4190
{ }={ } 𝑘𝑁
𝐹3 238,9451
𝐹4 219,3428
La fuerza sísmica elástica para

✓ Piso 1 será de aproximadamente 85 kN
✓ Piso 2 será aproximadamente de 172 kN
2. Fuerzas sísmicas de diseño en cada nivel de la edificación.
Para determinar las fuerzas sísmicas de diseño para cada nivel se requiere dividir los valores
obtenidos del vector de fuerza elástica entre un coeficiente de disipación de energía básico
cuyo valor es 5 y se obtuvo de la NSR-10. La ecuación con la que se calcula las fuerzas
sísmicas de diseño es la siguiente:
𝐹𝑠
𝐸=
𝑅0
22
Se toma el vector de fuerzas obtenidas en el numeral 1 y se divide cada valor entre el valor
del coeficiente de disipación de energía básico ya nombrado.
84,7891
171,4190
{ }
𝐹1 238,9451
𝐹 219,3428
{ 2} = 𝑘𝑁
𝐹3 5
𝐹4
Vector de fuerzas sísmicas de diseño:
𝐹1 16,9578
𝐹 34,2838
{ 2} = { } 𝑘𝑁
𝐹3 47,7890
𝐹4 43,8586
3. Cortante basal total elástico y de diseño.
El cortante basal de la estructura corresponde a la reacción que se presenta en la superficie

debido a la aplicación de las fuerzas sísmicas en cada uno de los pisos de la estructura, por
lo tanto, el cortante basal de cada modo de vibración serán las fuerzas aplicadas, se tiene que
el cortante basal para las fuerzas elásticas y de diseño serán:
23
➢ Cortante basal – Fuerzas elástico
➢ Modo 1
𝐹11 66,05974259
𝐹21 161,4898575
= [ ] = 66,0597 + 161,4899 + 237,0534 + 213,8734 = 678,4764 𝑘𝑁
𝐹31 237,0534019
1 213,873383
{𝐹4 }
➢ Modo 2
𝐹12 41,81331645
2
𝐹2 55,50427743
2 = [ 3,947813107 ] = 41,8133 + 55,5043 + 3,9478 − 45,2472 = 56,0182 𝑘𝑁
𝐹3
2 −45,24720204
{𝐹4 }
➢ Modo 3
𝐹13 30,32713825
𝐹23 −2,350844428
= [ ] = 30,3271 − 2,35084 − 27,3869 + 17,4347 = 18,02404 𝑘𝑁
𝐹33 −27,38690829
3 17,43465619
{𝐹4 }
➢ Modo 4
F14 12,538748
4
F2 −14,807374
4 = [ 11,611827 ] = 12,5388 − 14,8074 + 11,6118 − 4,2637 = 5,07947 𝑘𝑁
F3
4 −4,2637314
{ F4 }
Para hallar el cortante basal teniendo en cuenta los 4 modos de vibración correspondientes
se procede a hallar la magnitud total del cortante basal de las fuerzas elásticas. Por lo tanto,
se elevan al cuadrado los resultados obtenidos, se suman y se aplica raíz cuadrada al
resultado. Se obtiene entonces que el cortante basal debido a las fuerzas elásticas es:
𝑉𝑏 = √(678,4764)2 + (56,0182)2 + (18,02404)2 + (5,07947)2
𝑉𝑏 = 681,0425 𝑘𝑁
24
➢ Cortante basal – Fuerzas de diseño
Se tiene que:
➢ Modo 1
𝐹11 13,211949
𝐹21 32,297972
= [ ] = 13,211949 + 32,29797 + 47,4107 + 42,7747 = 135,69532 𝑘𝑁
𝐹31 47,410680
1 42,774677
{𝐹4 }
➢ Modo 2
𝐹12 8,362663289
2
𝐹2 11,10085549
2 = [ 0,789562621 ] = 8,3627 + 11,10086 + 0,78956 − 9,04944 = 11,20368 𝑘𝑁
𝐹3
2 −9,049440409
{𝐹4 }
➢ Modo 3
𝐹13 6,065427649
𝐹23 −0,470168886
= [ ] = 6,06543 − 0,470169 − 5,47739 + 3,4869 = 3,604771 𝑘𝑁
𝐹33 −5,477381658
3 3,486931238
{𝐹4 }
➢ Modo 4
F14 2,5077496
4
F2 −2,9614749
4 = [ 2,3223655 ] = 2,50775 − 2,961475 + 2,32237 − 0,852746 = 1,01589 𝑘𝑁
F3
4 −0,8527463
{ F4 }
Para obtener el cortante basal total se realiza el siguiente procedimiento con los valores de
cada modo obtenidos anteriormente:
𝑉𝑏 = √(135,69532)2 + (11,20368)2 + (3,604771)2 + (1,01589)2
𝑉𝑏 = 136,208549 𝑘𝑁
25
4. Desplazamiento lateral de cada piso y graficar el perfil de desplazamientos.
Para obtener los desplazamientos de cada uno de los pisos se emplea la siguiente ecuación:
Es necesario tener cada uno de los modos normalizados, el coeficiente de participación y el

desplazamiento obtenido por medio de las relaciones espectrales cuyos valores ya fueron
calculados en el numeral 1 del punto 3.
➢ Modo 1
𝑥11 0,02473483 0,007631349

𝑥21 0,06046684 0,01865562
1 = {0,08876019} 𝑥 12,3973 𝑥 0,024887 = {0,027384866}
𝑥3
1 0,10421191 0,032152128
{𝑥4 }
➢ Modo 2
𝑥12 0,06907433 0,000464309

𝑥22 0,09169138 0,000616338
2 = { 0,00652167 } 𝑥 4,51597 𝑥 0,001488 = {4,38379𝐸 − 05}
𝑥3
2 −0,09727077 −0,000653842
{𝑥4 }
➢ Modo 3
𝑥13 0,09539934 0,000104028

𝑥23 −0,00739499 −8,06387𝐸 − 06
3 = {−0,08615033} 𝑥 2,76685 𝑥 0,000394 = {−9,39427𝐸 − 05}
𝑥3
3 0,07137003 7,78255𝐸 − 05
{𝑥4 }
➢ Modo 4
26
𝑥14 0,07760315 2,1487𝐸 − 05

𝑥24 −0,09164383 −2,53746𝐸 − 05
= { } 𝑥 1,5341 𝑥 0,000180 = { }
𝑥34 0,07186637 1,98986𝐸 − 05
4 −0,03434026 −9,50823𝐸 − 06
{𝑥4 }
Para hallar el desplazamiento total se realiza el siguiente procedimiento teniendo en cuenta

los vectores para cada piso obtenidos anteriormente.
2
𝑥1 0,0076314 2 0,00046431 2 0,000104028 2 2,1487 ∗ 10−5
𝑥2 0,0186556 0,00061634 −8,06387 ∗ 10 −6
−2,53746 ∗ 10 −5
{𝑥 } = √{ } +{ } +{ } +{ }
3 0,0273849 4,38379 ∗ 10−5 −9,39427 ∗ 10−5 1,98986 ∗ 10−5
𝑥4 0,0321521 −0,000653842 7,78255 ∗ 10 −5 −9,50823 ∗ 10 −6
𝑥1 0,007646199
𝑥2 0,018665817
{𝑥 } = { } 𝑚
3 0,02738507
𝑥4 0,032158871
Graficando el respectivo grafico de perfil de desplazamiento se tiene que en el eje X el

desplazamiento resultante de cada piso y en el eje Y la elevación del piso.
Perfil de desplazamientos
14
12
10
Elevacion Piso
0
0 0,005 0,01 0,015 0,02 0,025 0,03 0,035
Desplazamiento resultante
5. Derivas por piso.
27
Las derivas son el desplazamiento relativo entre los pisos de la estructura. A continuación,
se realizan las respectivas diferencias entre desplazamientos para cada uno de los modos de
vibración:
➢ Modo 1
∆11 0,007631349 − 0 0,00763135

1
∆2 0,01865562 − 0,0076314 0,01102427
1 = { 0,027384866 − 0,01865562 } = {0,00872925} 𝑚
∆3
1 0,032152128 − 0,027384866 0,00476726
{∆4 }
➢ Modo 2
∆12 0,000464309 − 0 0,000464309

∆22 0,000616338 − 0,000464309 0,0000152029
={ }={ }𝑚
∆23 4,38379 ∗ 10−5 − 0,000616338 −0,0005725
2
{∆4 } −0,000653842 − 4,38379 ∗ 10−5 −0,00069768
➢ Modo 3
∆13 0,00010403 − 0 0,000104028

3 −6
∆2 −8,06387 ∗ 10 − 0,00010403 −0,000112092
3 = −5 −6 ={ }𝑚
∆3 −9,39427 ∗ 10 − (−8,06387 ∗ 10 ) −8,58788 ∗ 10−5
3
{∆4 } { 7,78255 ∗ 10−5 − (−9,39427 ∗ 10−5 ) } 0,000171768
➢ Modo 4
∆14 2,1487𝐸 − 05 − 0 2,1487 ∗ 10−5

4 −5 −5
∆2 −2,53746 ∗ 10 − 2,1487 ∗ 10 −4,68616 ∗ 10−5
4 = −5 −5 = 𝑚
∆3 1,98986 ∗ 10 − (−2,53746 ∗ 10 ) 4,52732 ∗ 10−5
4
{∆4 } { −9,50823 ∗ 10−6 + 1,98986 ∗ 10−5 } {−2,94068 ∗ 10−5 }
Para hallar las derivas totales por piso se realiza el siguiente procedimiento teniendo en
cuenta los vectores de derivas para cada piso obtenidos anteriormente.
2
∆1 0,00763135 2 0,000464309 2 0,000104028 2 2,1487 ∗ 10−5
∆ 0,01102427 0,0000152029 −0,000112092 −4,68616 ∗ 10−5
{ 2 } = √{ } +{ } +{ −5 } +{ }
∆3 0,00872925 −0,0005725 −8,58788 ∗ 10 4,52732 ∗ 10−5
∆4 0,00476726 −0,00069768 0,000171768 −2,94068 ∗ 10−5
28
∆1 0,007646199
∆ 0,011025989
{ 2} = { }𝑚
∆3 0,008748538
∆4 0,004821194
6. Verificar el cumplimiento de derivas de acuerdo con NSR-10.
Según la norma NSR-10, para una estructura de concreto reforzado, el valor de la deriva
correspondiente de cada uno de los pisos no debe exceder más del 1% de la altura de piso.
En la siguiente tabla se muestra el cumplimiento de las derivas obtenidas en el numeral 5
según como lo dice la norma.
Tomado: NSR-10.
Deriva
Altura Deriva Criterio de
Piso máxima
(m) obtenida (m) cumplimiento
permitida (m)
1 3 0,03 0,007646199 CUMPLE NSR-10
2 3 0,03 0,011025989 CUMPLE NSR-10
3 3 0,03 0,008748538 CUMPLE NSR-10
4 3 0,03 0,004821194 CUMPLE NSR-10
Como se puede observar en la tabla anterior las derivas obtenidas no exceden los valores
estipulados por la norma para este caso, por lo tanto, las derivas de la estructura si cumplen
con la normativa.
4. RESPUESTA SÍSMICA – PARTE 2
Usando el método de la fuerza horizontal equivalente, y el espectro de diseño de la

Microzonificación Sísmica de Bogotá asignado a cada grupo, calcular lo siguiente:
29
1. Fuerzas sísmicas elásticas en cada nivel de la edificación

2. Fuerzas sísmicas de diseño en cada nivel de la edificación
3. Cortante basal total elástico y de diseño
4. Desplazamiento lateral de cada piso y graficar el perfil de desplazamientos
5. Derivas por piso
6. Verificar el cumplimiento de derivas de acuerdo con NSR-10
Para el desarrollo de esta parte se tiene la siguiente tabla en la cual se tienen los datos de
entrada.
PARAMETRO SÍMBOLO REFERENCIA VALOR UNIDADES

Velocidad Pico Av A.2.3-2 0,2 [-]
Coeficiente de amp. Fv A.2.4-4 3,5 [-]
Altura Total Edif. h Guía Taller 12 [m]
Masa Total Mtot Guía Taller 183,900 [ton]
Ct A.4.2.1 0,047 [-]
α A.4.2.1 0,9 [-]
Periodo Aproximado Ta 0,44 [s]
Periodo Exacto T SAP2000 0,471 [s]
Cu 0,91 [-]
Periodo máximo Tmax 0,400 [s]
Coeficiente Básico Ro A.3.3 5 [-]
Aceleración espectral Sa(T) Esp. aceleraciones 0,45 [g]

Cortante basal total
Vb Mtot*Sa(T)*g 811,82655 [kN]
elástico
Cortante basal total diseño Vb/Ro 162,36531 [kN]
Coeficiente k k T < 0,5 1 [-]
A partir del espectro de aceleraciones se obtiene el valor de la aceleración espectral (Sa(T)),

lo cual permite calcular el cortante basal total con la formula presentada anteriormente, que
da como resultado un cortante de 811,827 kN.
Teniendo claro con el periodo que se va a trabajar se usa el coeficiente k=1 que cumple con
el requisito presentado en la tabla anterior.
𝑚𝑗 ∗ ℎ𝑗𝑘
𝐹𝑗 = 𝑝 ∗ 𝑉𝑏
∑𝑗=1(𝑚𝑗 ∗ ℎ𝑗𝑘 )
30
Se usa la fórmula presentada anteriormente para obtener el valor de fuerza sísmica elástica
(Fj), en la se cual utiliza los valores de las masas por piso, las alturas y el cortante basal total.
En la siguiente tabla se presentan los valores obtenidos para las fuerzas sísmicas elásticas por
piso, junto con sus respectivas fuerzas sísmicas de diseño obtenidas a partir la razón de
cambio entre la fuerza elástica y el coeficiente de disipación de energía básico, Ro.
Elástica Diseño
PISO mj[ton] hj[m] mj hjk fj [kN] E [kN]
1 48,8 3 146,4 89,4695927 17,89391854
2 48,8 6 292,8 178,939185 35,78783707
3 48,8 9 439,2 268,408778 53,68175561
4 37,5 12 450 275,008994 55,00179878
TOTAL 1328,4 811,82655
Con las fuerzas sísmicas calculadas anteriormente, se procede a asignarlas en el modelo de

SAP2000 para obtener el valor del desplazamiento de cada piso (𝛿x).
En la siguiente tabla se presentan los valores de desplazamientos laterales obtenidos junto

con el cálculo de las derivas por piso, las cuales están definidas por el desplazamiento de un
piso respecto al inmediatamente anterior.
REVISIÓN DERIVAS
PISO hk [m] 𝛿x [m] Δx [m] Δlim [m] Verificar NSR-10
Base - 0 - - -
1 3 0,009 0,009 0,03 Cumple
2 3 0,022 0,013 0,03 Cumple
3 3 0,0325 0,0105 0,03 Cumple
4 3 0,0386 0,0061 0,03 Cumple
Después de obtenidas las derivas por piso, se procede a realizar la verificación de

cumplimiento según el criterio de la NSR-10 presentado en la siguiente tabla.
31
Tomado: NSR-10.
Con lo anterior, el criterio establece que estas derivas no pueden superar más del 1% de la
altura de piso, que para este caso el valor límite es de 0,03 m.
Al comparar los valores obtenidos con el valor límite, se tiene que las derivas de todos los
pisos cumplen con el criterio establecido por la norma NSR-10.
PERFIL DE DESPLAZAMIENTO
14
12
10
8
hj [m]
0
0 0,005 0,01 0,015 0,02 0,025 0,03 0,035 0,04 0,045
𝛿x [m]
En la anterior gráfica se observa el perfil de desplazamiento obtenido al graficar los valores

de las alturas de cada piso contra el valor obtenido de los desplazamientos.
32

Anae Grupo3 TALLER3

Cargado por

Copyright:

Formatos disponibles

Anae Grupo3 TALLER3

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Anae Grupo3 TALLER3

Cargado por

Copyright:

Formatos disponibles

ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO

FACULTAD DE INGENIERÍA CIVIL

ING. HECTOR JOSÉ PEREZ BARRERA

ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO

a. Matriz de masa del sistema

b. Matriz de rigidez lateral usando el método 2 visto en clase.

Del programa se obtuvo que el desplazamiento para el piso 1 es de 0.0003106 m.

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = 𝐹11 → ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = 100𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎2 = 𝐹21 + 𝐹22 + 𝐹32 + 𝐹42

→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 = −14.86𝑘𝑁 − 14.2𝑘𝑁 − 14.42𝑘𝑁 − 13.68𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 = 𝐹31 + 𝐹32 + 𝐹33 + 𝐹34

→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 = 5.52𝑘𝑁 + 2.68𝑘𝑁 + 3.3𝑘𝑁 + 2.42𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = 𝐹41 + 𝐹42 + 𝐹43 + 𝐹44

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 100

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 −57.16

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 13.92

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 −1.7815

Del programa se obtuvo que el desplazamiento para el piso 2 es de 0.0003594 m.

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = 𝐹11 + 𝐹21 + 𝐹31 + 𝐹41

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 = 𝐹22 → ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 = 100𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 = 𝐹31 + 𝐹32 + 𝐹33 + 𝐹34

→ ∑ 𝑅𝑒𝑎𝑐𝑐𝑖𝑜𝑛𝑒𝑠3 = −13.66𝑘𝑁 − 16.31𝑘𝑁 − 15.67𝑘𝑁 − 16.73𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎𝑠 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = 𝐹41 + 𝐹42 + 𝐹43 + 𝐹44

→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = −0.72𝑘𝑁 + 3.33𝑘𝑁 + 2.48𝑘𝑁 + 7.51𝑘𝑁

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 −66.4

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 −100

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 −62.37

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 12.6

Del programa se obtuvo que el desplazamiento para el piso 3 es de 0.0003906 m.

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = 𝐹11 + 𝐹21 + 𝐹31 + 𝐹41

→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = 6.43𝑘𝑁 + 3.66𝑘𝑁 + 3.87𝑘𝑁 + 3.55𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 = 𝐹21 + 𝐹22 + 𝐹32 + 𝐹42

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 = 𝐹33 → ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 = 100𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = 𝐹41 + 𝐹42 + 𝐹43 + 𝐹44

→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = −14.52𝑘𝑁 − 12.68𝑘𝑁 − 13.51𝑘𝑁 − 6.26𝑘𝑁

ΣReacciones4 – R4 (kN) -46.97

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 17.51

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 −67.8

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 100

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 −46.97

Por último, en el Piso 4 en el programa se obtuvo que el desplazamiento es de 0.00111 m.

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = 𝐹11 + 𝐹21 + 𝐹31 + 𝐹41

→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 = −0.56𝑘𝑁 − 1.3𝑘𝑁 − 1.1𝑘𝑁 − 3.44𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 = 𝐹21 + 𝐹22 + 𝐹32 + 𝐹42

→ ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 = 9.91𝑘𝑁 + 9.1𝑘𝑁 + 8.87𝑘𝑁 + 11.08𝑘𝑁

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 = 𝐹31 + 𝐹32 + 𝐹33 + 𝐹34

∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = 𝐹44 → ∑ 𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎4 = 100𝑘𝑁

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎1 −6.4

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎2 38.96

𝐶𝑎𝑟𝑔𝑎 𝑎𝑝𝑙𝑖𝑐𝑎𝑑𝑎3 −133.59

𝑘11 𝑘12 𝑘13 𝑘14

321957.5 −184028.937 44828.469 −5765.7658

c. Matriz de rigidez lateral usando el método 3 (edificio cortante) visto en clase.

Teniendo en cuenta lo siguiente

𝐹11 𝐹12 𝐹13 𝐹14

Se obtiene la matriz de rigidez

d. Periodos de vibración de todos los modos de vibración de la estructura, usando

Mediante Microsoft Excel se obtiene las frecuencias de cada modo de vibración

e. Periodos de vibración de todos los modos de vibración de la estructura, usando