Tarea Tema 3

Taller de Dinámica Estructural e Ingeniería

Sísmica OAL- 31-08-2021
Respuesta Sísmica de un Pórtico de Tres
Pisos

La estructura indicada en la figura es un pórtico de concreto armado con vigas

idealizadas como infinitamente rígidas, columnas indeformables axialmente y masas
m1, m2 y m3 concentradas en cada nivel. Posee tres grados dinámicos de libertad
asociados a la traslación horizontal de cada piso. Las rigideces de cada entrepiso son K1,
K2 y K3. La altura de cada entrepiso es de 3 metros y la longitud entre columnas es de 6
metros. Use un concreto de 250 kgf/cm2 y su módulo de elasticidad correspondiente
(15100 f’c0,5).

Se trata de determinar la respuesta estructural (fuerzas, desplazamientos) ante el sismo

de diseño establecido en la Norma COVENIN 1756. El edificio está ubicado en Caracas,
sobre un terreno firme calificado como Clase BC duro con una profundidad a la roca de
30 metros y una condición topográfica Leve. Pertenece al Grupo B2. A la estructura de
concreto armado se le asigna un factor de reducción de R=6. Se puede adoptar un Factor
de Irregularidad FI=1 y un Factor de Redundancia =1.

Se pide lo siguiente:

1- Determine la matriz de rigidez y la matriz de masa de la estructura. Determine los

desplazamientos laterales para una carga lateral de 180 t en el nivel 3. Determine
los modos y los periodos de vibración. Normalice los modos tal que el
desplazamiento del último piso sea igual a 1. Dibuje la deformada de la estructura
en cada modo de vibración. Determine los factores de participación y las masas
participativas de cada modo para sismo horizontal.

2- Determine mediante el método de análisis espectral y utilizando cálculo manual

los siguientes valores de la respuesta máxima de la estructura ante el sismo
especificado por la norma 1756: a) Cortante en cada entrepiso; b) Desplazamiento
en cada nivel, elástico y total; c) Deriva lateral elástica y total en cada entrepiso;
d) Momento flector en la base de la columna derecha; d) Fuerza cortante en la
 base de la columna derecha; e) Fuerza axial en la columna derecha del tercer
entrepiso. Muestre los resultados parciales: aceleraciones espectrales, respuestas
en cada modo de vibración y coeficientes de correlación.

3- Selecciones las dimensiones de las columnas a fin de tener las rigideces del
pórtico indicadas en el enunciado del problema. Determine mediante el método
de análisis espectral y utilizando un programa computacional de su preferencia,
la respuesta máxima de la estructura ante el sismo indicado. a) Determine los
desplazamientos laterales para una carga estática de 180 t en el nivel 3, y compare
con el resultado del cálculo manual de la pregunta 1; b) Determine los períodos y
modos de vibración y compare con el resultado del cálculo manual de la pregunta
1; c) Determine los mismos valores de respuesta solicitados en la pregunta 2 y
compare con los resultados del cálculo manual. Comente.
1) Determine la matriz de rigidez y la matriz de masa de la estructura.

a) Determine los desplazamientos laterales para una carga lateral de 180 t en

el nivel 3.
b) Determine los modos y los periodos de vibración. Normalice los modos tal
que el desplazamiento del último piso sea igual a 1.
c) Dibuje la deformada de la estructura en cada modo de vibración.
d) Determine los factores de participación y las masas participativas de cada
modo para sismo horizontal.

Datos:

H= 300 centímetros (altura de cada entrepiso)

L= 600 centímetros (longitud entre columnas)
F’c= 250 kgf/cm²
E=15100 f’c0,5 = 238.751,96 kg/cm2 = 238,75 t/cm²

Determinando la matriz de masa:

m= t/cm/seg²

Determinando la matriz de rigidez a través de los desplazamientos por nivel

Quedando la matriz de rigidez:

a) Determine los desplazamientos laterales para una carga lateral de 180 t en el

nivel 3.

k=

A través del método de igualación hallaré el valor de los desplazamientos u1, u2 y u3:
Resultado los desplazamientos por nivel para la carga de 180 ton en el tercer nivel:

U1= 1 cm
U2= 2,5 cm
U3= 5,5 cm

Deriva piso 1= 1 cm – 0 cm= 1 cm

Deriva piso 2= 2,5 cm – 1 cm = 1,5 cm
Deriva piso 3= 5,5 cm – 2,5 cm = 3 cm

b) Determine los modos y los periodos de vibración. Normalice los modos tal que
el desplazamiento del último piso sea igual a 1.

Multiplicando la diagonal principal por el valor de landa y por la matriz de masa

obtengo los resultados para la ecuación cubica:
Ordenando los autovalores de menor a mayor:

Frecuencia 1= √188,3156 seg = 13,7228 seg

