Logo Studenta

Maldonado-Chio-capitulo3

¡Este material tiene más páginas!

Vista previa del material en texto

CAPÍTULO 3 
 
ESTRUCTURAS DE MÚLTIPLES 
 GRADOS DE LIBERTAD 
________________________________________________________________________________ 
 
Página 
 
3.1. ECUACIÓN DE MOVIMIENTO 111 
 
3.2. SISTEMA CON VIBRACIÓN LIBRE NO AMORTIGUADA 114 
3.2.1. FRECUENCIAS NATURALES Y MODOS DE VIBRACIÓN 114 
3.2.2. PROPIEDAD DE ORTOGONALIDAD DE LOS MODOS NORMALES 116 
3.2.3. CALCULO DE FRECUENCIAS NATURALES Y MODOS NORMALES UTILIZANDO 
 EL MÉTODO DE JACOBI 116 
 
3.3. MOVIMIENTO FORZADO 117 
3.3.1. MÉTODO DE SUPERPOSICIÓN MODAL 117 
 
3.4. RESPUESTA AL MOVIMIENTO DE LA BASE 120 
3.4.1. MÉTODO DE SUPERPOSICIÓN MODAL A LO LARGO DEL TIEMPO 120 
3.4.2. MÉTODO DE SUPERPOSICIÓN MODAL ESPECTRAL 125 
3.4.3. MÉTODO DE COMBINACIÓN DE LA RESPUESTA MODAL 127 
3.4.4. MOVIMIENTO AMORTIGUADO 128 
 
3.5. REDUCCIÓN DE MATRICES DINÁMICAS 128 
3.5.1. CONDENSACIÓN ESTÁTICA 128 
 
3.6. RESUMEN: SISTEMA DE VARIOS GRADOS DE LIBERTAD 130 
 
3.7. EJERCICIOS RESUELTOS 132 
 
3.8. EJERCICIOS PROPUESTOS 161 
116 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 117 
ESTRUCTURAS DE MÚLTIPLES 
 GRADOS DE LIBERTAD 
 
 
n aquellos casos donde el modelo matemático de un solo grado de libertad no provee la 
respuesta dinámica exacta, pues la estructura no puede asumir, durante su movimiento, una 
forma única de desplazamiento, se tiene que considerar el modelo como un sistema de 
varios grados de libertad. 
En este capítulo se estudia la respuesta dinámica de estos sistemas utilizando la solución modal 
o el método directo, bien sea cronológico o espectral. Se debe tener en cuenta que la solución 
modal, es aplicable, en general, sólo a sistemas que permanecen en el rango lineal de respuesta; en 
el caso contrario (rango inelástico de respuesta) se debe usar el método de integración de las 
ecuaciones de equilibrio dinámico. 
 
 
3.1. ECUACIÓN DE MOVIMIENTO 
 
Al plantear la ecuación de movimiento de un sistema de múltiples grados de libertad, es importante 
hablar del modelo estructural. Luego, un edificio se puede modelar como un sistema conformado por 
pórticos planos con una losa de entrepiso idealizado como un cuerpo infinitamente rígido para 
desplazamiento en su propio plano. Es decir, que la posición horizontal de cualquier punto dentro de 
la losa o diafragma es posible describirlo a partir de los desplazamientos ortogonales X, Y y un giro 
alrededor de un eje perpendicular al plano del diafragma, Z. Por consiguiente cada pórtico plano, a su 
vez, se puede modelar como una estructura con un solo grado de libertad horizontal por piso como se 
muestra en la Figura 3-1(b). 
 
Para este fin la estructura se debe suponer bajo estas condiciones: 
 
 1. La masa de la estructura está concentrada al nivel de los pisos. 
 2. Las vigas en los pisos no tienen deformaciones axiales debido a la rigidez infinita del 
diafragma; y por ello en cada piso existe un solo grado de libertad horizontal. 
 3. Se desprecia la deformación axial de las columnas, siempre y cuando la estructura sea poco 
esbelta, baja y larga. 
 
Bajo las anteriores condiciones al concentrar las masas a nivel de los pisos, el problema se 
transforma en uno de igual número de grados de libertad cuantas masas concentradas existan. Por 
ejemplo, un edificio de 3 pisos, tiene 3 grados de libertad, es decir, tres desplazamientos 
horizontales, uno a nivel de cada uno los pisos, ver Figura 3-1(a) y (b). El considerar las vigas 
rígidas, introduce el requisito de que las uniones entre las vigas y las columnas están fijas sin 
rotación. Bajo la tercera condición se garantiza que las vigas permanecen horizontales durante el 
movimiento de la estructura. 
 
Además, el edificio puede tener cualquier número de vanos y de igual manera se puede modelar 
como uno de un sólo vano; idealizando las columnas como una sola, con masas concentradas a la 
altura de los pisos, como se muestra en la Figura 3-1. 
E 
118 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
(a) Estructura con cuatro vanos (b) Grados de
masas concentradas
(c) Idealización con
m
horizontales y tres pisos
rigideces equivalentes
(d) Idealización con
libertad
a nivel de cada piso por piso
3
m2
m1
1
2
k
k
k3
 
 
Figura 3-1. Idealización del pórtico con masas concentradas a nivel de los pisos 
 
 
La masa mi corresponde a la masa concentrada del piso i. A su vez, recordando que el coeficiente de 
rigidez es la fuerza requerida para producir un desplazamiento relativo de magnitud unitaria entre dos 
pisos adyacentes. La rigidez para el modelo estructural corresponde a la rigidez equivalente por piso, 
es así como Ki es la rigidez equivalente del piso i. 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 119 
(a) Estructura de varios
grados de libertad
k
k
k
k
k
k
m
m
k
k
k
k
m
m
(b) Variables involucradas (c) Idealización
Estructural
(d) Diagrama de 
cuerpo libre
(d) Diagrama de 
cuerpo libreEstructural
(c) Idealización(b) Variables involucradas(a) Estructura de varios
grados de libertad
I. Pórtico de un solo vano
II. Pórtico de varios vanos
m u1 1
m u 2 2
12 2k (u - u )
k u1 1
F (t)2
F (t)1
2u
1u
m2
m1
2F (t)
1F (t)
2k
1k
m2
m1
2F (t)
1F (t)
2F (t)
F (t)
2k
1k
1m
1u
2m
2u
F (t)2
m u 2 2
m u1 1
1k u1
k (u - u )2 2 1
F (t)11
F (t)2
1F (t)
 
Figura 3-2. Sistemas de dos grados de libertad 
 
Luego, observando los sistemas estructurales, de dos pisos, mostrados en la Figura 3-2 analizamos 
sus diagramas de cuerpo libre y planteamos la ecuación de equilibrio para cada excitación arbitraria 
Fi(t), de la siguiente manera: 
 
Primer piso 
 0 = (t)F - )u - u( k - u k + u m 11221111  (3-1) 
0 = (t)F - u k - u )k + k( + u m 12212111  (3-1a) 
Segundo piso 
 0 = (t)F - )u - u( k + u m 212222  (3-2) 
0 = (t)F - u k + u k - u m 2221222  (3-2a) 
120 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
Ecuaciones que pueden escribirse usando matrices como: 
 
} F { = } U{ [K] + } U { [M]  (3-3) 
 
Es decir, 






=












−
−+
+












(t)F
(t)F
u
u
kk
kkk
u
u
m0
0m
2
1
2
1
22
221
2
1
2
1


 (3-3a) 
 
Es de resaltar como la matriz de masa es una matriz diagonal. 
 
 
3.2. SISTEMA CON VIBRACIÓN LIBRE NO AMORTIGUADA 
 
Es necesario e importante realizar el estudio de una edificación en vibración libre (aún cuando en muy 
pocas circunstancias este bajo esta condición) porque de esta forma se pueden determinar las 
propiedades dinámicas más importantes de la estructura, que son la frecuencia natural o período y los 
correspondientes modos de vibración. 
 
 
3.2.1. FRECUENCIAS NATURALES Y MODOS DE VIBRACIÓN 
 
Contamos, para vibración libre no amortiguada, con el siguiente sistema de ecuaciones de equilibrio: 
 
donde, las matrices [M] y [K] son las matrices de masa y rigidez respectivamente. 
 
Siendo la solución del sistema anterior de ecuaciones diferenciales simultáneas de la forma: 
donde: 
 ai = amplitud del movimiento de la coordenada i 
 n = número de grados de libertad 
 
 
 
 
 0 = {U} [K] + }U{ [M]  (3-4) 
 
n 1,2,3,..., = i α) -t (ωsen a = u ii (3-5) 
 
expresado en forma matricial: 
 
α) -t (ωsen {a} = {U} (3-6) 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 121 
Reemplazando la solución en la ecuación 3-4 tenemos: 
Eliminando sen (t - ) y reordenando términos tenemos: 
Como podemos ver estas ecuaciones de movimiento no dependen del tiempo, y corresponden a un 
sistema algebraico de n ecuaciones lineales homogéneas con n incógnitas, amplitudes ai, y 2 por 
determinar. La solución a esta formulación matemática conocidacomo problema característico, 
requiere que el determinante de la matriz del factor de {a} sea igual a cero, es decir, para que el 
sistema de ecuaciones simultáneas homogéneo (ecuación 3-7a) tenga solución no trivial, el 
determinante de la matriz de coeficientes debe ser cero: 
Donde: 
 : se denomina determinante característico del sistema de ecuaciones diferenciales simultáneas. 
 
Determinante que al resolver resulta ser una ecuación algebraica de grado n de la incógnita 2, que se 
satisface para n valores de 2. Esta ecuación corresponde a la ecuación característica del sistema o 
ecuación de frecuencias. Las n raíces de esta ecuación se denominan valores característicos o valores 
propios o eigenvalores. Estas son siempre reales y positivas, y su raíz cuadrada corresponde a las 
frecuencias naturales del sistema, en radianes por segundo. A la frecuencia más pequeña, se le denota 
por 1 y se le denomina frecuencia fundamental. 
Luego, con cada valor de p2 que satisface la ecuación característica del sistema, se reemplaza en la 
ecuación 3-7a y se obtiene un sistema de n ecuaciones del tipo: 
Resolviendo cada ecuación se obtiene para cada p2 un vector {a}p llamado vector característico, o 
modo de vibración o eigenvector; en términos de una constante de proporcionalidad arbitraria. Es 
decir, que los valores {a1, a2, a3… an}p son reales, pero no tienen un valor determinado, pues 
cualquier escalar p, i{a}p también es una solución del sistema de ecuaciones 3-9. Por lo tanto, 
para cada frecuencia p se tiene un vector {a}p que cuenta con una forma definida pero una 
amplitud arbitraria. 
 
 
 
 
 
 
 0 = α) -t (ωsen {a} [K] + α) -t (ωsen {a} [M] ω - 2 (3-6a) 
 
 0 = {a} [M] ω - {a} [K] 2 (3-7) 
o 
0 = {a} ] [M] ω - [k] [ 2 (3-7a) 
 
 0 = | [M] ω - [K] |Δ 2= (3-8) 
 
         n3...2,1,p0aMωK p2p ==− (3-9) 
122 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
3.2.2. PROPIEDAD DE ORTOGONALIDAD DE LOS MODOS NORMALES 
 
Las amplitudes de vibración {a} o desplazamientos en un modo normal de vibración corresponden 
solamente a valores relativos, como se mencionó anteriormente, que pueden determinarse asignando 
un valor de 1 a uno de los elementos del vector {a}p y definiendo los restantes n-1 términos en función 
de dicho valor, en forma única. Este proceso es llamado normalización y los vectores resultantes se 
denominan modos normales. Este no es el único proceso de normalización de los modos, en algunos 
casos es conveniente normalizar los modos con respecto a la masa del sistema: 
siendo φ la componente normalizada i del vector modal j. 
 
Con base en lo anterior la condición de ortogonalidad está dada, en general, por: 
siendo [ I ] la matriz identidad. 
 
Otra forma de expresar la ortogonalidad es: 
 
 
3.2.3. CÁLCULO DE FRECUENCIAS NATURALES Y MODOS NORMALES 
UTILIZANDO EL MÉTODO DE JACOBI 
 
Para un sistema con muchos grados de libertad, los cálculos numéricos requeridos por el método 
anteriormente expuesto (método directo) se vuelven inmanejables. Para esto existen métodos 
numéricos como el método de Jacobi, el cual calcula los valores característicos y los correspondientes 
vectores característicos de un sistema (García, 1998). 
 
a m 
a
 = Φ
2
 jk k
n
1=k
j i
j i

 (3-10) 
 
 
 1 = } { [M] } {
 ji donde 0 = } { [M] } {
i
T
i
j
T
i


 
  ]I[ = ][ [M] T  (3-11) 
 
 
 ω = } Φ { [K] } Φ {
 ji donde 0 = } Φ { [K] } Φ {
2
i
T
i
j
T
i

 
][ω = ][Φ [K] ][Φ 
2T
 (3-12) 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 123 
 
El método básico de Jacobi fue desarrollado para la solución del problema característico "estándar", 
esto es: 
 
 
3.3. MOVIMIENTO FORZADO 
 
En este numeral se estudiara el movimiento forzado en función de los modos normales de vibración, y 
su respuesta total mediante la superposición de las soluciones modales independientes, reduciendo el 
problema de encontrar la respuesta de un sistema con múltiples grados de libertad a la determinación 
de las respuestas de sistemas con un solo grado de libertad. 
 
 
3.3.1. MÉTODO DE SUPERPOSICIÓN MODAL 
 
Partiendo de las ecuaciones de equilibrio para el caso particular de un edificio con n niveles, 
modelado con n grados de libertad, como se muestra en la Figura 3-3, tenemos: 
 
 
(a) Edificio (c) Desplazamientos(b) Masas por piso
1m
mn-1
2m
nm
22m u 
u
1
k (u - u )
1k u
22
1
us
11m u
1
cuerpo libre
(d) Diagrama de 
k (u - u )
k (u - u )
k (u - u )
3 3
n-1n-1
n
u2
un-1
nu
2
n-2
n-1n-1m (u )
n n-1
nm (u )n
F (t)n
F (t)n-1
F (t)2
F (t)1
F (t)n
n-1F (t)
2F (t)
1F (t)
 
 
Figura 3-3. Sistema de n grados de libertad 
 
 
 
 } 0 { = } Φ { ] ω ]I[ - [K] [ 2 (3-13) 
 
124 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
Ecuaciones que expresadas en forma matricial corresponden a: 
es decir: 
Este sistema de ecuaciones dependientes (sistema acoplado) ecuaciones 3-3 se transforma en un 
sistema de ecuaciones independientes (sistema desacoplado), es decir, en un sistema en que cada 
ecuación contiene una sola incógnita en función del tiempo. Para ello es necesario expresar los 
desplazamientos u en función de cada modo de la siguiente manera: 
siendo: 
{U} = la respuesta dinámica de la estructura en desplazamientos. 
{z} = coeficientes que determinan la contribución de cada modo, los cuales corresponden a funciones 
generales del tiempo o a coordenadas generalizadas. 
[φ] =matriz modal del sistema. 
 
Introduciendo la anterior transformación de coordenadas en la ecuación de movimiento se obtiene: 
 
 (t)F )u - u( k - u k + u m 11221111 = (3-14) 
(t)F)u - (u k - )u - (uk + u m 223312222 = (3-15) 
(t)F )u - (uk + u m
 
n1-nnnnn =

 (3-16) 
 
Reordenando: 
(t)Fuk - u)k k( + u m 22 112111 =+ (3-14a) 
(t)Fuk - u )k(k + u k - u m 2332321222 =+ (3-15a) 
(t)F uku k - u m
 
nnn1-nnnn =+

 (3-16a) 
 














=


























−
−
+−
−+
+


























(t)F
(t)F
(t)F
u
u
u
kk0
k
kkk
kkk
u
u
u
m0
m0
00m
n
2
1
n
2
1
nn
n
322
221
n
2
1
n
2
1








 (3-17) 
 
 } F(t) { = } U{ [K] + } U { [M]  (3-3) 
 
     ZΦU = (3-18) 
 
          F(t)ZΦKZΦM =+ (3-19) 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 125 
 
Multiplicando lo anterior por el vector modal i transpuesto, [φ]Ti y aplicando las relaciones de 
ortogonalidad que dicen: 
Se obtiene: 
 
Tanto [ I ] como [ 2 ] son matrices diagonales, y por esto el sistema se desacopla, lo cual indica que 
se ha pasado de un sistema de n ecuaciones simultáneas a n ecuaciones independientes de un grado de 
libertad, del tipo: 
 
        F(t)Φzωz T2 =+ (3-23) 
 
 
Por lo tanto, las ecuaciones de movimiento pueden ser escritas por n ecuaciones desacopladas de la 
forma: 
 
donde: 
 
siendo n el número de grados de libertad. 
 
La ecuación 3-24 se soluciona aplicando metodologías como la expuesta en el capítulo 2 como la 
integral de Duhamel. 
 
Una vez se obtiene la respuesta a lo largo del tiempo de los grados de libertad generalizados, zi, la 
respuesta total es la superposición de la contribución de cada uno de los modos: 
 
       IΦM i
T
i = (3-20) 
 
       i2i
T
i ωΦKΦ = (3-21) 
                F(t)ΦZΦKΦZΦMΦ Ti
T
i
T
i =+
 (3-22) 
 
 [ I ] [ 2 ] 
 n),1,2,3,.... = (i (t)P = z ω + z ii
2
ii (3-24) 
 (t)F + ... + (t)F + (t)F = (t)P nni22i11ii  (3-25) 
 
=
==
n
1i
ii )(t)z}{Φ({z}][Φ{u}(3-26) 
 
126 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
3.4. RESPUESTA AL MOVIMIENTO DE LA BASE 
 
Un sistema puede estar sometido a un movimiento en la base, como es el caso de una edificación 
bajo la acción de un sismo. En dichas circunstancias, la respuesta del sistema puede evaluarse 
representando la acción por medio de un acelerograma, o mediante un espectro de respuesta. En la 
primera opción la respuesta se podrá evaluar a lo largo del tiempo, pero en la segunda, la respuesta 
corresponderá a los valores máximos. A continuación se expone en detalle cada una de estas 
opciones. 
 
3.4.1. MÉTODO DE SUPERPOSICIÓN MODAL A LO LARGO DEL TIEMPO 
 
Un edificio, modelado como una estructura simple, de n pisos excitado en su base se representa como 
se muestra en la Figura 3-4: 
 
(c) Desplazamientos
us
(a) Edificio
k2
k1
m1
(b) Idealización
Estructural
kn
n-1k
m2
mn-1
nm
22u
m u 
1
k (u - u )
1
22
1
k u
us
(d) Diagrama de 
cuerpo libre
11
1
m u
k (u - u )
k (u - u )
k (u - u )
33
n-1n-1
n
u2
un-1
un
2
n-2
n-1n-1
n-1n
nnm (u )
m (u )
 
 Figura 3-4. Sistema de múltiples grados de libertad sometido a un movimiento en la base 
 
Las ecuaciones de equilibrio sin tener en cuenta el amortiguamiento son: 
 Piso 
 
(n) 0 = )u - u( k + u m
. . . 
. . . 
(2) 0 = )u - u( k - )u - u( k + u m
 
(1) 0 = )u - u( k - )u - u( k + u m
1-nnnnn
23312222
122s1111



 (3-27) 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 127 
Reemplazando los desplazamientos por expresiones en función del desplazamiento del suelo us 
tenemos: 
donde y corresponde al desplazamiento relativo a la base de la estructura. 
 
Luego, introduciendo la ecuación 3-28 en 3-27 se obtiene: 
 
Las anteriores ecuaciones pueden ser escritas en la forma: 
 
siendo: 
{} = matriz cuyo tamaño depende de las componentes del acelerograma, (t)us ; es decir si este tiene 
2 componentes (x,y) el vector  tendrá un tamaño de 2xn, siendo n los grados de libertad de la 
estructura. En otras palabras, sus términos son iguales a 1 si la aceleración del terreno es colineal 
con la aceleración del grado de libertad que representa la matriz. Si no es colineal se coloca 0 a 
la componente correspondiente. 
(t)us = aceleración aplicada a la fundación del edificio, representada por un acelerograma. 
 
Determinando la matriz modal [φ] de la solución del problema característico: 
 
Introduciendo la transformación: 
y derivando dos veces contra el tiempo: 
 
 n) 3... 2, 1,=i (siendo u - u = y sii (3-28) 
 
u + y = u sii (3-29) 
 Piso 
 
(n) u m - = )y - y( k + y m 
. . . 
. . . 
(2) u m - = )y - y( k - )y - y( k + y m
 
(1) u m - = )y - y( k - y k + y m
sn1-nnnnn
s223312222
s11221111



 (3-30) 
 
 (t)u } γ{ [M] - = }y { [K] + } y { [M] s (3-31) 
 } 0 { = }y { |[M] ω - [K] | 2 (3-32) 
 
     zΦy = (3-33) 
 
128 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
Reemplazando las ecuaciones 3-33 y 3-34 en las ecuaciones diferenciales 3-30 y multiplicando por 
[]iT, se llega a: 
Aplicando la propiedad de ortogonalidad de los modos normales, obtenemos: 
 
donde  corresponde a los factores de participación o coeficientes de participación. 
 
Por consiguiente, se obtienen n ecuaciones independientes de un grado de libertad del tipo: 
donde i es el factor de participación del modo i, y corresponde a la fila i de la matriz [  ] obtenida, 
como se mencionó anteriormente por medio de: 
 
La solución de las ecuaciones 3-36 se realiza por medio de la integral de Duhamel o por métodos 
numéricos (Garcia, 1998; Oller, 1997). Una vez se obtienen los valores de zi, para cualquier tiempo t, 
se pueden obtener los desplazamientos de la estructura para ese instante. Luego la respuesta total es la 
superposición de las respuestas individuales de cada uno de los modos: 
Las fuerzas dinámicas inerciales que se producen en la estructura, en cada modo de vibración, pueden 
obtenerse multiplicando los desplazamientos de cada modo por la matriz de rigidez de la estructura: 
 
Luego, para el modo i las fuerzas en la estructura corresponden a la indicada en la Figura 3-5a. 
 
     zΦy  = (3-34) 
 
      (t)u } γ{ [M] Φ - = } {Z ][Φ [K] Φ + } Z {][Φ [M] Φ s
T
i
T
i
T
i 
 (3-35) 
 
      (t)u } γ{ [M] Φ - = } {Z ][Φ [K] Φ + } Z {][Φ [M] Φ s
T
i
T
i
T
i 
 (3-35a) 
 
 [I] [2] [] 
 n)1,2,...., = (i (t)u Γ = z ω + z sii
2
ii  (3-36) 
     γMΦΓ T−= (3-37) 
 
   (t)z}{φ{z}Φ{y}
n
1i
ii
=
== (3-38) 
 
   ii {y}k{F} = (3-39) 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 129 
(b) Cortante basal y
Fn i
n-1F i
3F i
2F i
1F i
momento de vuelco
(a) Fuerzas dinámicas
2 i
F
inerciales
1 i
F
n-1F
3 i
F
i
nF i
u2
u1
n-1u
u
un
3
h
h
h
h
h
n
n-1
3
2
1
 
 
Figura 3-5. Fuerzas del modo i 
 
Con las fuerzas dinámicas inerciales se define el cortante basal del modo i en el instante ti como: 
Ahora si se retoma el concepto de coeficiente de participación {} y si lo multiplicamos por []T se 
obtiene: 
Aplicando el principio de ( [A] [B] )T = [B]T [A]T a []T [M], se obtiene [M]T [] = [M] [], dado 
que [M] es simétrica, con lo cual la ecuación 3-42 se puede expresar como: 
 
=
=
n
1j
ii FjV (3-40) 
y el momento de vuelco 
 
=
=
n
1j
jji FhM (3-41) 
 
           γMΦΦΓΦ TTT = (3-42) 
 
130 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
A su vez, la masa total de la estructura esta dada por la contribución de cada masa individual: 
 
Nuevamente, aplicando el principio de ( [A] [B] )T = [B]T [A]T 
 
Entonces, la masa total en cada dirección principal corresponde a la suma de los cuadrados de los 
coeficientes de participación modal, i, en esa dirección. El valor de 2 de cada modo es llamada la 
masa efectiva modal y puede entenderse como la fracción de la masa total que se activa en ese modo 
al vibrar debido a la excitación en la base. 
 
En los casos donde los modos de vibración no son ortogonales, los coeficientes de participación se 
definen de la siguiente manera: 
 
Y la masa efectiva modal en este caso se obtiene, para cada modo, por medio de: 
            γγΦMΦΓΦ TT == (3-43) 
 
 [I] 
 
       γMγM TTOTAL= (3-44) 
 
y reemplazando [] en la ecuación anterior, se obtiene: 
      ( )     ΓΦMΓΦM TTTTOTAL= (3-45) 
          ΓΦMΦΓM TTTOTAL= (3-46) 
 
 [I] 
Luego 
         == 2iTTOTAL ΓΓIΓM (3-47) 
 
    
   


=
=
=−=
n
1j
2i
jj
n
1j
j
i
j
ii
Ti
i
)(Φm
)m(Φ
ΦMΦ
γMΦ
Γ (3-48) 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 131 
El concepto de masa efectiva modal se utiliza para definir el número mínimo de modos necesarios 
para describir la respuesta, ya que en los casos donde la estructura cuente con muchos grados de 
libertad los modos superiores se hacen muy pequeños y poco aportan a la respuesta. 
 
 
3.4.2. MÉTODO DE SUPERPOSICIÓN MODAL ESPECTRAL 
 
De acuerdo a lo visto anteriormente, las ecuaciones de movimiento para un sistema sometido a una 
excitación en su base se expresan mediante: 
Dado que se pueden definir los modos y frecuencias de la estructura, con base en sus propiedades para 
vibración libre, la solución de las ecuaciones diferenciales simultáneas se obtiene desacoplando el 
sistema mediante la utilización de la siguiente transformación de coordenadas: 
 
y derivando dosveces contra el tiempo: 
 
Reemplazando las ecuaciones 3-33 y 3-34 en 3-50 y premultiplicando por []iT, se obtiene: 
 
Dado que [I] y [2] son matrices diagonales el sistema se desacopla, lo cual indica que se tienen n 
ecuaciones independientes de un solo grado de libertad del tipo: 
 
    ( )
   


=
=






=−=
n
1j
2i
jj
2
n
1j
j
i
j
ii
2Ti
iefec
)(Φm
)m(Φ
ΦMΦ
γMΦ
m (3-49) 
 
          (t)sUγMyKyM  −=+ (3-50) 
 
     zΦy = (3-33) 
 
     zΦy  = (3-34) 
 
       (t)u } γ{ [M] Φ - = } {Z ][Φ [K] Φ + } Z {][Φ [M] Φ s
TTT
 (3-35a) 
 
 [I] [2] [] 
 n)1,2,...., = (i (t)u Γ = z ω + z sii
2
ii  (3-36) 
132 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
La solución de estas ecuaciones se puede realizar por medio de las metodologías propuestas para 
sistemas de un grado de libertad. Obtenidos los valores de {z(t)} se pueden calcular los 
desplazamientos a lo largo del tiempo. 
Luego, la diferencia entre el análisis modal a lo largo del tiempo y el análisis modal espectral se 
establece a partir de la definición de la acción sísmica. En un análisis a lo largo del tiempo la acción 
está representada por el acelerograma o registro del movimiento en aceleración a lo largo del tiempo, 
y en un análisis espectral la acción esta representada por un espectro de respuesta bien sea de 
desplazamiento o aceleración. 
De esta forma, como su mismo nombre lo dice, un análisis a lo largo del tiempo, presenta sus 
resultados cronológicamente, pero en un análisis espectral los resultados obtenidos se reducen a los 
máximos. 
Por consiguiente, el valor máximo que puede alcanzar el desplazamiento relativo y, entre la base 
y la masa de un sistema de un grado de libertad sometido a un acelerograma en su base (t)us , es el 
valor que se lee del espectro de desplazamiento Sd(T,), calculado para el mismo acelerograma. 
Por ello, el máximo valor que puede tener zi corresponde al valor leído del espectro de 
desplazamientos de la excitación amplificado por el coeficiente de participación i, luego: 
En los casos que se disponga de un espectro de aceleraciones el zi máx es igual a: 
 
Luego los desplazamientos dinámicos modales máximos correspondientes al modo i, pueden 
obtenerse mediante: 
 
A su vez, las fuerzas dinámicas inerciales modales máximas que se presenta en la estructura se pueden 
obtener multiplicando la matriz de rigidez de la estructura por el vector de desplazamientos máximos 
del modo correspondiente: 
 
 
Estas fuerzas pueden utilizarse como fuerzas estáticas y mediante un análisis estático tradicional, se 
pueden encontrar las fuerzas internas causadas por el modo i en cada uno de los elementos. 
 
 
   )ξ,(TSΓz iiidimáxi = (3-51) 
  
2
i
iiia
imáxi
ω
)ξ,(TS
Γz = (3-52) 
 
        )ξ,(TSΓΦzΦy iiidiimáxiii == (3-53) 
     ii yKF = (3-54) 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 133 
3.4.3. MÉTODO DE COMBINACIÓN DE LA RESPUESTA MODAL 
 
Una vez se han calculado las respuestas máximas, para cada modo, bien sea por el método espectral, 
se observa que éstas no ocurren al mismo tiempo. Luego para determinar la respuesta máxima 
tomando la contribución de todos los modos existen varios métodos, a continuación, se presentan dos 
de ellos (asumiendo un sistema de dos grados de libertad): 
 
1. Sumando los valores absolutos de las contribuciones modales máximas 
Como se mencionó anteriormente, generalmente cuando la respuesta de un modo llega a su máximo, 
las otras respuestas modales, en ese instante, no están en su máximo valor. Luego, es obvio, que este 
método de combinación ofrece el límite superior de respuesta. 
 
2. Por medio de la raíz cuadrada de la suma de los cuadrados de las contribuciones modales (RCSC) 
Los anteriores métodos de combinación modal pueden ser aplicados para calcular no solamente la 
respuesta máxima, sino también, para determinar la respuesta en un instante de tiempo determinado 
cuando el análisis se ha realizado a lo largo del tiempo: por ejemplo para calcular el desplazamiento 
del grado de libertad i en el instante t. 
 

=
=
n
1j
jiji (t)zΦ(t)U (3-55a) 
o 
( )
=
=
n
1j
2
jiji (tzΦ(t)U (3-36a) 
 
Estos métodos de combinación modal son aplicables para calcular cualquier tipo de respuesta, bien 
sea velocidad, aceleración o fuerzas. 
 
Existe otro método que toma en consideración la correlación que existe entre los diferentes 
modos y es el método de combinación cuadrática completa (CCC) (Rosenblueth y Elorduy, 1969). 
 
 
 
 
 
| z φ | + | z φ | = u
 
| z φ | + | z φ | = u
max 222max 121max 2
max 212max 111max 1
 (3-55) 
 
 
)z φ( + )z φ( = u
 
)z φ( + )z φ( = u
2
max 222
2
max 121max 2
2
max 212
2
max 111max 1
 (3-56) 
 
134 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
3.4.4. MOVIMIENTO AMORTIGUADO 
 
Cuando un sistema es amortiguado el procedimiento para desacoplar las ecuaciones diferenciales no 
es tan sencillo como hasta el momento se ha planteado, la amortiguación que existe generalmente en 
una estructura es muy pequeña, razón por la cual es la mayoría de los casos no es tenida en cuenta ni 
para el cálculo de las frecuencias y los modos, ni para el resto del planteamiento, reduciendo el 
problema al de una estructura sin amortiguamiento. Sin embargo en casos que sea conocido las 
ecuaciones de movimiento serán: 
 
 
3.5. REDUCCIÓN DE MATRICES DINÁMICAS 
 
En sistema donde se cuenta con un gran número de grados de libertad, las matrices de masa y rigidez 
saldrán de grandes dimensiones dificultando los procesos necesarios para llegar a la respuesta. En 
estos casos, es recomendable utilizar un procedimiento que reduce el tamaño de las matrices llamado 
condensación. Uno de los métodos sencillos utilizado en problemas dinámicos es la condensación 
estática, aunque es sólo una metodología aproximada. 
 
 
3.5.1. CONDENSACIÓN ESTÁTICA 
 
El método para reducir la matriz de rigidez parte de definir los grados de libertad que se desean 
condensar (dependientes o secundarios) y de expresarlos en función de los grados de libertad 
definidos como independientes o primarios. Quedando la matriz de rigidez dentro de la relación 
fuerza-desplazamiento, de la siguiente forma: 
Donde el subíndice s indica grados secundarios por reducir y el subíndice p designa los grados 
primarios. Es de mencionar que en la anterior ecuación se ha supuesto que no existen fuerzas externas 
aplicadas a los grados de libertad secundarios. 
Expandiendo las anteriores ecuaciones tenemos: 
De esta última relación (3-61) obtenemos: 
 (t)u Γ = z ω + z ω ξ 2 + z sii
2
iiii  (3-57) 
[F] = [U] [K] (3-58) 
 
   
   
 
 
 
 





=












0
F
U
U
KK
KK P
S
P
SSSP
PSPP
 (3-59) 
 
 } F { = } U { ]K[ + } U { ]K[ PSPSPPP (3-60) 
 
} 0 { = } U { ]K[ + } U { ]K[ SSSPSP (3-61) 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 135 
y reemplazando en (3-60)tenemos: 
Luego, la matriz de rigidez reducida o condensada será igual a: 
 } U { ]K[ ]K[ - = } U { PSP
-1
SSS (3-62) 
 
][K
 
} F { = } U { ] ]K[ ]K[ ]K[ - ]K[ [
 
} F { = } U { ]K[ ]K[ ]K[ - } U { ]K[
condensada
PPSP
1-
SSPSPP
PPSP
-1
SSPSPPP
 (3-63) 
 
 ]K[ ]K[ ]K[ - ]K[ = ][K SP
-1
SSPSPP
c
 (3-63a) 
 
136 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
3.6. RESUMEN: SISTEMAS DE VARIOS GRADOS DE LIBERTAD 
 
A continuación se presenta un resumen de la definición de la respuesta dinámica de una estructura de 
múltiples grados de libertad bajo la acción sísmica. 
us
 
Figura 3-6. Sistema sometido a una acción sísmicaSiendo la ecuación de movimiento: 
y si determinamos una matriz [] tal que cumpla lo siguiente: 
Podemos plantear la siguiente transformación de coordenadas: 
La cual reemplazando en la ecuación de movimiento: 
 } u { } γ{ [M] - = }y { [K] + } y { [M] s (3-31) 
       IΦM i
T
i = (3-11) 
 
       i2i
T
i ωΦKΦ = (3-12) 
  } {z = }y {  (3-33) 
 
  } z{ = } y {   (3-33a) 
 
  } z { = } y {   (3-34) 
 
   }u{ }{γ [M] - = {z} ][Φ [K] + }z{Φ [M] s (3-64) 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 137 
 
Premultiplicando por []T: 
la cual se reduce a: 
De aquí tenemos n ecuaciones independientes de la forma: 
Ecuaciones que incorporándoles los efectos del amortiguamiento quedan: 
El valor de zi se puede obtener por los métodos mencionados en el capítulo 2, pero al máximo valor se 
puede leer en el espectro de respuesta del sismo. 
donde Sd(Ti,i) es el valor del espectro de desplazamientos para un período Ti y un amortiguamiento i. 
 
Siendo la respuesta para cualquier tiempo: 
y la máxima respuesta: 
}u{ }{γ [M] ][ - = {z} ][ [K] ][ + }z { ][ [M] ][ s
TTT
  (3-35a) 
 
 [I] [2] 
 }u{ }{γ [M] ][Φ - = {z} ]ω[ + }z{ [I] s
T2  (3-35b) 
 
 [i] 
 u Γ = z ω + z sii
2
i  (3-36) 
 u Γ = z ω + z ω ξ 2 + z sii
2
ii  (3-65) 
 )ξ , T( S | Γ | = )z( iidimaxi (3-51) 
 (t)z | Γ | }{φ = {y(t)} iii
n
1=i
 (3-66) 
 
 )ξ , T( S | Γ | }{φ = }{y iidii
n
=1i
max  (3-67) 
ó 
] )ξ , T( S | Γ | }{φ [ = }{y
2
iidii
n
1=i
max  (3-67a) 
138 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
3.7. EJERCICIOS RESUELTOS 
 
En este numeral se analizará un edificio de dos niveles modelándolo como un pórtico plano, donde 
primero se trabaja con la masa y rigidez equivalente de dicho pórtico y segundo con la rigidez 
condensada. 
 
 
3.7.1. Ejemplo N° 1 
 
La estructura mostrada en la Figura 3-7 corresponde a una edificación apóticada construida en hormigón 
armado. Las vigas de la edificación tienen 0.35 m de altura por 0.30 m de ancho y, las columnas son 
cuadradas de 0.35 m de lado. 
Las masas de la estructura se estiman en 200 kg/m2 en la placa del primer nivel y de 100 kg/m2 
en la cubierta. La estructura cuenta con un módulo de elasticidad E = 20000 MPa. Es de interés 
determinar las frecuencias, los modos de vibración del edificio y la respuesta dinámica de la estructura 
cuando vibra libremente en la dirección X. 
 
Y
X
6.00 m
6.00 m
3.00 m
4.00 m
4.00 m 4.00 mZ
su
 
 
Figura 3-7. Estructura conformada por pórticos simétricos 
 
Solución. 
 
Para la determinación de las frecuencias y los modos de vibración se estudiará la estructura en la 
dirección X solamente, para ello se desarrollaron los siguientes pasos: 
 
 
 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 139 
1. IDEALIZACIÓN ESTRUCTURAL 
 
Debido a la uniformidad de las propiedades estructurales a lo largo del edificio y a que un mismo 
pórtico se repite cada 6 metros en la dirección X, el análisis de un pórtico interior puede representar la 
respuesta para todo el edificio. 
 
Luego, la estructura la podemos modelar como un edificio simple, es decir, suponemos las siguientes 
condiciones: 1) Toda la masa de la estructura se concentra a nivel de los pisos, 2) Las vigas en los 
pisos son infinitamente rígidas respecto a las columnas y 3) La deformación de la estructura es 
independiente de las fuerzas axiales en las columnas, ver Figura 3-8. 
 