Frecuencia 2= √900 seg = 30 seg
Frecuencia 3= √1911,68 seg = 43,7228 seg

T1= (2 x pi)/13,7228 seg

T1= 0,4578 seg

T2= (2 x pi)/30 seg

T2= 0,209 seg

T3= (2 x pi)/43,7228 seg

T3= 0,1437 seg

Calculando el primer modo: para landa 1= 188,3156 seg

Multiplicando la matriz de rigidez por el escalar de landa 1=

Quedando la resta de la matriz de rigidez menos la matriz de masa multiplicada por el

landa1:

Matriz “A”= (matriz de rigidez) menos matriz “B” (matriz de masa multiplicada por
188,3156) y matriz “C” (matriz resultante de restar “A” – “B”)
A través del método de reducción hayo los valores para la matriz “C”

Φx1= 0,313863 cm
Φy1= 0,68615 cm
Φ Φ1= 1 cm

Calculando los valores para el segundo modo: landa 2= 900 seg

Multiplicando la matriz de rigidez por el escalar de landa 2=

Quedando la resta de la matriz de rigidez menos la matriz de masa multiplicada por el

landa 2:

Matriz A menos matriz B resultando la matriz:

Φx2= - 0,5 cm
Φy2= - 0,5 cm
Φ Φ2= 1 cm
Calculando el tercer modo: landa 3= 1911,684 seg

Restando la matriz A menos la matriz B

A través del método de igualación obtengo los valores para el modo 3:

Φx3= 3,1850 cm
Φy3= - 2,18613 cm
Φ Φ3= 1 cm

c) Dibuje la deformada de la estructura en cada modo de vibración.

 d) Determine los factores de participación y las masas participativas de cada
modo para sismo horizontal.

Modo 1:

X =

x =

x =
 x =

Participación del modo 1 = 0,299999/ 0,2138544525 = 1,402823254

Modo 2:

x =

x =

x =
Participación del modo 2 = -0,1 / 0,2 = -0,5

Modo 3:

x =

x =

x =

Participación del modo 3 = 0,299774 / 3,18467787538 = 0,094130085

T1= 0,2999999 / 0,2138544525 = 1,402823212
T2= 0,1 / 0,2 = 0,5
T3= 0,299774 / 3,18467787538 = 0,094130085

Cálculo de las masas participativas:

Bx1= 0,299999 ² / 0,2138544525 = 0,420846977 t/seg²/cm

Bx2= -0,1 ² / 0,2 = 0,05 t/seg²/cm
Bx3= = 0,299774 ² / 3,18467787538 = 0,028217752 t/seg²/cm

Bx total= 0,499067429 = 0,50

Porcentajes de participación:

Modo 1: 0,420849699 / 0,50 = 84,1699398 %

Modo 2: 0,05 / 0,50 = 10 %
Modo 3: 0,028217752 / 0,50 = 5,6435504 %

2) Determine mediante el método de análisis espectral y utilizando cálculo manual

los siguientes valores de la respuesta máxima de la estructura ante el sismo
especificado por la norma 1756:

a) Cortante en cada entrepiso;

b) Desplazamiento en cada nivel, elástico y total;
c) Deriva lateral elástica y total en cada entrepiso;
d) Momento flector en la base de la columna derecha;
d) Fuerza cortante en la base de la columna derecha;
e) Fuerza axial en la columna derecha del tercer entrepiso. Muestre los
resultados parciales: aceleraciones espectrales, respuestas en cada modo de
vibración y coeficientes de correlación.

Valor fijado para A0= 0,30

Valor fijado para A1= 0,24

Valor fijado para Tl= 3,9

Factor de importancia B2: 1

 Clase de sitio BC

Profundidad al basamento rocoso= 30 m

Condición Topográfica: Leve
Datos a considerar para calcular el espectro elástico

A0 (aceleración del terreno, periodo 0 seg) = 0,30 seg