2
Horizontales
(a) Pórtico con 2 vanos
ba
e
k
k
k
dk
u1
2u
m
22m u 
cuerpo libre
(c) Diagrama de 
1
(b) Idealización
Estructural
c
f
k
k
1m
u
1k
2k
m1
k u1 1
12 2k (u - u )
11m u 
u
m2
2
 
Figura 3-8. Fuerzas involucradas e idealización de la estructura 
 
Por lo tanto, la estructura puede representarse por un sistema de masas y resortes: 
 
u21u
2k1k
k u1 1
k (u - u )2 2 1
22m u m u1 1
m1 m2
 
 
Figura 3-9. Representación de la estructura como un sistema masa - resorte 
 
 
 
140 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
2. DETERMINACIÓN DE LA MASA 
 
Conocidas las masas distribuidas en el área de la placa intermedia y de la cubierta, se determinan los 
valores correspondientes al pórtico interior: 
 
m1 = 200x8x6 = 9600 kg = 9.60 Mg1 
m2 = 100x8x6 = 4800 kg = 4.8 Mg 
 
3. DETERMINACIÓN DE LA RIGIDEZ DEL SISTEMA 
 
Al suponer las vigas rígidas, la rigidez de las columnas entre pisos estará dada por 12EI/L3, luego 
asumiendo un módulo de elasticidad del hormigón E = 20000 MPa. 
 
4. PLANTEAMIENTO DE LAS ECUACIONES DE MOVIMIENTO 
 
Planteando la ecuación de movimiento de cada masa en vibración libre: 
Ordenando en forma matricial las ecuaciones de movimiento se obtiene: 
Siendo la respuesta para este sistema de ecuaciones la siguiente: 
 
 
1 Un megagramo (Mg) es igual a 1000 kg 
 
( )
( )
( )
( ) m
KN
33347.22
m3.00
m0.35
12
1
m
N2x1012
(3)k
m
KN
14068.36
m4.00
m0.35
12
1
m
N2x1012
(3)k
33
44
2
10
2
33
44
2
10
1
=






=
=






=
 
 
Primer piso 
 
0 = u k - u )k + k( + u m 
 0 = )u - u( k - u k + u m
2212111
1221111


 (1) 
Segundo piso 
 
0 = u k + u k - u m 
 0 = )u - u( k + u m
221222
12212


 (2) 
 






=












−
−+
+












0
0
u
u
kk
kkk
u
u
m0
0m
2
1
22
221
2
1
2
1


 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 141 
 
Derivando dos veces se obtiene: 
Reemplazando la respuesta en las ecuaciones de equilibrio se obtiene: 
 
 
Reordenando en forma matricial: 
 
Correspondiendo a un sistema algebraico de n ecuaciones lineales homogéneas con n incógnitas, y 
amplitudes ai y 2 por determinar. Para una solución no trivial de estas ecuaciones, se requiere que el 
determinante de la matriz de coeficientes sea igual a cero, esto corresponde a la ecuación característica 
del sistema. Para cada valor de 2 que satisface la ecuación característica del sistema, podemos 
 
α) - t (ωsen a = u 
 α) - t (ωsen a = u
22
11
 
 
α) - t (ωsen ω a - = u 
 α) - t (ωsen ω a - = u
2
22
2
11


 
Primer piso 
 0 = α) - t (ωsen a k - α) - t (ωsen a )k+ k( + α) - t (ωsen ω a m - 22121
2
11 
 0 = a k - a )k+ k( + ω a m - 22121
2
11 
 0 = a k - a )]k+ k( + ω m [- 22121
2
1 (3) 
Segundo piso 
 0 = α) - t (ωsen a k+ α) - t (ωsen a k - α) - t (ωsen ω a m - 2212
2
22 
 0 = akak ω a m - 2212
2
12 +− 
 0 = a )k + ω m (- + a k - 22
2
212 (4) 
 
       0aMωK 
0
0
a
a
m0
0m
ω
kk
kkk
0
0
a
a
kωmk
k)k(kωm
2
2
1
2
12
22
221
2
1
2
2
22
221
2
1
=−






=


















−





−
−+






=












+−−
−++−
 
 
142 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
resolver la ecuación de movimiento para a1, a2 en términos de una constante de proporcionalidad 
arbitraria. 
 
5. DETERMINACIÓN DE LAS FRECUENCIAS DE VIBRACIÓN 
 
Las frecuencias se pueden obtener haciendo que el determinante de la matriz del factor {a} sea igual a 
cero, es decir: 
Reemplazando los valores de masas y de coeficientes de rigidez respectivos tenemos: 
Las raíces de la ecuación cuadrática son: 
Por lo tanto las frecuencias naturales de la estructura son: 
 1 = 30.48 rad/seg (modo 1) 
 2 = 104.68 rad/seg (modo 2) 
 
y en ciclos por segundo: 
 f1 = 1/2 = 4.85 cps (modo 1) 
 f2 = 2/2 = 16.66 cps (modo 2) 
 
y los períodos naturales: 
 T1 = 1/f1 = 0.206 seg (modo 1) 
 
0 = k k + ω ]m )k + k( - k m [- + ω m m 
 
 0 = k k - k k + k k + ωk m - ωm )k + k( - ω m m 
 
 0 = )}k (- ) k {(- - ]} k + ω m [- )] k + k( + ω m {[(- 
 
0 
k + ω m-k-
k-)k + k( + ω m-
21
2
22121
4
21
222221
2
21
2
221
4
21
222
2
221
2
1
2
2
22
221
2
1
=
 
 
 0 = 469140696 + ω 547728.10 - ω 46.08
24
 
 
10957.31 = ω
 
 929.15 = ω
 
 
46.08*2
469140696*46.08*4 - 10)(547728. 547728.10+
 = ω
2
2
2
1
2
2
1,2

 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 143 
 T2 = 1/f2 = 0.060 seg (modo 2) 
 
6. DETERMINACIÓN DE LOS MODOS DE VIBRACIÓN 
 
Con las frecuencias 1 y 2 se vuelve a la ecuación (3) o (4) y se determinan los valores de a1 y a2. 2 
 
Para la primera frecuencia natural, reemplazando 1 en (3) obtenemos: 
 
 38496.89 a11 - 33347.22 a21 = 0 
 
Para la segunda frecuencia natural, reemplazando 2 en (3) obtenemos: 
 
 -57780.28 a12 - 33347.22 a22 = 0 
 
Siendo aij el valor de a en el grado de libertad i y en el modo de vibración j. 
 
Las ecuaciones anteriores solo pueden resolversen para valores relativos de a11 y a21 ó a12 y a22. 
Luego, se acostumbra asignar un valor definido a uno de ellos, por ejemplo a11=1, y se determina a21. 
Haciéndose el mismo procedimiento para los valores del modo 2. Es así, como: 
 
Para el primer modo: 
 a11 = 1 
 a21 = 1.15 
 
Para el segundo modo: 
 a12 = 1 
 a22 = -1.73 
 
Los modos obtenidos son: 
 
 
 
2 En este caso se reemplaza la frecuencia 1 en la ecuación (3), pero se hubiese podido reemplazar en la 
ecuación (4). 
1Para
a = 111 a = 112
2Para 
a = -1.73a = 1.1521 22
(a) Primer modo o
Modo Fundamental Modo Armónico
(b) Segundo Modo o
 
Figura 3-10. Formas modales para la estructura 
144 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
Por consiguiente, el movimiento total es la superposición de los modos: 
donde C1`, C2` y α1, α2 son constantes de integración, que se determinan de condiciones iniciales de 
movimiento, tales como: 
 
3.7.2. Ejemplo N° 2 
 
En el edificio anteriormente trabajado, dar la respuesta de los modos en forma normalizada y verificar 
la condición de ortogonalidad. 
 
Solución. 
 
Definiendo los modos normalizados como: 
 
Los factores normalizadores son: 
 
Luego, los modos normalizados son: 
 
Primer piso 
Segundo piso 
)α - t ω(sen a C + )α - t ω(sen a C = (t)u
 
 )α - t ω(sen a C + )α - t ω(sen a C = (t)u
2222
‘
21121
‘
12
2212
‘
21111
‘
11
 
 
 
u = (0)u u = (0)u
u = (0)u u = (0)u
022022
011011


 
 
 
a m 
a
 = 
} a { [M] } a {
a
 = φ
2
kjk
n
1=k
ij
j
T
j
ij
ij

 
 
 
4.895 = )3(4.8)(-1.7 + )(9.6)(1.0 :2Modo
3.994 = )(4.8)(1.15 + )(9.6)(1.0 :1Modo
22
22
 
 
 
0.3534- = 
4.895
1.73-
 = 0.2879 = 
3.994
1.15
 = 
0.2043 = 
4.895
1
 = 0.2504 = 
3.994
1
 = 
2221
1211


 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 145 
 
Siendo, la matriz modal del sistema: 
 
Verificando la ortogonalidad: 
 
 
 
   
















0.3534-0.2879
0.20430.2504
 = =
2221
1211
 
 
 
























10
01
 = 
0.3534-0.2879
0.20430.2504
 
4.80
09.6
 
0.3534-0.2043
0.28790.2504
 
[I] =][Φ [M] ][Φ
T
 
 
 
























10957.310
0929.15
 = 
0.3534-0.2879
0.20430.2504
 
33347.220
014068.36
 
0.3534-0.2043
0.28790.2504
 
][ω =][Φ [K] ][Φ
2T
 
 
146 
 
 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
MALDONADO & CHIO 
 
3.7.3. Ejemplo N° 3. 
 
Encontrar la respuesta dinámica en la 
dirección X, al sistema apórticado ilustrado 
en la Figura 3-11 si se asume como acción 
sísmica el espectro de respuesta propuesto 
para el territorio colombiano (NSR-98). La 
edificación se encuentra en la ciudad de 
Bucaramanga (zona de riesgo sísmico alto). 
Es una edificación normal, con un coeficiente 
de amortiguamiento con respecto al crítico 
igual a 5%. 
 
Figura 3-11. Estructura tipo Pórtico simétrica
 
Solución. 
 
Las unidades utilizadas en la realización de este problema son: para la masa el kg (kilogramo), para 
las fuerzas el N (Newton) y para las distancias el m (metro). 
 
Cada pórtico en la dirección X cuenta con las siguientes dimensiones, Ver Figura 3-12. 
 
 
 
 
 
La estructura cuenta con una masa de: 
✓ Cubierta = 550 kg/m2 
✓ Piso intermedio = 750 kg/m2 
 
 
 
 
 
 
 
 
 
 
Figura 3-12. Dimensiones pórtico interior en la dirección X 
 
3.00 m
3.00 m
4.00 m
0.3 m
0.3 m
0.25 m
0.3 m
0.3 m
0.3 m
0.3 m
0.25 m
X
Y
4.00 m
5.00 m
5.00 m
3.00 m
3.00 m
Z
X
Y
A
B
C
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 147 
1. IDEALIZACIÓN ESTRUCTURAL 
 
▪ El estudio o análisis del pórtico interior puede 
ser representativo para estudiar la respuesta 
en todo un edificio, esto debido a la 
uniformidad de las propiedades estructurales 
a lo largo de la edificación. 
 
 
[K]PORTICO * 3 = [K]ESTRUCTURA 
 
 
▪ Se asume la losa infinitamente rígida en su 
propio plano, esto hace a las vigas infinitamente rígidas axialmente. Luego, el modelo de 12 
grados de libertad se reduce a uno de 10. 
[K]
12x12
[K]
10x10
(b) Modelo de 10 GDL(a) Modelo de 12 GDL
 
Figura 3-14. Eliminación de los grados de libertad horizontales 
 
▪ Al no tener gran altura la edificación con respecto a su ancho, es decir, tener poca esbeltez, se 
pueden despreciar las deformaciones axiales en las columnas, reduciendo el modelo a uno de 6 
grados de libertad. 
[K]
6x6
(b) Modelo de 6 GDL
10x10
[K]
(a) Modelo de 10 GDL
 
Figura 3-15. Eliminación de los grados de libertad verticales 
(b) Grados de libertad(a) Modelo
del modelo
12 Grados
de libertad
 
Figura 3-13. Pórtico en la dirección X h 
148 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
▪ Al no existir cargas aplicadas en los grados de libertad rotacionales estos grados se pueden 
condensar. Por consiguiente, la matriz de rigidez de la estructura pasa de 6x6 a una condensada 
de 2x2. 
[K]6 X 6 → [K]C2 X 2 
 
 
2. DETERMINACIÓN DE LAS MASAS 
 
Trabajando con masas traslacionales tenemos: 
▪ Cubierta: m = área x masa = (4)(10) m2 (550) kg/m2 = 22000 kg 
▪ Entrepiso: m = área x masa = (4)(10) m2 (750) kg/m2 = 30000 kg 
 
Expresando en forma matricial: 
 
  kg
300000
022000
m0
0m
M
1
2






=





= 
 
Siendo: m2 = masa del nivel superior 
 m1 = masa del nivel inferior 
 
 
3. DETERMINACIÓN DE LA MATRIZ DE RIGIDEZ 
 
Para determinar la matriz de rigidez de toda la estructura se calcula la matriz de rigidez para un 
pórtico y luego se multiplica por el número de pórticos en la dirección de interés (en este caso X). 
 
Para determinar la matriz de rigidez de un pórtico tenemos: 
 
3.1 Numeración de nodos y elementos 
 
 
 
 
Dimensiones de las vigas 0.25x0.30 m 
 
Dimensiones de las columnas 0.30x0.30 m 
 
Módulo de elasticidad del concreto 20000 
MPa 
 
 
 
 
Figura 3-16. Numeración de nodos y elementos 
 
 
3.00 m
3.00 m
4.00 m
2
1
3
5 6
4
1 2
43
65
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 149 
3.2 Determinación de la matriz de rigidez de los elementos horizontales en coordenadas globales 
(elementos 1,2): 
 
Inercia = (1/12)(0.25)(0.30)3 = 5.625x10-4 m4 
Area = (0.25)(0.30) = 0.075 m2 
Longitud = 4 m 
 
m
N
1125000042187500562500042187500
421875021093750421875021093750
0037500000000375000000
5625000421875001125000042187500
421875021093750421875021093750
0037500000000375000000




















−
−−−
−
−
−
−
 
 
 
 
3.3 Determinación de la matriz de rigidez de los elementos verticales en coordenadas globales 
(elementos 3,4,5,6): 
 
Inercia = (1/12)(0.30)4 = 6.75x10-4m4 
Area = (0.30)(0.30) = 0.090 m2 
Longitud = 3 m 
 
m
N
1800000009000000900000009000000
0600000000006000000000
900000006000000900000006000000
9000000090000001800000009000000
0600000000006000000000
900000006000000900000006000000




















−
−
−−−
−
−
−
 
 
 
 
 
3.4 Ensamble de la matriz de rigidez de la estructura (Pórtico) 
 
Primero partimos de las expresiones de las matrices de cada elemento en función de sus 
submatrices: 
150 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
 
 
Elemento 1: 
6x62221
1211
kk
kk






 
 
Elemento 2: 
6x64443
3433
kk
kk






 
 
Elemento 3: 
6x63331
1311
kk
kk






 
 
Elemento 4: 
6x64442
2422
kk
kk






 
 
Elemento 5: 
6x65553
3533
kk
kk






 
 
Elemento 6: 
6x66664
4644
kk
kk






 
 
 
 
Luego, ensamblando las anteriores submatrices tenemos: (el exponente indica a cual elemento se 
refiere) 
 
m
N
k0k000
0k0k00
k0kkkkk0
0kkkkk0k
00k0kkk
000kkkk
66
6
64
6
55
5
53
5
46
6
44
6
44
4
44
2
43
2
42
4
35
5
34
2
33
5
33
3
33
2
31
3
24
4
22
4
22
1
21
1
13
3
12
1
11
3
11
1




















++
++
+
+
 
 
 
 
Reemplazando los valores de las submatrices, e introduciendo las condiciones de borde, la matriz de 
rigidez es: 
Estas submatrices son función de la numeración de 
los nodos i y j del elemento 
 
21
5 6
3
5
3
2
4
4
6
1
 Figura 3-17. Numeración de nodos y elementos 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 151 
 
152 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
3.5 Planteamiento de las ecuaciones de ligadura 
 
En la Figura 3-18 se representa el pórtico con todos los grados de libertad yij, correspondiendo i al 
nodo y j a la dirección del grado de libertad. Estos 12 grados de libertad se deben reducir a seis, dos 
horizontales (uno por nivel) y cuatro rotaciones. Para esto, se aplican las siguientes ecuaciones de 
ligadura: 
 
 
➢ Para las vigas: 
y11 = y21 
y31 = y41 
 
➢ Para las columnas: 
y12 = y22 = y32 = y42 = 0 
 
De las ecuaciones de ligadura definimos unos grados de 
libertad como independientes y otros como 
dependientes: 
▪ Grados de libertad independientes: y11 , y31 
▪ Grados de libertad dependientes : y21 , y41 
 
 
Figura 3-18. Pórtico con todos sus grados de libertad 
 
Para indicar lo anterior en la matriz de rigidez total, se suma la fila y luego la columna del grado de 
libertad dependiente en el independiente. Para el caso de los grados de libertad iguales a cero, 
sencillamente se anulan su fila y columna correspondiente. 
 
Finalmente el pórtico se reduce a: 
Figura 3-19. Pórtico con los grados de libertad reducidos 
 
A continuación se aplican las ecuaciones de ligadura, en la forma anteriormente expuesta: 
 ➢ 
y21
22y
y23
12y
y11
y13
32y
y31
y33
42y
y41
y43
y
y
y
31y
33
11y
13y
43
23
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 153 
Sumando las filas y21 a y1-1 y y41 a y31; después de haber eliminado las respectivas filas y columnas 
y12 , y22 , y32 , y42 se obtiene: 
 
m
N
y
y
y
y
y
y
472500000562500009000000900000000
562500004725000000090000009000000
0120000000120000009000000600000090000006000000
900000090000000029250000900000056250000
009000000900000056250000292500009000000
90000006000000900000060000009000000600000090000006000000
yyyyyyyy
43
33
31
23
13
11
4341333123211311




















−−−−
−
−
−−
 
 
 
Sumando las columnas y21 a y11 y y41 a y31, se reduce la matriz a: 
 
m
N
y
y
y
y
y
y
4725000056250000900000009000000
5625000472500000090000009000000
00240000009000000900000012000000
9000000090000002925000056250009000000
0900000090000005625000292500009000000
90000009000000120000009000000900000012000000
yyyyyy
43
33
31
23
13
11
433331231311




















−−−
−
−
−
 
 
 
 
Reordenando se obtiene la matriz correspondiente al sistema mostrado en la Figura 3-19: 
 
m
N
y
y
y
y
y
y
4725000056250009000000009000000
5625000472500000900000009000000
9000000029250000562500090000009000000
0900000056250002925000090000009000000
00900000090000002400000012000000
90000009000000900000090000001200000012000000
yyyyyy
43
33
23
13
31
11
433323133111




















−
−
−−−
−
 
 
154 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
 
3.6 Condensación de la matriz de rigidez para un pórtico. 
 
La anterior matriz la podemos nomenclar de la siguiente forma: 
 
 
 
   
   
 
 











=






S
P
4x4SS4x2SP
2x4PS2x2PPP
y
y
KK
KK
0
F
 
 
Despejando: 
 
{yS} = - [KSS]-1 [KSP] {yP} 
 
Y retomando el concepto de matriz de rigidez condensada: 
 
[K]C = [KPP] – [KPS] [KSS]-1 [KSP] 
 
Se desarrollan los pasos necesarios para su cálculo: 
 
 
m
N
2.29809x103.19951x107.53920x102.4343x10
3.19951x102.29809x102.4343x107.53920x10
7.53920x102.4343x103.80594x108.06812x10
2.4343x107.53920x108.06812x103.80594x10
K
8999
9899
9989
9998
1
SS














−−
−−
−−
−−
=
−−−−
−−−−
−−−−
−−−−
−
 
 
   
m
N
57430.0459440065110.13208901
57430.0459440065110.13208901
3740.2699210379970.22397702
3740.2699210379970.22397702
KK SP
1
SS












−−
−−
−
−
=−
−
 
 
    
m
N
7324858578.603954031586.50
3954031586.501146409188.80
KKK SP
1
SSPS 





−
−
=−
−
 
 
 
m
N
92719141421.36057968413.49
6057968413.498865590811.19
K PORTICO
C






−
−
= 
 
 
3.7 Determinación de la matriz condensada del sistema (matriz de rigidez de toda la estructura): 
 
De la matriz de rigidez condensada del pórtico llagamos a la matriz de rigidez de toda la estructura 
multiplicando por tres (3) la del pórtico. 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 155 
 
m
N
78157424264.188223905240.4
88223905240.496616772433.5
K ESTRUCTURA
C






−
−
= 
 
4. DETERMINACIÓN DE LAS ECUACIONES DE MOVIMIENTO 
 
Siendo la ecuación de movimiento de un sistema igual a: 
 
         s(t)UγMyKyM  −=+ 
 
Reemplazando los valores de masa y rigidez calculados, la ecuación de movimiento se expresa 
como: 
 
s(t)U
1
1
300000
022000
y
y
78157424264.188223905240.4
88223905240.496616772433.5
y
y
300000
022000
1
2
1
2 














−=












−
−
+












 
Siendo: 
311311
112112
yyyy
yyyy
==
==


 
 
5. DETERMINACIÓN DE LOS MODOS Y LAS FRECUENCIAS 
 
5.1. Determinación de las frecuencias de vibración 
 
Haciendo el determinante de la ecuación característica igual a cero, obtenemos las frecuencias. 
 
   
0885x105.71460522
ω660000000ω196x101.26333381ω005031730080988x109.63144657
0
ω30000857424264.1923905240.4
923905240.4ω22000016772433.6
0
300000
022000
ω
78157424264.188223905240.4
88223905240.496616772433.5
0MωK
14
4212214
2
2
2
2
=−
+−−
=
−−
−−
=





−





−
−
=−
 
 
Ordenando: 
010x103.91684135ω996x101.76650681ω660000000 142124 =+− 
 
Cuyas raíces son: 
156 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
rad/s2432.5599ω
rad/s243.9655ω
2
2
2
1
=
=
 
 
de donde las frecuencias y los periodos son iguales a: 
 
Modo 2  (rad/s) f (Hertz) T (s) 
1 243.97 15.62 2.49 0.40 
2 2432.55 49.32 7.85 0.13 
 
Siendo: 
f = /2 
T = 2/ 
 
 
5.2. Determinación de los modos de vibración 
 
Reemplazando los valores de 2 en la ecuación característica,hallamos los modos de vibración: 
 
       






=












−−
−−
=−
0
0
a
a
ω30000857424264.1923905240.4
923905240.4ω2200060167724.33.
0aMωK
1
2
2
2
2
 
 
 
✓ Para el modo 1: 
0a850105164.1a923905240.4
o
0a923905240.4a11405093.6
1121
1121
=+−
=−
 
 
 
Haciendo a21=1 y reemplazando en cualquiera de las dos expresiones tenemos: 
a21 = 1 
a11 = 0.4771 
 
✓ Para el modo 2: 
0a215552535.8a923905240.4
o
0a923905240.4a36743886.4
1222
1222
=−−
=−−
 
 
Haciendo a21=1 y reemplazando en cualquiera de las dos expresiones tenemos: 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 157 
a21 = 1 
a11 = -1.5371 
 
Por lo tanto los vectores de las amplitudes de vibración son: 
 
✓ Modo 1:  






=
0.4771
1
a 1 
 
✓ Modo 2:  






−
=
1.5371
1
a 2 
 
Siendo las formas de vibración para cada modo las indicadas en la Figura 3-20. 
 
(a) MODO 1 (b) MODO 2
 
 
Figura 3-20. Formas de vibración para cada modo 
 
5.3. Normalización de los modos 
 
Sabiendo que  es igual a: 

=
=
n
jK
2
KJK
ij
ij
am
a
Φ 
 
Se calculan los denominadores: 
✓ Modo 1: ( ) ( ) 169.790.477130000122000amam 22111
2
212 =+=+ 
✓ Modo 2: ( ) ( ) 304.761.537130000122000amam 22121
2
222 =−+=+ 
158 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
Luego: 
21 = 1/169.79 = 0.0058896 
11 = 0.4771/169.79 = 0.0028099 
 
22 = 1/304.76 = 0.0032813 
12 = -1.5371/304.76 = -0.0050436 
 
En forma vectorial los modos de vibración son: 
 
   






−
=






=
0.0050436
0.0032813
Φ
0.0028099
0.0058896
Φ 21 
 
Y en forma matricial éstos quedan: 
 
 
Piso1
Piso2
ΦΦ
ΦΦ
Φ
2Modo1Modo
er
do
1211
2221






=
 
Piso1
Piso2
0.00504360.0028099
0.00328130.0058896
2Modo1Modo
er
do






−
=
 
 
 
Siendo ij el valor modal en el grado de libertad i y en el modo j. 
 
 
 
6. DETERMINACIÓN DE LAS ECUACIONES DESACOPLADAS 
 
Vamos a transformar éste sistema de ecuaciones dependientes (sistema acoplado) en un sistema de 
ecuaciones independientes (sistema desacoplado), es decir, que cada ecuación contenga una sola 
incógnita en función del tiempo. Para esto, partimos de la siguiente transformación de coordenadas: 
 
 
 z][Φ }y{
z][Φ {y}
 =
=
 
 
Donde: {z} = coeficientes que determinan la contribución de cada modo 
 
Reemplazando en la ecuación de movimiento: 
 
           s(t)U1MZΦKZΦM  −=+ 
 
Multiplicando por []T: 
 
                s(t)U1MΦZΦKΦZΦMΦ TTT  −=+ 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 159 
 
 
Donde de acuerdo con las propiedades de ortogonalidad: 
 
      
      2T
T
ω
2432.550
0243.97
ΦKΦ
I
10
01
ΦMΦ
=





=
=





=
 
A su vez, de la ecuación de movimiento se determinan los coeficientes de participación modal de la 
siguiente forma: 
 
     
2Modo
1Modo
79.1194
213.8682
1
1
MΦ
T





−
=






−= 
 
y con ellos se calculan las masas modales efectivas como 2 
 
Modo i 

i 
2 % Masa total 
% Masa total 
(acumulada) 
1 - 213.8682 45739.61 87.961 87.961 
2 79.1194 6259.88 12.038 99.999 
 
 
Por lo tanto las ecuaciones desacopladas quedan de la forma: 
 
       
2)(Modos(t)U0.791194z2432.56z
1)(Modos(t)U213.8682z243.97z
s(t)U
79.1194
213.8682
z
z
2432.550
0243.97
z
z
s(t)UΓZωZI
22
11
2
1
2
1
2






=+
−=+





−
=












+






=+
 
 
donde s(t)U es el registro del acelerograma. 
 
Introduciendo el amortiguamiento en el análisis como: 
 
M
C
ωξ2 = 
 
160 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
Para el caso particular de una edificación, donde  = 5% los valores (2  ) son: 
 
( )( )
( )( ) 4.9320892432.560.052:2Modo
1.561954243.970.052:1Modo
=
=
 
 
Por lo tanto las ecuaciones quedan: 
 
s(t)U79.1194z2432.55z4.932089z
s(t)U213.8682z243.97z1.561964z
222
111


=++
−=++
 
 
Si la acción sísmica se representa con un acelerograma, la respuesta dinámica se determina a través 
de la solución de las anteriores ecuaciones, y esta puede evaluarse a lo largo del tiempo. En aquellos 
casos que se desea conocer las respuestas máximas, se puede utilizar como acción sísmica el 
espectro de respuesta y la forma de evaluarlas es la siguiente: 
 
 
7. DETERMINACIÓN DE LOS MÁXIMOS VALORES DE zi 
 
El valor máximo de zi para cada ecuación se puede obtener del espectro de desplazamiento, en 
función de cada frecuencia, siendo: 
 
 {y} = [] {z} 
 
donde: {y} = respuesta 
 [] = modos de vibración 
 {z} = es el valor leído del espectro de desplazamiento o del de aceleraciones por el coeficiente de 
participación, dicho de otra forma: 
 
( )
( )
i2
i
iia
máxiiiidmáxi
r*
ω
ξ,TS
zorξ,TSz == 
 
donde Sd corresponde al valor espectral de desplazamiento para el período Ti y el amortiguamiento 
i. 
 
Tomando el espectro elástico de aceleraciones, para un coeficiente de amortiguamiento crítico del 
5% dado por la NSR-98 (Norma Sismo Resistente de 1998), ver Figura 3-21, y de acuerdo con los 
siguientes parámetros: 
 
▪ Zona de riesgo sísmico alto (Bucaramanga) 
  Aa = 0.25 (Aa = aceleración pico efectiva para diseño) 
▪ Perfil del suelo sobre el cual se apoya la edificación S2 
  S = 1.20 (S = coeficiente de sitio) 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 161 
▪ Edificación de ocupación normal 
  I = 1.0 (I = coeficiente de importancia) 
 
T = 0.3 s
en cada dirección principal en planta
Para análisis dinámico, solo
modos diferentes al fundamental
0 Lc
aS = A I(1.0+5.0T)a
T = 0.48 S
igual al 5 por ciento del crítico
un coeficiente de amortiguamiento
Nota:
Este espectro está definido para
T = 2.4 S
T (s)
S =
T
S = 2.5 A I
S = A I
(g)
2
A I
S =
1.2 A S I
a
a
a
a
a a
aS
a a
 
Figura 3-21. Espectro elástico de Diseño NSR-98 
 
 
donde: 
 
L
a
a
LC
a
a
Caa
TTpara
2
IA
S
TTTpara
T
ISA1.2
S
TTparaIA2.5S
=
=
=
 
 
Además, para T menores a To en modos diferentes al fundamental 
Sa = Aa I (1.0 + 5.0 T) 
 
▪ Determinación de los periodos corto y largo 
 TC = 0.48 S = 0.48 (1.20) = 0.576 s 
 TL = 2.40 S = 2.40 (1.20) = 2.88 s 
 
162 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
✓ Para el modo 1 (Fundamental) 
12 = 243.97 rad2/s2 → T = 0.40 s < TC → Sa = 2.5 Aa I = 2.5(0.25)(1) = 0.625 (g) 
Sa1 = Sa x (g) = 0.625 (10) = 6.25 m/s2 
Sd1 = Sa1/12 = (6.25 m/s2) / (243.97 rad2/s2) = 0.025618 m 
 
✓ Para el modo 2 (Armónico) 
22 = 2432.55 rad2/s2 → T = 0.13s < TC → Sa = Aa I(1+5T) = (0.25)(1)[1+5(0.13)]= 0.4125 (g) 
Sa2 = Sa x (g) = 0.4125 (10) = 4.125 m/s2 
Sd2 = Sa2/22 = (4.125 m/s2) / (2432.55 rad2/s2) = 0.001696 m 
 
Luego, los valores máximos de zi son: 
 
zi máx = Sdi ri 
 
✓ Para el modo 1 (Fundamental) 
 z1 máx = (0.025618 m) (213.8682) = 5.47888 m 
 
✓ Para el modo 2 (Armónico) 
 z2 máx = (0.001696 m) (79.1194) = 0.13417 m 
 
 
8.DETERMINACIÓN DE LOS DESPLAZAMIENTOS MÁXIMOS DE LA ESTRUCTURA PARA 
CADA MODO 
 
Los desplazamientos máximos se determinan en función de zi máx, determinando primero los 
primarios o transnacionales y con ellos los secundarios o rotacionales. 
 
▪ Desplazamientos traslacionales 
 
    máxiiiP zΦy = 
 
✓ Para el modo 1 (Fundamental) 
 
( ) m
0.015395
0.032268
m5.47888
0.0028099
0.0058896
y
y
31
11






=






=






 
 
✓ Para el modo 2 (Armónico) 
 
( ) m
0.0006767
0.0004403
m0.13417
0.0050436
0.0032813
y
y
31
11






−
=






−
=





MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 163 
▪ Desplazamientos rotacionales o giros (desplazamientos condensados) 
 
      PSP
1
SSS yKKy
−
−= 
 
✓ Para el modo 1 (Fundamental) 
 
rad
4.9696x10
4.9696x10
3.0719x10
3.0719x10
0.015395
0.032268
57430.0459440065110.13208901
57430.0459440065110.13208901
3740.2699210379970.22397702
3740.2699210379970.22397702
y
y
y
y
3
3
3
3
*43
33
23
13














−
−
−
−
=


















−−
−−
−
−
=














−
−
−
−
 
 
(*) expresado en los grados de libertad 
 
 
✓ Para el modo 2 (Armónico) 
 
rad
2.7068x10
2.7068x10
2.8127x10
2.8127x10
0.0006767
0.0004403
57430.0459440065110.13208901
57430.0459440065110.13208901
3740.2699210379970.22397702
3740.2699210379970.22397702
y
y
y
y
5
5
4
4
*43
33
23
13














−
−
−
−
=






−












−−
−−
−
−
=














−
−
−
−
 
 
(*) expresado en los grados de libertad 
 
▪ Los desplazamientos totales son: 
 
 






=
S
P
y
y
y 
 
✓ Para el modo 1 (Fundamental) 
 
rad
rad
rad
rad
m
m
4.9696x10
4.9696x10
3.0719x10
3.0719x10
0.015395
0.032268
y
y
y
y
y
y
3
3
3
3
43
33
23
13
31
11




















−
−
−
−
=




















−
−
−
−
 
 
164 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
✓ Para el modo 2 (Armónico) 
 
rad
rad
rad
rad
m
m
2.7068x10
2.7068x10
2.8127x10
2.8127x10
0.0006767
0.0004403
y
y
y
y
y
y
5
5
4
4
43
33
23
13
31
11




















−
−
−
−
−
=




















−
−
−
−
 
 
 
En la Figura 3-22 podemos ver la deformada de la estructura para cada uno de los modos de 
vibración: 
 
(a) MODO 1
(Fundamental) (Armónico)
(b) MODO 2
 
 
Figura 3-22. Deformaciones totales de la estructura para cada modo 
 
 
9. DETERMINACIÓN DE LAS FUERZAS INTERNAS POR ELEMENTO 
 
Con los desplazamientos totales del pórtico calculados, se pueden calcular las fuerzas en el pórtico 
de la siguiente manera: 
 
    yKF = 
 
Y si se desean calcular las fuerzas internas por elemento: 
 
     
eee
yKF = 
 
 
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 165 
Por ejemplo, para el elemento 1 se procederá a calcular sus fuerzas internas: 
 
 
✓ Para el modo 1 (Fundamental) 
 
23
22
21
13
12
11
3
3
16
15
14
13
12
11
y
y
y
y
y
y
3.0719X10
0
0.032268
3.0719X10
0
0.032268
1125000042187500562500042187500
421875021093750421875021093750
0037500000000375000000
5625000421875001125000042187500
421875021093750421875021093750
0037500000000375000000
F
F
F
F
F
F




















−
−




















−
−−−
−
−
−
−
=




















−
−
 
 
mN
N
N
mN
N
N
51838.31
25919.16
0
51838.31
25919.16
0
F
F
F
F
F
F
16
15
14
13
12
11






















−
−
−
=




















 
 
1
25919.16 N
2
25919.16 N
51838.31 N m 51838.31 N m
 
Figura 3-23. Fuerzas internas del elemento 1 para el modo 1 
 
 
✓ Para el modo 2 (Armónico) 
 
23
22
21
13
12
11
4
4
16
15
14
13
12
11
y
y
y
y
y
y
2.8127x10
0
0.0004403
2.8127x10
0
0.0004403
1125000042187500562500042187500
421875021093750421875021093750
0037500000000375000000
5625000421875001125000042187500
421875021093750421875021093750
0037500000000375000000
F
F
F
F
F
F




















−
−




















−
−−−
−
−
−
−
=




















−
−
 
 
166 
 
 
MALDONADO & CHIO 
ESTRUCTURAS DE MÚLTIPLES GRADOS DE LIBERTAD 
mN
N
N
mN
N
N
4746.43
2373.22
0
4746.43
2373.22
0
F
F
F
F
F
F
16
15
14
13
12
11






















−
−
−
=




















 
 
4746.43 N m4746.43 N m
2373.22 N
1
2373.22 N
2
 
Figura 3-24. Fuerzas internas del elemento 1 para el modo 2 
Las fuerzas elásticas de diseño son las que consideran la contribución de todos los modos y se 
calculan de la siguiente manera: 
 
( )
=
=
n
1j
2
ijiT
FF 
 
 
( ) ( )
( ) ( )
( ) ( )
( ) ( )
mN
N
N
mN
N
N
52054.67
26027.39
0
52054.67
26027.39
0
4746.4351837.82
2373.2225918.97
00
4746.4351837.82
2373.2225918.97
00
F
F
F
F
F
F
22
2
2
22
22
22
22
16
15
14
13
12
11






















=






















−+−
+
+
−+−
−+−
+
=




















 
 
Para determinar las fuerzas internas en todos los otros elementos, se debe realizar iterativamente lo 
planteado en el paso 9 de este ejercicio.
 
 
MALDONADO & CHIO Capitulo3 – v1 13/08/2020 10:21 
ANÁLISIS SÍSMICO DE EDIFICACIONES 167 
 
3.8. EJERCICIOS PROPUESTOS 
 
La estructura consiste en pórticos no arriostrados de concreto reforzado, con placas aligeradas y 
armadas en dirección Z. La estructura está construida con un concreto cuyo módulo de elasticidad es 
EC = 20000 Mpa. Las dimensiones de las vigas y columnas son: 
 Viga 1 : 0.30x0.45 m 
 Viga 2 : 0.40x0.45 m 
 Columna :0. 90x0.45 m 
 Carga en los pisos : 1000 kg/m2 
 Carga en la cubierta : 700 kg/m2 
 
Z
5.00 m
4.00 m
3.00 m
5.00 m
X
6.00 m
Y
3.00 m
6.00 m
6.00 m
6.00 m
Vi
ga
 1
C
o
lu
m
n
a
Vi
ga
 2
 
 
Figura 3-25. Pórtico simétrico de concreto reforzado 
 
 
 
Determine: 
1) Las propiedades dinámicas de la edificación 
2) Las fuerzas sísmicas y el correspondiente efecto en cada elemento de la estructura

Continuar navegando