A1 (periodo para estructuras de periodo 1 seg= 0,24 seg
Tl (periodo de transición para periodos largos) = 3,9 seg
Clase BC Vs= 850 m/s
El edificio posee una topografía leve o sea < a 40% y para una profundidad al basamento
rocoso de 30 metros Fat= 1,00 ; Fvt=1,00; Fdt= 1,00
Factor de importancia= B2= 2,2
aA0= 1 x 0,30 = 0,30
aA1= 1 x 0,24 = 0,24
Los valores obtenidos para la clase de sitio son: Fac= 1,00 ; Fcv= 1,00 ; Fdc= 1,00
hallando Fca
Obteniendo los valores para profundidad del basamento rocoso de 30 m Fah= 1,00 ;
Fvh= 1,00 ; Fdh= 1,00

Obteniendo los valores para profundidad del basamento rocoso de 30 m Fah= 1,00 ;
Fvh= 1,00 ; Fdh= 1,00

Calculando Fa
Fa= Fac x Fah x Fat
Fa= 1,00 x 1,00 x 1,00 = 1,00

Calculando Fv
Fv= Fvc x Fvh x Fvt
Fv= 1,00 x 1,00 x 1,00 = 1,00

Calculando Fd
Fd= Fdc x Fdh x Fdt
Fd= 1,00 x 1,00 x 1,00 = 1,00

Calculando el parámetro Aa (aceleración para periodo 0 o aceleración del terreno)

Aa= fa x a x A0
Aa= 1,00 x 1 x 0,30
Aa= 0,30 seg

Calculando el parámetro Av (aceleración para periodo 0 o aceleración del terreno)

Av= fv x a x A1
Av= 1,00 x 1 x 0,24
Av= 0,24 seg

Siendo en amortiguamiento 5% = 0,05

El parámetro “B” (beta) o factor de amplificación para amortiguamiento de 5% es de 2,4
“B*” será el mayor entre 2,4 y (Av/Aa= 0,24/0,30= 0,8), tomaré el valor de 2,4
Calculando los periodos para hacer el espectro

Tc= 1/B x Av/Aa

Tc= ( 1 / 2,4) x (0,24 / 0,30)
Tc= 0,416 x 0,8
Tc= 0,333333333

Tb= 0,25 x Tc
Tb= 0,25 x 0,33333333
Tb= 0,083333333

Ta= 0,20 x Tb
Ta= 0,20 x 0,33333333
Ta= 0,01666666666

Td= Tl x (Fd/Fv)
Td= 3,9 x (1,00 / 1,00)
Td= 3,9 x 1,00
Td= 3,9 seg

Factor de irregularidad = 1 (es regular)

Factor de redundancia = 1
Factor de reducción = 6

Hallando el valor de T+ para el puente siendo R= 3, obtengo

T+= para R > 5, seria 0,4

T+= 0,4 seg
Acotado por la expresión:

Entre (0,0833333 y 0,3333333)

Asumiré T+ = Tc = 0,3333333 seg

Hallando el valor para periodo 0

Ad (T)= (1 x 1 x 0,30) / 1,5

Ad (T)= 0,20 seg
Hallando el espectro de aceleración inelástico
Ad (T)= ((espectro de aceleración elástico x ro x F1) / R)
Ad (T)= (0,72 x 1 x 1) / 6
Ad (T)= 0,12

Hallando los valores de ad1, ad2 y ad3 a través de periodos obtenidos en los 3 modos
para el espectro inelástico:

Modo 1: 0,4578 seg (entre t+ y td)

Ad1= 1 x 1 x ((2,40 x 0,30) / 6) x (0,333333 / 0,4578)

Ad1= 0,12 x 0,7281
Ad1= 0,087372

Modo 2: 0,209 seg (entre tb y t+)

Ad2= 1 x 1 x 0,30 ( 1 / 1,5 + ( 2,40 / 6 – 1 / 1,5) (0,209 – 0,02 / 0,333 – 0,02)
Ad2= 0,30 x ( 0,6666 + (- 0,26666 x 0,603833865 )
Ad2= 0,151692498

Modo 3: 0,1437 seg (t < ta)

Ad3= 1 x 1 x 0,30 ( 1 / 1,5 + ( 2,40 / 6 – 1 / 1,5) (0,1437 – 0,02 / 0,333 – 0,02)

Ad3= 0,30 x ( 0,6666 + (- 0,26666 x 0,395207667 )
Ad3= 0,168383392

Resultando:

Ad1= 0,087372
Ad2= 0,151692498
Ad3= 0,168383392

Cálculo de las aceleraciones:

A1= 0,087372 x 981 = 85,711932 cm/s²

A2= 0,151692498 x 981 = 148,810341 cm/s²
A3= 0,168383392 x 981 = 165,1841076 cm/s²

Cálculo de los desplazamientos espectrales:

Frecuencia circular 1= √188,3156 seg = 13,7228 seg

Frecuencia circular 2= √900 seg = 30 seg
Frecuencia circular 3= √1911,68 seg = 43,7228 seg

D1 = 85,711932 cm/s² / 188,3152398 s² = 0,45515133 cm

D2 = 148,810341 cm/s² / 900 s² = 0,165344823 cm
D3 = 165,1841076 cm/s² / 1.911,68324 s² = 0,086407676 cm

Cálculo de los vectores de desplazamiento para los 3 modos:

Modo 1:
Modo 2:

Modo 3:

El vector de desplazamiento para el modo 1 queda:

( 0,200400523 ; 0,438104585 ; 0,638496808 )

El vector de desplazamiento para el modo 2 queda:

( 0,041336205 ; 0,041336205 ; - 0,082672411 )

El vector de desplazamiento para el modo 3 queda:

( 0,025905394 ; - 0,017781023 ; 0,082672411 )

Cálculo del vector de fuerzas en cada piso:

Modo 1:

F1: ( 7,556500937 ; 16,51960606 ; 12,03789701) ton

Modo 2:

F2 ( 7,44051705 ; 7,44051705 ; - 7,44051705) ton

Modo 3:

F3: ( 9,904581835 ; - 6,798337044 ; 1,554879409) ton

Cálculo del cortante basal:

Vb1= 36,07160748 ton

Vb2= 7,44051705 ton
Vb3= 4,661124183 ton
Cálculo de los valores máximos probables:

T1= 13,7228 seg

T2= 30 seg
T3= 43,7228 seg
 a) Cortante por entrepiso por modo y valores maximos probables (usando hoja
de excel obtengo los valores):

Modo 1:

Vbasal: 36,071 ton

V1= 36,110 ton
V2= 28,557 ton
V3= 12,037 ton

Modo 2:

Vbasal: 7,440 ton

V1= 7,440 ton
V2= 0 ton
V3= -7,440 ton

Modo 3:

Vbasal: 4,661
V1= 4,661
V2= -5,243
V3= 1,554

Valores máximos probables:

Vbasal: 37,3104674 ton

V1= 14,914 ton
V2= 19,241 ton
V3= 14,192 ton

b) Desplazamiento en cada nivel, elástico y total:

u1= 0,20728538 cm
u2= 0,44077817 cm
u3= 0,6476303 cm

Multiplicando por el valor cd= 4,25 x u, (valor obtenido en la tabla 15 de la norma):

U1m= 0,880962865 cm
U2m= 1,873307223 cm
U3m= 2,72428775 cm
 c) Deriva lateral elástica y total en cada entrepiso

Los resultados son idénticos en la deriva lateral y elástica

Deriva piso 1= 0,880962865 – 0 = 0,8809 cm

Deriva piso 2= 1,873307223 – 0,880962865 = 0,9923 cm
Deriva piso 3= 2,72428775 – 1,873307223 = 0,850980527 cm

d) Momento flector en la base de la columna derecha

Calculando el momento flector a través del cortante basal:

Modo 1: 36,07/2= 18,035 ton

Modo 2: 7,44/2= 3,72 ton
Modo 3: 4,66/2= 2,33 ton

Introduciendo esos valores en la hoja de Excel antes mostrada obtengo el valor de

correlación de los 3 modos como respuesta máxima probable:

R=18,6552337 ton
Mmax = (18,6552337 x 3 ) / 2 = 27,98 ton-m

e) Fuerza cortante en la base de la columna derecha

Obtendré el valor a través del cortante basal del primer entrepiso dividido entre 2 para
cada columna, este valor lo obtuve haciendo uso de la hoja de Excel con el valor de la
respuesta máxima probable del cortante basal de los 3 modos:

R= 37,3104674 ton / 2 = 18,6552 ton

f) Fuerza axial en la columna derecha del tercer entrepiso

Obtendré el valor de la fuerza axial a través de la fuerza cortante del tercer nivel que
obtuve en la hoja de Excel como valor de respuesta máxima probable de cortante para
ese nivel

R= 14,095033 ton / 2 = 7,0475165 ton

Momento para la columna del tercer nivel: (7,0475165 x 3) / 2 = 10,5712 ton – m

Fuerza axial= ( 2 x m ) / L

Fuerza axial= ( 2 x 10,5712 ) / 6 = 3,52375 ton

 3 Seleccione las dimensiones de las columnas a fin de tener las rigideces del
pórtico indicadas en el enunciado del problema. Determine mediante el
método de análisis espectral y utilizando un programa computacional de su
preferencia, la respuesta máxima de la estructura ante el sismo indicado.

a) Determine los desplazamientos laterales para una carga estática de 180 t

en el nivel 3, y compare con el resultado del cálculo manual de la pregunta 1;
b) Determine los períodos y modos de vibración y compare con el resultado
del cálculo manual de la pregunta 1;
c) Determine los mismos valores de respuesta solicitados en la pregunta 2 y
compare con los resultados del cálculo manual. Comente.

a) Determine los desplazamientos laterales para una carga estática de 180 t en

el nivel 3, y compare con el resultado del cálculo manual de la pregunta 1;

A través del programa SAP 2000 procederé a encontrar los desplazamientos horizontales
ante la carga de 180 toneladas en el tercer nivel:

Para el cálculo de los valores de las masas de las losas multipliqué la diagonal principal
de la matriz de masa por 981 m/s²:

Ejemplo: 0,20 ton/cm/s x 981cm/s =196,2 ton

Comentario:

Los valores del desplazamiento U1 en los 3 niveles coincide con los valores obtenidos
por el método manual en la pregunta uno