9

ANÁLISIS SÍSMICO POR

SUPERPOSICIÓN MODAL ESPECTRAL

9.1 ANÁLISIS SÍSMICO

Para lograr el objetivo del diseño estructural asísmico o antisísmico es indispensable

atravesar la etapa del análisis. Esta es, a su vez, posterior a la de estructuración y
determinación de las características elásticas y geométricas de la estructura,
incluyendo la distribución de sus masas. En general el análisis estructural consiste en
la determinación de los efectos que la solicitación aplicada demande de la estructura.
En el caso de los sismos hablamos del análisis sísmico. En este caso la solicitación o
carga sísmica está caracterizada por la norma local correspondiente y viene expresada
en términos de un espectro de diseño. Los efectos que se desean determinar consisten
las en fuerzas y deformaciones resultantes de la carga sísmica. Por fuerzas se entiende
de modo general, tanto fuerzas de distinto tipo: axiales, cortantes, como también
momentos flectores. Por deformaciones se entiende principalmente desplazamientos y
rotaciones de los entrepisos así como distorsiones relativas entre piso y piso.
La práctica actual mundialmente aceptada del diseño antisísmico considera que
las solicitaciones sísmicas sobre la estructura se determinan por medio de un análisis
elástico. Si bien la tendencia moderna incorpora criterios de comportamiento
inelástico como herramientas de disipación de energía, el análisis se hace sobre la
base de que la estructura y sus elementos no exceden su resistencia y mantienen su
forma inicial, hipótesis implícitas en el análisis estructural en el rango elástico. Desde
este punto de vista entonces, se cuenta con dos caminos contemplados en los códigos
de diseño: análisis estático o análisis dinámico.
2 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

El análisis estático reduce las acciones sísmicas a fuerzas estáticas equivalentes y

todo el análisis se hace considerando un sólo juego de fuerzas aplicado a la estructura
estáticamente. El edificio puede analizarse tri- o bi- dimensionalmente pero el análisis
sigue siendo estático y único. Por otro lado el análisis dinámico, también contemplado
en los códigos modernos de diseño sísmico, considera las características o
propiedades dinámicas de la estructura en la determinación de las fuerzas sísmicas y
en cada efecto particular que desee calcularse. Su aplicación, sin embargo, no ha
estado tan difundida hasta la década de los 80’s en vista de la complejidad del
cómputo involucrado y en la necesidad de disponer de máquinas para el cómputo y
procedimientos para la determinación de las propiedades dinámicas de la estructura
misma, sin mencionar el trabajo posterior para determinar y combinar los efectos
modales.
Con la disponibilidad y potencia de las computadoras modernas, principalmente
las personales, el análisis dinámico de edificios es la herramienta apropiada para la
determinación de las fuerzas sísmicas, v.g.: el análisis dinámico.
En edificaciones particularmente elevadas el análisis dinámico viene a ser la
única herramienta racional de análisis pues los métodos estáticos equivalentes se
tornan demasiado conservadores. La distribución de fuerzas máximas resultante a lo
alto del edificio es bastante diferente de la triangular supuesta en los códigos (ver la
Fig. 9.1.b).

Fig. 9.1 Resultados de un análisis dinámico para un edificio de 10 pisos

INGENIERÍA SISMORRESISTENTE
SECC. 9.2: SISTEMAS DE VARIOS GRADOS DE LIBERTAD: VIBRACIÓN FORZADA 3
En la Fig. 9.1.a se puede apreciar también, que los desplazamientos máximos de
cada piso tienen configuraciones que no responden a la de la hipótesis simplificatorias
del análisis estático equivalente. Asimismo cuando las características de la estructura
estimulan la contribución de modos adicionales al fundamental en la respuesta, se
puede estar subestimando peligrosamente efectos locales en los pisos bajos y en los
más altos. En realidad con la facilidad para realizar este tipo de análisis, tan difundidos
actualmente, tiene poco asidero el seguir utilizando procedimientos estáticos
equivalentes.

9.2 SISTEMAS DE VARIOS GRADOS DE LIBERTAD: VIBRACIÓN

FORZADA
Un edificio real es un sistema de varios grados de libertad (esto se vio en el
Cap..8). El establecimiento de las ecuaciones de equilibrio se desarrollará más
adelante [Podría consultar también Ref. 11-Cáp.3 y 4 ]. Estas ecuaciones de
movimiento para el sistema de varios grados de libertad, como se vio en el capítulo
anterior, tienen la siguiente forma:

MU  CU  KU  F (t ) (9.1)

El vector de fuerzas F  t  puede tener distintas variaciones en el tiempo, pero para
sistemas lineales elásticos siempre es posible expresar estas variaciones como una
superposición de términos de la forma F f (t ) . Por lo tanto la Ec. (9.1) puede
reemplazarse por una más simple:

M U + C U + K U = F f (t ) (9.2)
Donde F representa un vector independiente del tiempo que contiene las
magnitudes de las fuerzas aplicadas en correspondencia con cada grado de libertad (o
en cada piso si se trata de un pórtico plano)

9.3 MÉTODOS DE ANÁLISIS DINÁMICO

La respuesta dinámica de una estructura a una excitación sísmica (caracterizada
usualmente por un movimiento de la base) puede ser obtenida por cualquiera de los
tres métodos generales usados en la solución de sistemas de varios grados de libertad.
1) Integración directa en el tiempo de las ecuaciones de movimiento, resolviendo
simultáneamente las n ecuaciones diferenciales a través de un procedimiento
de integración paso a paso.
2) Solución directa en el campo de frecuencias, resolviendo nuevamente n
ecuaciones simultáneas.
3) Análisis Modal.

Dr. JAVIER PIQUÉ DEL POZO

4 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

De todos estos procedimientos el primero es el único medio riguroso para tomar en

cuenta comportamiento no lineal. Sin embargo, si se efectúa un análisis lineal será
necesario, definir una matriz C de amortiguamiento (ver Cáp.8). Si el análisis es
estrictamente no lineal, entonces la mayor parte de la disipación de energía será
automáticamente incorporada y la matriz C se hace innecesaria, ya que sólo
representará una pequeña cantidad de amortiguamiento a pequeñas amplitudes debida
a otras causas.
En el tercer procedimiento la solución en cada modo puede nuevamente llevarse a
cabo en el dominio del tiempo o en el de las frecuencias. Las soluciones en el campo
de frecuencias están siempre limitadas a sistemas lineales pero tienen la ventaja que
permiten considerar propiedades dependientes de la frecuencia (una condición
deseable en el caso de los suelos). El comportamiento nolineal puede ser simulado a
través de un procedimiento iterativo en que los valores de la rigidez y el
amortiguamiento son recalculados al final de cada análisis para igualar el nivel de
deformaciones obtenido. [Ref. 2]
El análisis modal es de lejos el procedimiento más usado en dinámica estructural.
Permite desacoplar las 3n ecuaciones diferenciales de movimiento, reduciendo el
problema a la solución de n ecuaciones independientes de 1 grado de libertad. En la
mayoría de los casos sólo algunos modos contribuyen significativamente a la respuesta
y por lo tanto ni siquiera tienen que resolverse los n sistemas simples.

9.4 DESCOMPOSICIÓN MODAL DE LAS ECUACIONES DE

MOVIMIENTO
La existencia de los modos como un espacio vectorial es extremadamente importante
ya que permite reducir la solución de un sistema de n grados de libertad a la solución
de n sistemas independientes de 1 GDL, desacoplando de esa manera las ecuaciones
de movimiento.
Suponiendo que al inicio se ha resuelto el problema de valores propios o
característicos para determinar las frecuencias naturales ωi y las correspondientes
formas de modo. Asimismo se supondrá que las formas de modo X i han sido
normalizadas con respecto a la matriz de masas de manera que el producto
X i M X i  1 (véase Cap. 8).
T

9.4.1 Descomposición Modal sin considerar Amortiguamiento

Al no considerar el amortiguamiento la Ec. (9.2) quedaría reducida de la siguiente
manera:

M U + K U = F f (t ) (9.3)

INGENIERÍA SISMORRESISTENTE
SECC. 9.4.1: DESCOMPOSICIÓN MODAL SIN CONSIDERAR AMORTIGUAMIENTO 5
Usando la propiedad de los modos, que permite expresar cualquier vector del
espacio vectorial, por ellos definido, como una combinación lineal de las formas
modales y ciertos coeficientes, se supondrá que la solución de las ecuaciones de
movimiento viene dada por:

n
U =  ai (t ) X i (9.4)
i 1

derivándola dos veces obtendríamos:

n
U =  ai (t ) X i (9.5)
i=1

Al sustituir las Ecs. (9.4) y (9.5) en (9.3) se tendría:

 n   n 
M   ai (t ) X i  + K   ai (t ) X i  = F f (t )
 i 1   i 1 
n
 [M X i ai (t ) + K X i ai (t ) ] = F f (t )
i=1

Al premultiplicar esta última ecuación por X Tj (para j  1,2,, n ), el cual es

independiente de “i ” , obtendríamos:

n
X Tj  [M X i ai (t ) + K X i a i (t ) ] = X Tj F f (t )
i=1
n (9.6)
 [X Tj M X i ai (t ) + X Tj K X i a i (t ) ] = X Tj F f (t )
i=1

Al aplicar las condiciones de ortogonalidad:

X Tj M X i = 0 para j  i pero si j i X iT MX i = 1 (9.7)

X Tj K X i = 0 para j  i pero si j i X iT K X i = ωi 2 (9.8)

en la Ec. (9.6), para “ j = i ” y teniendo además en cuenta las condiciones de

ortogonalidad, quedaría reducida así:

X iT M X i ai (t ) + X iT K X i ai (t ) = X iT F f (t ) (9.9)

Dr. JAVIER PIQUÉ DEL POZO

6 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

Al dividir la Ec. (9.9) entre X iT M X i :

X iT K X i X iT F
ai (t ) + ai (t ) = f (t ) (9.10)
X iT M X i X iT M X i

Como podrá notarse la Ec. (9.10) aún no esta simplificada del todo. Sin embargo al
observar que hay un término que involucra los modos y las matrices K y M, podríamos
pensar en hacer uso de una expresión ya demostrada en el capítulo anterior, dada por:

K X i   i2 M X i la cual al ser premultiplicada por X Tj , con j = i, queda de la

siguiente forma:

X iT K X i   i2 X iT M X i

realizando el despeje de la frecuencia se tendría:

X iT K X i
 i2  (9.11)
X iT M X i

Además, definiendo como factor de participación estática “  i ”, al término que

relaciona los modos y las matrices F y M, según la Ec. (9.10) éste sería:

X FT  Fj x ji
j 1
i  i
 (9.12)
X iT M X i  
n
 m j x ji
2

j 1

Cabe señalar que las Ecs. (9.10) ,(9.11) y (9.12) podrían reducirse aún más puesto
que se normalizaron los modos respecto a la matriz de masas, es decir X iT M X i  1 ,
según esto se tendría:
T T
ai (t ) + X i K X i ai (t ) = X i F f (t ) . (9.13)

 i2  X iT K X i (9.14)

i  X iT F (9.15)

Finalmente de las Ecs. (9.14) y (9.15) en (9.13) se obtiene:

INGENIERÍA SISMORRESISTENTE
SECC. 9.4.: DESCOMPOSICIÓN MODAL CONSIDERANDO AMORTIGUAMIENTO 7

ai (t ) +  i ai (t ) =  i f (t )
2
(9.16)

que representa las “ n ” ecuaciones modales de movimiento para un sistema

forzado sin amortiguamiento.
Es conveniente señalar que con frecuencia también se suele expresar “ U ” como
sigue:
n
U   d i (t )  i X i (9.17)
i=1

entonces, si la relacionamos con la Ec. (9.4), podremos apreciar con claridad que
ai )  d i (t ) i . Donde “ d i (t ) ” es el factor de participación dinámica (dependiente
(t
del tiempo) y “  i ” es el factor de participación estática (independiente del tiempo).
Luego la Ec. (9.16) en función de “ d i (t ) ”quedaría expresada como:

di (t ) +  i d i (t ) = f (t )
2
(9.18)

9.4.2 Descomposición Modal considerando Amortiguamiento

Considerando amortiguamiento la Ec. (9.19) sería la misma que la Ec. (9.2), es decir:

M U + C U + K U = F f(t) (9.19)

De manera similar a la sección anterior, en ésta, se hará uso de la propiedad de los

modos, que permite expresar cualquier vector del espacio vectorial, por ellos definido,
como una combinación lineal de las formas modales y ciertos coeficientes, se
supondrá para ello que la solución de las ecuaciones de movimiento viene dada por:

n
U =  ai (t ) X i (9.20)
i 1

derivando una vez obtendríamos:

n
U =  a i (t ) X i (9.21)
i=1

derivando dos veces obtendríamos:

n
U =  ai (t ) X i (9.22)
i 1

Dr. JAVIER PIQUÉ DEL POZO

8 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

Al sustituir las Ecs. (9.20), (9.21) y (9.22) en (9.19), es decir, sustituyendo este
vector U y sus derivadas U y U , expresadas en función de las formas modales X i ,
las ecuaciones de movimiento que obtendríamos serían:

 n   n   n 
M   ai (t ) X i  + C   a i (t ) X i   K   a i (t ) X i  = F f (t )
 i 1   i 1   i 1 

n
 [M X i ai (t ) + C X i a i (t ) + K X i a i (t ) ] = F f (t )
i=1

Luego al pre multiplicar esta última ecuación por el vector X Tj ( para

j  1,2, , n ), el cual es independiente de “ i ” , obtendríamos:

n
X Tj  [M X i ai (t ) + C X i a i (t ) + K X i a i (t ) ] = X Tj F f (t )
i=1
n
(9.23)
 [X Tj M X i ai (t ) + X Tj C X i a i (t ) + X Tj K X i ai (t ) ] = X Tj F f (t )
i=1

Aplicando las condiciones de ortogonalidad:

X Tj M X i = 0 para j  i pero si j i X iT MX i = 1 (9.24)

X Tj CX i = 0 para j  i pero si j i X iT C X i = 2  i  i (9.25)

( Si C tiene una forma especial )

X Tj K X i = 0 para j  i pero X iT K X i = ωi 2 (9.26)

en la Ec. (9.23), para “ j = i ”, ésta quedaría reducida así:

X iT M X i ai (t ) + X iT C X i a i (t ) + X iT K X i ai (t ) = X iT F f (t ) (9.27)

Escribiendo de otra manera la Ec. (9.27) se tiene:

M ie ai (t ) + Cie a i (t ) + K ie ai (t ) = Fi e f (t ) (9.28)

siendo M ie , Cie , K ie y Fi e escalares, correspondientes a cada modo de vibración

“ i ”. Luego, al dividir la Ec. (9.28) entre M ie resulta:

INGENIERÍA SISMORRESISTENTE
SECC. 9.4.2: DESCOMPOSICIÓN MODAL CONSIDERANDO AMORTIGUAMIENTO 9
C ie K ie Fi e
ai (t ) + a i (t ) + ai (t ) = f (t ) (9.29)
M ie M ie M ie

Como podrá observarse, a la Ec. (9.29), que representa las “n” ecuaciones modales
del movimiento, se le puede hacer una analogía para el caso en el que solo se tiene 1
GDL. Entonces tendríamos lo siguiente:

Cie Cie c c
i (%)   es equivalente a  (%)  
C e
i ( crítico ) 2M iei ccrítico 2m
Ce X T C Xi X iT C X i (9.30)
  Ti  2 ii  2 ii 
Me Xi M Xi X iT M X i

También, recordando que se demostró en la sección anterior:

X iT K X i
 i2  (9.31)
X iT M X i

y que además, el término que involucra a los modos y a las matrices F y M,

llamado factor de participación estática “  i ”, estaba dado por:

X iT F
i  (9.32)
X iT M X i

Y debido a que se normalizaron los modos respecto a la matriz de masas, o sea

X iT M X i  1 , las Ecs. (9.27), (9.30), (9.31) y (9.32) podrían reducirse, según ello
estas quedarían así:
T T T
ai (t ) + X i C X i a i (t ) + X i K X i ai (t ) = X i F f (t ) . (9.33)

2 i i  X iT C X i (9.34)

 i2  X iT K X i (9.35)

i  X iT F (9.36)

Entonces, finalmente de las Ecs. (9.34), (9.35) y (9.36) en (9.33) se obtiene:

ai (t) + 2 i  i a i (t)   i ai (t) =  i f (t)

2
(9.37)

Dr. JAVIER PIQUÉ DEL POZO

10 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

que representa las “ n ” ecuaciones modales de movimiento para un sistema

forzado considerando amortiguamiento.
Análogamente a la sección anterior, “ U ” se suele expresar como:

n
U   d i (t)  i X i (9.38)
i=1

se puede apreciar con claridad que ai (t)  d i (t)i al relacionar la Ec. (9.38) con la
Ec. (9.20). Donde “ d i (t) ” , como ya se indico, es el factor de participación dinámica
(dependiente del tiempo) y “  i ” el factor de participación estática (independiente del
tiempo). Entonces la Ec. (9.37) quedaría expresada en función de “ d i (t) ”como sigue:

di (t) + 2 i  i d i (t)   i d i (t) = f (t)

2
(9.39)

De las Ecs. (9.38) y (9.39) observamos que la contribución de cada modo X i a la

respuesta está afectada por el factor de participación estática  i y un factor de
participación dinámica d i (t) que resulta de la solución de una ecuación de un sistema
de un grado de libertad con la frecuencia natural  i sometida a la función del tiempo
f t  .

Si la distribución de fuerzas dinámicas F (o para fuerzas estáticas) es

proporcional en cada masa al producto de la masa por su desplazamiento en el modo
j , es decir F  M X j , únicamente  j  X Tj F no será igual a cero y por
consiguiente sólo el modo j será excitado. El sistema vibrará manteniendo constante
la forma del modo j , o sea X j , variando solamente su amplitud, que dependerá de la
función f t  . En la mayoría de los casos prácticos el factor de participación estática
 i tiende a decrecer para los modos más altos, es decir aquellos con valores altos de
frecuencias.
La importancia relativa del factor de participación dinámica d i (t ) para cada modo
será una función de la variación de f t  con el tiempo en relación con la frecuencia
natural  i . Nuevamente, en general, las frecuencias más altas tendrán menor
amplificación y como resultado, la contribución de los modos altos en la respuesta no
será tan significativa. En la mayoría de casos prácticos solamente algunos modos (3 a
5 a lo más) serán suficientes para obtener una respuesta apropiada. (ello considerando
el problema plano; sin embargo si el problema se modela tridimensionalmente habrá
que triplicar este número).

INGENIERÍA SISMORRESISTENTE
SECC. 9.5: ANÁLISIS MODAL PARA EXCITACIONES SÍSIMICAS 11
La determinación de d i ( t ) requiere la solución de la ecuación de movimiento
para un sistema de 1 GDL. Esta puede efectuarse tanto en el campo del tiempo como
en el campo de frecuencias.
Debe tenerse en cuenta que la aplicación del análisis modal requiere no solamente
que el problema sea lineal (ya que está basado en la superposición) sino también la
existencia de una matriz de amortiguamiento C apropiada que satisfaga la condición
de ortogonalidad. Si se usa un modelo de acoplamiento cercano (ver Cáp. 8) cada
masa estará conectada a la superior e inferior por un amortiguador y la matriz de
amortiguamiento tendría una forma similar a la de la matriz de rigidez:

c1  c2  c2 0 0 .. 0 
 c c2  c3  c3 0 .. 0 
 2
 0  c3 c3  c4  c4 .. 0 
C 
 : : : : : 
 0 0  cn 1 cn 1  cn  cn 
 
 0 0 0  cn cn 
(9.40)

Otra forma de hacerlo es calcular una matriz de amortiguamiento con la siguiente

expresión

C = M Q B QT M (9.41)

donde M es la matriz de masas, Q la matriz modal (conteniendo todas las formas

modales como columnas) y B es una matriz diagonal cuyo término iésimo es igual a
2  ii (ver Cáp. 8).

En la mayoría de los casos, cuando se usa análisis modal, la matriz de

amortiguamiento ni siquiera se ensambla, sino que se comienza definiendo los
porcentajes modales de amortiguamiento  i y se los incorpora directamente en las
ecuaciones modales, Ec. (9.39).

9.5 ANÁLISIS MODAL PARA EXCITACIONES SÍSMICAS

En esta sección veremos cuando un sistema de varios grados de libertad está sometido
a una excitación sísmica, la que es representada usualmente como una aceleración
horizontal en la base.
Por simplicidad, para un mejor entendimiento de la expresión general,
demostraremos la expresión general basándonos en un sistema de vibración libre de

Dr. JAVIER PIQUÉ DEL POZO

12 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

2 GDL dinámicos en el que no se considerará el amortiguamiento. Además, en dicho

sistema se indicarán los desplazamientos absolutos “ u ” y relativos “ y ” (Fig. 9.2.a).

u2
u1

y2
y1

1 2

m 2 u2
m2 m2 m2
k 2 2  k2 ( y2  y1 )
k2 k2
k 2 2  k2 ( y2  y1 )
m1u1
m1 m1 m1
k11  k1 y1
k1 k1

u G (t )
Fig.9.2.a Fig.9.2.b
Sistema Movimiento
simplificado de la Base

Fig.9.2 (a) Sistema simplificado no forzado de 2 GDL dinámicos

(b) Movimiento de la base debido a una excitación sísmica.

De la Fig. 9.2.b aplicando equilibrio dinámico para el primer y segundo nivel,

en ese orden, resulta:

m1u1  k1 y1  k 2 ( y 2  y1 )  0  m1u1  (k1  k 2 ) y1  k 2 y 2  0 (9.42)

m2 u2  k 2 ( y 2  y1 )  0  m2 u2  k 2 y1  k 2 y 2  0 (9.43)

Como se pude observar las Ecs. (9.42) y (9.43) están en función de

desplazamientos absolutos “ u ” y desplazamientos relativos a la base “ y ”. Entre lo
absoluto y relativo podría optarse por escoger cualquiera de los dos. Sin embargo

INGENIERÍA SISMORRESISTENTE
SECC. 9.5: ANÁLISIS MODAL PARA EXCITACIONES SÍSMICAS 13
optaremos por trabajar con los desplazamientos relativos a la base, lo cual es
conveniente, puesto que las Ecs. (9.42) y (9.43) quedarán expresadas de una forma
ya tratada en el Cap. 8. El fundamento de lo dicho nuevamente será dado en breve
una vez que se obtengan las ecuaciones de movimiento en función de “ y ”.
Observando la Fig. 9.2.b, vemos que los desplazamientos absolutos y lo
relativos a la base están relacionados mediante:

ui  uG (t )  yi (9.44)

donde, para nuestro caso, “ i ” va de 1 a 2, puesto que estamos analizando un

sistema de 2 GDL dinámicos. Derivando dos veces la Ec. (9.44) tenemos:

ui  uG (t )  yi (9.45)

Al reemplazar las Ecs. (9.44) y (9.45) en (9.42) y también en (9.43) se tiene:

m1 (uG (t )  y1 )  k1 y1  k 2 ( y 2  y1 )  0

m2 (uG (t )  y2 )  k 2 ( y 2  y1 )  0

luego, al reordenar estas ecuaciones se tiene:

m1 y1  (k1  k 2 ) y1  k 2 y2  m1uG (t ) (9.46)

m2 y2  k 2 y1  k 2 y 2  m2 uG (t ) (9.47)

Ya reordenadas es fácil darse cuenta ahora que, como ya se dijo, fue

conveniente colocar las ecuaciones en función de los desplazamientos relativos a la
base ya que ecuaciones similares fueron tratadas en el Cap. 8, solo que en este caso
la fuerza, es decir el término  m i uG (t ) , depende de la masa “ m i ”y de la
aceleración del suelo o de la base “ uG (t ) ” . Expresado de otra manera, podemos
decir que las Ecs. (9.46) y (9.47) tienen por vector fuerza a “ F(t) ”dado por:

 P f (t )    m1 uG (t ) 
F (t )   1  
P2 f (t )  m2 uG (t )

Entonces, un sistema equivalente al sistema libre de la Fig. 9.2 vendría a estar

dado por el sistema forzado que se muestra en la Fig. 9.3:

Dr. JAVIER PIQUÉ DEL POZO

14 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

y2
y1

1 2
m 2 y2
m2 P2 f (t )   m 2 uG (t )

k2 k 2 2  k2 ( y2  y1 )
k 2 2  k2 ( y2  y1 )
m1 y1
m1 P1 f (t )   m1 uG (t )
k11  k1 y1
k1

Fig. 9.3 Diagrama de cuerpo libre ( DCL ) del sistema simplificado forzado sin
amortiguamiento, expresado en desplazamientos relativos a la base “ y ”

Ordenando matricialmente las Ecs. ( 6.46) y ( 6.47 ) se tiene:

m1 0   y1  k1  k 2  k 2   y1   m1 uG (t ) 

0     
 m2   y2    k 2 k 2   y 2  m 2 uG (t )

La ecuación anterior se suele escribir de la siguiente manera:

m1 0   y1  k1  k 2  k 2   y1  m1 0  1

0        uG (t )
m2   y2    k 2  m2  1
(9.48)
 k 2  y 2  0

su notación matricial de una manera mas concisa sería:

M Y  K Y   M I uG (t ) (9.49)

donde:

INGENIERÍA SISMORRESISTENTE
SECC. 9.5: ANÁLISIS MODAL PARA EXCITACIONES SÍSMICAS 15

 y  y  1
Y   1  , Y  1 e I  
 y2  y2  1

Son los vectores aceleración y desplazamiento relativos a la base, y el vector

columna 1 , en ese orden; además:

m 0 k  k  k2 
M  1 K  1 2
m2  k 2 
y
0   k2

son la matriz masa y de rigidez respectivamente.

Una expresión más general, para el sistema forzado con amortiguamiento de
2 GDL dinámicos que se muestra a continuación:

P2 f (t )   m2 uG (t )
m2

c2 k2

P1 f (t )  m1uG (t )
m1

c1 k1

Fig. 9.4 Sistema forzado con amortiguamiento de 2 GDL dinámico (con fuerzas
que dependen de las masas y de la acleración de la base), el cual es la
equivalencia del problema original mostrado en la “ Fig. 9.2 ”(un
sistema de vibración Libre de “ n ” GDL con amortiguamiento cuando
está sometida a una aceleración en el suelo o la base).

1
Se debe de tener bien claro, para lo concerniente al tema, que I es un vector columna cuyos elementos
son todos unos. No debemos confundirlo con la matriz identidad.

Dr. JAVIER PIQUÉ DEL POZO

16 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

es la siguiente expresión:

M Y  C Y  K Y =  M I uG ( t ) (9.50)

Se debe enfatizar que la Ec. (9.50), expresada en desplazamientos relativos,

representa la ecuación para cuando se tiene un sistema de vibración libre de “ n ” GDL
con amortiguamiento cuando está sometido a una aceleración en la base (ver Fig. 9.2 ,
9.3 y 9.4).
Además, se debe recordar al lector que en el Cap. 8 se explicó como es que se
forma la matriz de amortiguamiento “ C ”. Además se muestran las condiciones que
debe cumplir dicha matriz para poder ser incluida en la Ec. (9.49) para finalmente
tomar la forma de la Ec. (9.50).

De la Ec. (9.50) con Y , Y e Y los vectores de desplazamiento, velocidad y

aceleración relativos a la base ( Y  U  I uG ) e I vector cuyos elementos son
todos iguales a la unidad, y uG (t) la aceleración del suelo, procederemos a aplicar
descomposición modal presentada en la sección anterior.
Basados en la forma que tiene la Ec. (9.50), de forma análoga que en secciones
anteriores, se hará uso de la propiedad de los modos, puesto que, como ya se sabe, nos
permite expresar cualquier vector del espacio vectorial, por ellos definido, como una
combinación lineal de las formas modales y ciertos coeficientes. Para tal propósito se
supondrá que la solución de las ecuaciones de movimiento viene dada por:
n
Y =  ai (t ) X i
i 1 (9.51)
derivando una vez obtendríamos:
n
Y =  a i (t ) X i
i=1 (9.52)
derivando dos veces obtendríamos:
n
Y =  ai (t ) X i
i 1 (9.53)
al ser sustituidas las Ecs. (9.51), (9.52) y (9.53) en la Ec. (9.50), es decir,
sustituyendo este vector Y y sus derivadas Y y Y , expresadas en función de las
formas modales X i , las ecuaciones de movimiento que obtendríamos serían:

INGENIERÍA SISMORRESISTENTE
SECC. 9.5: ANÁLISIS MODAL PARA EXCITACIONES SÍSMICAS 17
 n   n   n 
M   ai (t) X i  + C   a i (t) X i   K   ai (t) X i  =  M I uG ( t )
 i 1   i 1   i 1 
n
 [M X i ai (t) + C X i a i (t) + K X i ai (t) ] =  M I uG ( t )
i=1

Luego al premultiplicar esta última ecuación por el vector X Tj ( para

j  1,2,, n ), el cual es independiente de “ i ” , obtendríamos:

n
X Tj  [M X i ai ( t ) + C X i a i ( t ) + K X i a i ( t ) ] =  X Tj M I uG ( t )
i=1
n (9.54)
 [X Tj M X i ai ( t ) + X Tj C X i a i ( t ) + X Tj K X i a i ( t ) ] =  X Tj M I uG ( t )
i=1

Al aplicar las condiciones de ortogonalidad:

X Tj M X i = 0 para j  i pero si j i X iT M X i = 1 (9.55)

X Tj CX i = 0 para j  i pero si j i X iT C X i = 2  i  i (9.56)

( Si C tiene una forma especial )

X Tj K X i = 0 para j  i pero X iT K X i = ωi 2 (9.57)

en la Ec. (9.54), para “ j = i ” , ahora, esta sería:

X iT M X i ai (t) + X iT C X i a i (t) + X iT K X i ai (t) =  X iT M I uG ( t ) (9.58)

De manera similar a lo hecho en secciones anteriores, al escribir la Ec. (9.58) de

otra manera se tiene:

M ie ai (t) + Cie a i (t) + K ie ai (t) =  Fi e uG (t ) (9.59)

siendo M ie , Cie , K ie y Fi e escalares, correspondientes a cada modo de vibración

“ i ”. Luego, dividiendo la Ec. (9.59) entre M ie resulta :

C ie K ie Fi e
ai (t) + a i (t) + ai (t) =  uG (t ) (9.60)
M ie M ie M ie

Dr. JAVIER PIQUÉ DEL POZO

18 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

Siendo la Ec. (9.60) la que representa las “ n ” ecuaciones modales del

movimiento. Al observarla, vemos que posible realizar una analogía de ésta con el
caso cuando solo se tenía 1 GDL, o sea:

e
C C ie c c
 i (%)  i
 es equivalent e a  (%)  
C e
i ( crítico ) 2M  i
i
e
c crítico 2m

(9.61)
Ce X T C Xi X iT C X i
 Ti  2 i i  2 i i 
 Me Xi M Xi X iT M X i

Se demostró también en la Secc. 9.4.1 que:

X iT K X i
 i2  (9.62)
X iT M X i

y que además, según la Secc. 9.4.2 , el factor de participación estática “  i ”,

término que relaciona los modos y las matrices F y M, estaba dado por:

X iT F
i  solo que en este caso F = M I
X iT M X i

quedando entonces, expresado como:

n

XTM I
m j 1
j x ji
i  Ti  (9.63)
m j x ji 
n


Xi M Xi 2

j 1

Teniendo en cuenta que las Ecs. (9.58), (9.59), (9.60) y (9.61) podrían reducirse
debido a que los modos fueron normalizados respecto a la matriz de masas, o sea
X iT M X i  1 , según esto, dichas ecuaciones se escribirán así:

ai (t) + X i C X i ai (t) + X i K X i ai (t) =  X i MI uG (t ) .

T T T
(9.64)

2 i i  X iT C X i (9.65)

 i2  X iT K X i (9.66)

 i  X iT M I (9.67)

Entonces en la Ec. (9.64) al reemplazar las Ecs. (9.65), (9.66) y (9.67), se tiene:
INGENIERÍA SISMORRESISTENTE
SECC. 9.5: ANÁLISIS MODAL PARA EXCITACIONES SÍSMICAS 19

ai (t) + 2 i i a i (t)   i ai (t) =  i uG (t )

2
(9.68)

que representa las “ n ” ecuaciones modales de movimiento para un sistema

forzado considerando amortiguamiento.
“ Y ” se suele expresar, al igual que en las secciones anteriores , como:

n
Y   d i (t )  i X i (9.69)
i=1

Al relacionar la Ec. (9.69) con la Ec. (9.51), vemos que a i ( t )  d i ( t ) i . Donde

“ i (t ) ” , es el factor de participación dinámica (dependiente del tiempo) y “  i ” el
d
factor de participación estática (independiente del tiempo). Luego la Ec. (9.68) en
función de “ d i (t ) ”quedaría expresada como:

di (t ) + 2 i i d i (t )   i d i (t ) =  uG (t )
2
(9.70)
En resumen la respuesta estará dada por:
n
Y =  d i (t)  i X i (9.71)
i=1

di (t) + 2  i i di (t) + i d i (t) = -uG (t) para i = 1,2,...,n (9.72)
2

 i = X Ti MI ó X i MJ
T
(9.73)
1
(Si el modelaje es tridimensional )
Hay dos formas de realizar el análisis modal:
a) Se puede resolver cada ecuación modal tanto en el dominio del tiempo como
en el de frecuencias es decir integrada directamente o haciendo un cambio de
variables de t a  y resuelta en ese campo mediante el uso de las
transformadas de Fourier.
Es más usual lo primero en que la solución de la ecuación modal o sea toda
la historia en el tiempo de d i (t) es almacenada. Luego los modos se
superponen apropiadamente en cada intervalo de tiempo y el tiempo-historia
para cada efecto se revisa para encontrar su máximo valor. Esta superposición
tiene que ser repetida independientemente para cada efecto ya que los
coeficientes que afectan las respuestas modales (o sea las contribuciones de
cada modo a cada respuesta en particular) variarán de un efecto al otro. Por
ejemplo los desplazamientos de un piso relativo al terreno, la aceleración
absoluta de una masa o la fuerza cortante en una columna. Por consiguiente

1
Esto se verá mas adelante en la Secc. 9.7.1.

Dr. JAVIER PIQUÉ DEL POZO

20 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

este proceso es tedioso si se desea determinar muchas respuestas, como por

ejemplo todas las fuerzas en los elementos de un edificio. Nótese, sin
embargo, que la determinación de todas las fuerzas, momentos, cortes en una
estructura es un problema estático una vez que se aplica a cada pórtico un
juego de desplazamientos iguales a la forma modal X i . Estos valores
modales después serán multiplicados por d i (t) y  i .

b) El análisis modal puede también llevarse a cabo manteniendo para cada

modo sólo la máxima respuesta d i,máx . Esto es particularmente conveniente
cuando se usa un espectro de respuesta para representar el movimiento, en
vez de un registro -que es precisamente el caso de los análisis sísmicos
especificados en los códigos de diseño- ya que el valor d i,máx se lee
directamente del espectro para el amortiguamiento deseado,
d i,máx = S d (  i ,  i ) ; véase Cap. 5. Este procedimiento es el que se conoce
precisamente como análisis modal espectral.
Para ilustrar el primer procedimiento supongamos que el edificio de la Fig. 9.11
está siendo sometido a una aceleración de la base de -1.0 m/s² que actúa durante
medio segundo.

Para simplificar el análisis supongamos que no hay amortiguamiento. La solución

de la Ec. (9.72) para el caso sin amortiguamiento y en que üG = -1 está dada por:
1
d i (t) = (1 - cos i t) para t  t d (9.74)
i 2
1
di (t)= [ cosi (t - t d ) - cosit] para t > t d (9.75)
i
2

Luego el desplazamiento para cada modo está dado por las siguientes expresiones,
aplicando la Ec. (9.71):

0.017000 (1  cos 9.074 t )

 0.000600 (1  cos 33.05 t ) para t  0.5 (9.76)
0.000017 (1  cos 65.84 t )

Cuando t > 0.5 las expresiones se modifican de acuerdo a la Ec. (9.75).

La respuesta del desplazamiento del piso superior debido a los dos primeros modos
puede observarse en la Fig. 9.5. El modo 3 no tiene significación práctica.

INGENIERÍA SISMORRESISTENTE
SECC. 9.6: ANÁLISIS DINÁMICO MODAL ESPECTRAL 21

Fig. 9.5 Desplazamiento del piso superior. 1er y 2do modo

Pulso de -1m/s2 y td = 0.50 s en la base

9.6 ANÁLISIS DINÁMICO MODAL ESPECTRAL

En este caso particular del análisis modal la respuesta correspondiente al modo i estará
expresada como sigue:

Y i = S di  i X i (9.77)

donde S di es el valor leído del espectro de respuesta que se está usando y que
puede ser el valor máximo de la solución de la ecuación modal:

di (t) + 2  i  i d i (t) +  i d i (t) = - uG (t) con i = 1,2,..., n

2
(9.78)

ó el valor leído de un espectro teórico suavizado como los que se consignan en las
normas de diseño. Lo cierto es que en ambos casos del espectro se obtienen los
valores máximos de la aceleración, desplazamiento o velocidad para una frecuencia
determinada y un amortiguamiento fijo que son el dato de entrada para la expresión
(9.77)
El factor de participación estática tiene la expresión presentada anteriormente para
el caso de una excitación sísmica:

 i = X Ti MI ó X Ti MJ (9.79)
(si el modelaje es tridimensional)

Dr. JAVIER PIQUÉ DEL POZO

22 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

Nótese que los máximos para otros efectos como fuerzas en los elementos se
determinan para cada modo de un análisis estático, obteniendo los valores del juego
de desplazamientos Xi y multiplicándolos por S d y  i .

9.6.1 Combinación Modal

En el análisis modal espectral la determinación del efecto debido a la
superposición de todos los modos sólo puede ser hecha de forma aproximada
combinando (ya no superponiendo) las respuestas o participaciones modales. Como es
poco probable que todas las respuestas máximas de los modos coincidan en el tiempo,
sumar los valores absolutos de los valores modales máximos sería demasiado
conservador. El procedimiento establece que se deben calcular los efectos modales
para la respuesta que se desee: desplazamientos, fuerzas globales, efectos locales en
los elementos, y combinarlos siguiendo diversos criterios. Tradicionalmente se
calculaba la raíz cuadrada de la suma de los cuadrados de los máximos efectos (RCSC)
pero modernamente se están usando otras aproximaciones, cada una tratando de
acercarse a la respuesta predicha por el análisis tiempo-historia.
El Reglamento Nacional de Construcciones a través de su Norma de Diseño Sismo
Resistente [Ref. 10] prescribe (ordena) para el caso en que se use análisis dinámico
modal espectral que los modos se combinen usando el promedio ponderado de la raíz
cuadrada de la suma de las respuestas al cuadrado (RCSC) con la suma de los valores
absolutos (Σ ABS).
0.25  ABS + 0.75 RCSC (9.80)

Tradicionalmente se había usado sólo la RCSC pero se ha demostrado que es

insegura para edificios de más de 8 pisos [Ref. 7]. La de la Norma Peruana de
1997 [ Ref. 8 ] sin embargo, se ha evaluado que es conservadora. Por otro lado, en la
Norma E-030 de 1997: Diseño Sismorresistente [Ref. 10], se ha adoptado:

0.25  ABS + 0.75 RCSC (9.83)

También se está usando en otros países la llamada "Combinación Cuadrática

Completa", (CQC) (del Inglés: Complete Quadratic Combination) que es más
laboriosa de implementar pero que según sus promotores es más precisa que todas las
conocidas [Ref. 9]. Esta expresada de la siguiente forma:
Rk =   Rki  ij Rkj (9.84)

donde R representa las respuestas modales, desplazamientos o fuerzas. Y los

coeficientes de correlación están dados por:
8 β 2 ( 1+r) r 3 / 2 ωj
ρij = r= (9.85)
( 1+ r 2 )2 + 4 β 2 r( 1+r )2 ωi

INGENIERÍA SISMORRESISTENTE
SECC 6.6.1: COMBINACIÓN MODAL 23
Para concluir quizás debería recordarse que como el análisis modal espectral
involucra claramente una aproximación en la combinación de los efectos modales, el
grado de precisión que se tiene en el uso del espectro suavizado o de un registro
sísmico en particular no parece justificar la necesidad de mayor precisión.
En forma sintetizada, un esquema que contiene los pasos para realizar el análisis
dinámico modal espectral estaría dado por:

Análisis Dinámico Modal Espectral

Modelación de la Estructura
Definición de las matrices de masas y rigidez

Solución del Problema de Valores Característicos

Determinación de las frecuencias y periodos, modos

Cálculo de los Factores de Participación Estática

Leer Espectros de Diseño:

Aceleraciones o Desplazamientos

Cálculo de las Respuestas Modales Máximas

Combinación de las respuestas Modales para la

determinación de cada efecto
Fuerzas y deformaciones

Dr. JAVIER PIQUÉ DEL POZO

24 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

9.7 ANÁLISIS DINÁMICO SEUDO TRIDIMENSIONAL

Consideremos la edificación aporticada de 4 niveles mostrada en la Fig. 9.6. En
dicha figura se puede observar los elementos diafragmas rígidos, que juegan un papel
de suma importancia, puesto que distribuyen la fuerza horizontal, producida por el
movimiento en la base (excitación sísmica), sobre los elementos verticales. La
magnitud de la fracción de dicha fuerza sísmica que será tomada por cada elemento
vertical es función directa de sus rigideces, o sea, mientras mayor sea la rigidez del
elemento vertical tomará mayor fracción de la fuerza sísmica. Como ejemplo se podría
citar a los muros de corte, los cuales absorben, por no decir toda, gran parte de la
fuerza sísmica. Ello debido a su gran rigidez (esto se vio en el capítulo de rigideces).

Diafragmas Rígidos

uG, y (t )

z y
uG, x (t )
x
Fig. 9.6 Vista tridimensional de una edificación de 4 niveles . Se indican los
diafragmas rígidos y la aceleración de la base en ambas direcciones.

La vista de planta del piso “ p ”de la Fig. 9.6 se muestra en la Fig. 9.7. En dicha
vista no sólo se aprecia 1 GDL como en el caso del análisis en el plano (como se
vio el Cap. 8) sino mas bien se observa 3 GDL que corresponden a los
desplazamientos “ up ” y “ vp ”, y al giro del diafragma “  p ” (con respecto a
un eje perpendicular al plano que lo contiene). Además, el centro de masas
(C.M. ) debido a la excentricidad accidental “ eacc ” se mueve en ambas sentidos tal
como se muestra en la Fig. 9.7. Por ello, para obtener los mayores valores se
deberá resolver para dichos sentidos y escoger aquel que produzca los
mayores efectos.
También se debe acotar que la excentricidad no solo se da en la dirección “ x ”
(en el cual se moverá en los dos sentidos ya indicados) sino también en la dirección

INGENIERÍA SISMORRESISTENTE
SECC. 9.7: ANÁLISIS DINÁMICO SEUDO TRIDIMENSIONAL 25
del eje “ y ” (en él que, de manera análoga, se moverá en ambos sentidos). Lo dicho
significa que realmente de debe de hacerse 4 análisis: 2 en la dirección “ x ” y 2 en
la dirección “ y ”.

vp
p uG, y (t )
up
C.M

eacc eacc
y
uG, x (t )

x
Fig. 9.7 Vista de planta del piso “ p ” de la edificación de 4 niveles , donde se
muestra los sentidos de la excentricidad “ eacc ” en dirección del eje
“ x ” (de manera similar se da en la dirección y).

La irregularidad torsional es el principal problema debido a que los giros en

planta tienen resultados adversos. Por esto debe tratarse que el giro en planta tienda
a cero. Magnitudes de giros tales como 10 -2 rad son malas [Ref. 13].
Para realizar el análisis de un modelo seudo tridimensional se supone a la
estructura como un ensamble de pórticos planos, los cuales se encuentran
interconectados por un diafragma rígido. Lo que importa es el desplazamiento
horizontal a lo largo del alineamiento del pórtico (no el perpendicular a su
plano). Dicho desplazamiento puede describirse en función en de las tres
componentes de desplazamiento ( u o , vo y  o ) que definen el movimiento del
diafragma [ Ref. 11]. Para explicar tal relación entre dichos desplazamientos nos
basaremos en la Fig. 9.8:

vo
o
xo , y o  uo

y  xi , y i  ui
i
x
Fig. 9.8.a Planta genérica

Dr. JAVIER PIQUÉ DEL POZO

26 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

vo
alineamien to

xo , y o  uo
o
y ri
 xi , y i  ui
i

x
xo , y o 

ui ' ' o
T ri
ui ' vo
i i ui' ' '

uo
 xi , y i 
i n
(1) (2) (3)

Fig. 9.8.b Relación entre las coordenadas que describen el desplazamiento del
diafragma y el desplazamiento ( ui ) del pórtico “i”.

Debido a que no importa el desplazamiento perpendicular a su plano sino más

bien el desplazamiento a lo largo de su alineamiento (ui) de la Fig. 9.8.b se tiene
que:
u i  u i ' u i ' ' u i ' ' '

ui  uo cosi  vo seni  o ri (9.86)

La Ec. (9.86) escrita vectorialmente es:

u o 
 
ui  cos  i , sen i , ri vo   G i .u o (9.87)
 
 o
De la Fig. 9.8.b.(3) usando relaciones vectoriales se tiene que:

ri  T .n  x i  xo ,y i  yo 
. sen i , cos  i 
(9.88)
ri  x i  xo  sen i   y i  yo cos  i

La fuerza “Vi” producida en el pórtico, función de su rigidez lateral (kLi) y de su

desplazamiento (ui), también puede ser calculada en función de las coordenadas
que definen el desplazamiento del diafragma:
Vi  k Li ui (9.89)
usando la Ec. (9.87):

INGENIERÍA SISMORRESISTENTE
SECC. 9.7: ANÁLISIS DINÁMICO SEUDO TRIDIMENSIONAL 27
Vi  k Li G i .u o (9.90)

Para una explicación más detallada véase Ref. 11 y 13.

Luego de haber visto que es posible realizar el análisis considerando solo 3

GDL, los cuales definen el desplazamiento del diafragma, proseguiremos a definir
el momento polar de inercia “J”, el cual representa una medida de inercia
rotacional. Para efecto del análisis se descompondrá en la suma de momentos
polares de inercia. Veamos primeramente la expresión general para una placa de
masa “ Mp ” situada en e plano xy (ver Fig. 9.9), luego veremos lo concerniente a
nuestro caso.

J Mp  m p ( I x  I y )
J Mp  m p ( J )
r
dM p M
y J Mp  (J )
A

x
Fig. 9.9 Esquema para el hallar la expresión del Momento Polar de Inercia de
una Placa de masa “Mp” situada en el plano xy.

Dicha expresión general correspondiente a la Fig. 9.9 es:

J Mp   r 2 dM p

Para nuestro caso, como se dijo, se obtendrá el momento polar de inercia sumando
aquellos según como se muestra en la Fig. 9.10:
y

dM p
x
r  

J Mp   r 2 dM p  Jx  Jy

Fig. 9.10 Esquema para el hallar la expresión del momento polar de inercia de
una placa de masa “Mp” situada en el plano xy.

Dr. JAVIER PIQUÉ DEL POZO

28 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

En el análisis dinámico seudo tridimensional la dirección del sismo se toma a

través del factor de participación. O sea:

Para cada i , x  Sismo actuando en " x"

modo i , y  Sismo actuando en " y"
i ,

9.7.1 Análisis Modal Seudo Tridimensional

Para este análisis los vectores desplazamientos y aceleraciones son:

Yi  S di  i X i
(9.91)
U  S  X
i ai i i

Sabemos por lo visto en las secciones iniciales del presente capítulo que para un
análisis plano el factor de participación es:

XTM I
i  i

X TM Xi
i

Para el análisis seudotridimensional se tiene:

X TM Jx XTM Jy
 i ,x  i
y i , y  i
(9.92)
X T M Xi
i
X T M Xi
i

siendo definida la matriz de masas M como:

M x 0 0
M   0 My 0
 0 0 J o 

m1 0 0
..
0 m2 0 
..
y Mx  My  
: : :. : 
 
0 0 0 mn 

además Jo es la matriz de momentos polares de inercia de la masa.

INGENIERÍA SISMORRESISTENTE
SECC. 9.8: EJEMPLO DE APLICACIÓN 29
Al trabajar por separado en cada dirección se tiene para:

fuerzas solo en “ x ” : M Y  C Y  K Y   M J x uG , x (9.93)

fuerzas solo en “ y ” : M Y  C Y  K Y   M J y u G , y (9.94)

 u1  1 0
: : :  “n“
      componentes
u n  1 0
     
 v1  0 1
      “n“
donde: Y   : , J x  :  y J y  :  componentes
“3n“
v  0 1
 n    
 1  0 0
      “n“
: : :  componentes
 n  0 0

9.8 EJEMPLO DE APLICACIÓN

El edificio cuya planta y elevaciones se muestran en la Fig. 9.11 ha sido modelado
considerando tres grados de libertad por piso: dos desplazamientos horizontales u ,
v , y un giro en planta  . Esto da origen a 12 grados de libertad dinámicos: 4 pisos x 3
GDL por piso.
Usando el programa "A3s" 1, versión 4 (1991) desarrollado por el Dr. Hugo
Scaletti de la Universidad Nacional de Ingeniería se ha efectuado un análisis dinámico
seudo-tridimensional, modelando el edificio a base de pórticos y muros o placas.
La solución del problema de valores propios o característicos da como resultado
las frecuencias (períodos ya ordenados de mayor a menor), formas de modo y factores
de participación. A continuación se presentan los valores numéricos para los 8
primeros modos, de un total de 12. Los modos son vectores con tres componentes por
piso o nivel, cada una correspondiendo a los grados de libertad dinámicos. Su
dimensión real es (12x1) en este caso ó (3nx1). La primera columna (u) corresponde a
sus componentes en la dirección X , la segunda columna (v) a las componentes en la
dirección Y y la tercera a las componentes de giro (  ). Las filas corresponden a cada
piso. Las formas de modo están normalizadas, es decir X Ti M X i = 1 .

1
Este programa usa como separador decimal el punto y no la coma.

Dr. JAVIER PIQUÉ DEL POZO

30 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

Fig. 9.11 Edificio de 4 pisos. Ejemplo de aplicación.

(a) Planta

INGENIERÍA SISMORRESISTENTE
SECC. 9.8: EJEMPLO DE APLICACIÓN 31

Fig. 9.11 Edificio de 4 pisos. Ejemplo de aplicación

(b) Elevaciones

Dr. JAVIER PIQUÉ DEL POZO

32 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

MODO 1 MODO 2
T = .925 seg f = 1.4441 Hertz, T = .1901 seg f=5.2606 Hertz
 = 9.074 rad/seg  = 33.053 rad/seg
Nivel u v  u v 
4 .203504 -.000071 -.000023 .238364 .000998 .000319
3 .163405 -.000070 -.000028 .021132 .000981 .000294
2 .107715 -.000059 -.000023 -.146275 .000940 .000297
1 .045146 -.000040 -.000015 -.125810 .000697 .000230
Factores de Participación:
6.959185 -.003430 -.043786 -2.639089 .052041 .546903

MODO 3 MODO 4
T = .0954 seg f = 10.4788 Hertz T = .0945 seg f=10.5837 Hertz
 = 65.840 rad/seg  = 66.500 rad/seg

Nivel u v  u v 
4 -.101434 .140142 .007967 .187837 .079009 .004222
3 .054347 .130545 .007733 -.107718 .073480 .004060
2 .013816 .102608 .006220 -.009803 .056968 .002995
1 -.086816 .058343 .003793 .148186 .031241 .001410
Factores de Participación:
-.679863 6.050405 12.184690 1.232911 3.368265 5.948236

MODO 5 MODO 6
T = .0670 seg f = 14.9365 Hertz T = .0581 seg f=17.1984 Hertz
 = 93.849 rad/seg  = 108.061 rad/seg
Nivel u v  u v 
4 -.109810 .017752 -.004659 -.019704 -.058716 .027787
3 .110824 .015356 -.004346 .018992 -.048260 .025470
2 -.142864 .011724 -.003359 -.025454 -.036213 .019961
1 .138948 .004663 -.002505 .029646 -.018933 .011124
Factores de Participación:
.764115 .680903 -6.989503 .193751 -2.210295 39.458440

INGENIERÍA SISMORRESISTENTE
SECC. 9.8: EJEMPLO DE APLICACIÓN 33

MODO 7 MODO 8
T = .0319 seg f = 31.3972 Hertz T = .0216 seg f= 46.3189 Hertz
 = 197.274 rad/seg  = 291.030 rad/seg
Nivel u v  u v 
4 -.000582 .187761 .007484 .000295 .236363 -.016025
3 .000247 .074525 .005383 .000050 -.062996 -.006327
2 .000776 -.118627 -.005540 -.1114 -.080886 .007818
1 -.002348 -.149932 -.007965 .002601 .101583 .018901
Factores de Participación:
-.019667 -1.828819 -2.461876 .020819 .368742 6.904037

Las masas y el momento polar de inercia de las mismas son las siguientes: (t - s 2 /m)

Nivel xo yo Masa (x ó y) Jo
4 8.50 4.10 6.40E+00 2.16E+02
3 8.50 4.87 1.84E+01 6.14E+02
2 8.50 4.87 1.84E+01 6.14E+02
1 8.50 6.50 1.48E+01 5.00E+02

donde x o e y o corresponden a las coordenadas del centro de masas. (m)

Con estos valores se ha formado la matriz (M) de masas, que es una matriz
diagonal de 12x12. Los cuatro primeros términos corresponden a la masa de cada piso
en la dirección X , o sea M x , los siguientes 4 son las mismas masas que
corresponden a la dirección Y , o sea M y y los últimos 4 son los momentos polares
de inercia de la masa, J 0 .

M x 0 0
M   0 M y 0
 0 0 J o 
Los factores de participación  i se han calculado aplicando la expresión (9.16)

 i = X Ti MJ (9.16)

donde J es un vector con unos y ceros dependiendo de donde proviene el sismo que
se está considerando. Recuérdese que al momento de desacoplar las ecuaciones de

Dr. JAVIER PIQUÉ DEL POZO

34 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

movimiento mediante la descomposición modal, éstas (las ecuaciones de movimiento)

se plantean por separado para cada dirección de la aceleración de la base. Por
consiguiente si el sismo es en la dirección X habrá que colocar 1s en los primeros
n componentes de J y el resto cero. Si el sismo es en la dirección Y , serán 1s desde
n + 1 hasta 2n y el resto cero. Veamos:
 1 0 0
     
J x  0, J y   1 y J   0
0 0  1
     
Por ejemplo  1x , el factor de participación estática para el primer modo cuando el
sismo actúa en la dirección X es:

 1x = X T1 M J x (9.95)

T = [(1111) (0 0 0 0) (0 0 0 0)] (9.96)

J x

Al efectuar el producto matricial debido a que sólo los cuatro primeros términos de
J son 1s, el resto no contribuye, por lo que (9.25) equivale a multiplicar solamente
los cuatro primeros términos de la matriz de masa por las componentes de la forma
modal correspondientes a la dirección X .
4
Γ 1x =  u i mi
i 1

Γ1x = 6.4 x.203504+18.4 x.163405+18.4 x.107715+14.8x.045146 = 6.95919

9.8.1 Cálculo de Respuestas Modales

Para proseguir con el análisis sísmico usando el análisis dinámico modal espectral, es
necesario considerar el sismo mediante un espectro de diseño. En este caso se ha
usado el espectro de las Normas Peruanas (8) RNC, con un factor de reducción por
ductilidad Rd = 3. Este espectro considera un porcentaje de amortiguamiento (  ) del
5%.

9.8.2 Desplazamientos
Los desplazamientos correspondientes a cada modo se obtienen aplicando la Ec. (9.77)
Y i = S di  i X i (9.77)

y luego combinando estas contribuciones usando la combinación estipulada en el

RNC.

INGENIERÍA SISMORRESISTENTE
SECC. 9.8.2: DESPLAZAMIENTOS 35
Por ejemplo, para determinar el desplazamiento en el piso superior cuando el
sismo actúa en la dirección X se aplicará la expresión (9.77) para cada modo.
Los valores espectrales del desplazamiento S d se leen del espectro de
desplazamientos, o como en este caso en que se tiene como dato el espectro de
aceleraciones S a , se calcula S d usando la relación que existe entre ellos mediante
S d = S a /  (véase Cap. 5).
2

Las valores para cada modo son:

Modo  Sa Sd i Xi yi
2
rad/s m/ s m m
1 9.074 8.060E-01 9.790E-03 6.959185 .203504 .013865
2 33.053 1.330E+00 1.217E-03 -2.639089 .238364 -.000766
3 65.840 1.330E+00 3.068E-04 -.679863 -.101434 .000021
4 66.500 1.330E+00 3.008E-04 1.232911 .187837 .000070
5 93.849 1.330E+00 1.510E-04 .764115 -.109810 .000041
6 108.061 1.330E+00 1.139E-04 .193751 -.019704 .000000
7 197.274 1.330E+00 3.418E-05 -.019667 -.000582 .000000
8 291.030 1.330E+00 1.570E-05 .020819 .000295 .000000
 ABS = .014763
RCSC = .013886

Aplicando la combinación del RNC, o sea el promedio de ambos valores se tiene:

1/ 2 (  ABS + RCSC) = 0.0143 m (9.97)

De manera análoga se obtiene, para los restantes pisos (niveles), los

correspondientes desplazamientos para el sismo en la dirección X (en metros):
Nivel x y r
4 1.431E-02 5.598E-05 5.033E-06
3 1.120E-02 5.216E-05 5.077E-06
2 7.593E-03 4.118E-05 4.166E-06
1 3.336E-03 2.351E-05 2.648E-06
Como puede observarse, los desplazamientos en la dirección Y debidos al sismo
en la dirección X son muy pequeños, igualmente los giros en planta. Indicando que
hay poca influencia de la torsión. Los desplazamientos reales, de acuerdo a la Norma,
serán los calculados en el análisis anterior multiplicados por 0.75 Rd , o sea por 2.25 .

Dr. JAVIER PIQUÉ DEL POZO

36 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

En este caso el máximo desplazamiento del piso superior será de

2.25 x 1.43 cm = 3.22 cm .
9.8.3 Fuerzas
Se pueden determinar las fuerzas globales, como cortes y fuerzas aplicadas en cada
piso, para todo el edificio o para cada pórtico, o también los efectos locales, o sea
momentos, cortes, fuerzas axiales en cada viga y columna. En cada caso, como se ha
visto previamente en la teoría, es necesario efectuar la combinación de las
contribuciones modales para cada efecto por separado. Es teóricamente incorrecto
obtener efectos modales de cortantes, calculándolos a partir de las fuerzas aplicadas ya
combinadas. Los resultados son muy distintos, como puede apreciarse de los valores
que se presentan a continuación.
Efectos Globales - Sismo actuando según la Dirección X
Fuerzas Concentradas (en toneladas)
Nivel x y r
4 1.264E+01 1.561E+00 4.443E+00
3 2.096E+01 4.085E+00 1.187E+01
2 1.936E+01 3.211E+00 9.324E+00
1 1.297E+01 1.458E+00 4.627E+00
  65.93 t
Cortantes en Cada Nivel (en toneladas)
Nivel x y r
4 1.264E+01 1.561E+00 4.443E+00
3 2.953E+01 5.613E+00 1.631E+01
2 3.849E+01 8.757E+00 2.541E+01
1 4.597E+01 1.017E+01 2.990E+01

Nótese primeramente que la distribución de las fuerzas en altura no es triangular.

Asimismo se puede observar que el cortante en la base calculado superponiendo
directamente los cortantes que se obtienen en cada modo (como debe ser), son
menores que si se calcularan sumando las fuerzas resultantes en cada piso, como se
haría en el análisis estático. La fuerza cortante en la base es de 45.97 toneladas,
mientras que si se suman las fuerzas de cada piso se obtiene 65.93 t. Recuerde que la
combinación es la última operación que se realiza para obtener cualquier efecto.
Las fuerzas pueden obtenerse de dos maneras:
a) Determinando las aceleraciones modales para cada modo y multiplicando por
las masas o,
b) Determinando los desplazamientos de cada modo y multiplicando por la
matriz de rigidez lateral. En cada caso se usará la que corresponda, la de todo
el edificio si se desean valores globales, o la de cada pórtico -con los
desplazamientos de cada pórtico- si se desean los efectos por pórtico. Para

INGENIERÍA SISMORRESISTENTE
SECC. 9.8.3: FUERZAS 37
determinar los cortantes habrá que calcularlos para cada modo, con
cualquiera de los procedimientos mencionados. Recuérdese que la solución
del pórtico o del edificio para cada modo es un problema estático, de manera
que aplicando al pórtico los desplazamientos de un modo se pueden
determinar todos los efectos, tanto globales como locales y luego combinar la
contribución de cada modo para cada efecto por separado.
A continuación se ilustra el primer procedimiento; es decir, determinando las
fuerzas en función de las aceleraciones para cada modo. Las aceleraciones modales se
obtienen aplicando la siguiente expresión:
Yi = S ai  i X i (9.98)
y luego combinando estas contribuciones usando el procedimiento estipulado en el
RNC (8).
Por ejemplo, para determinar la fuerza global que se presenta en el piso superior de
todo el edificio cuando el sismo actúa en la dirección X se aplicará la expresión
(9.28) para determinar la contribución de la aceleración en cada modo y luego la
fuerza correspondiente.
Los valores espectrales de la aceleración Sa se leen directamente del espectro de
aceleraciones Sa .
Los valores para cada modo son:

Modo  Sa Sd i Xi yi
2
rad/s m/ s m m
1 9.074 .806 6.959185 .203504 1.141475 7.305439
2 33.053 1.330 -2.639089 .238364 -.836655 -5.354591
3 65.840 1.330 -.679863 -.101434 .091718 .586998
4 66.500 1.330 1.232911 .187837 .308010 1.971263
5 93.849 1.330 .764115 -.109810 -.111597 -.714220
6 108.061 1.330 .193751 -.019704 -.005078 -.032496
7 197.274 1.330 -.019667 -.000582 .000015 .000097
8 291.030 1.330 .020819 .000295 .000008 .000052
 ABS = 15.9650
RCSC = 9.3167

Aplicando la combinación de la Norma Peruana (8), o sea el promedio de ambos

valores se tiene

1/2(  ABS + RCSC)= 12.64toneladas (9.99)

Dr. JAVIER PIQUÉ DEL POZO

38 CAP. 9:ANÁLISIS SÍSMICO POR SUPERPOSICIÓN MODAL ESPECTRAL

Obsérvese que esta fuerza está constituida por contribuciones importantes del
primero, segundo y cuarto modos. La fuerza debida al 2° modo es el 70% del 1° Y la
del 4° modo es el 27% de la del 1o. Por lo tanto se aprecia que para la determinación
de fuerzas que actúan sobre la estructura el análisis dinámico es una herramienta más
adecuada considerando acciones que un análisis estático no puede representar.
En el caso de la determinación de desplazamientos sin embargo la contribución de
los modos superiores es prácticamente despreciable. En la primera parte de este
ejemplo para determinar el desplazamiento del piso superior, se puede observar que el
desplazamiento del 2o. modo solo representa el 5% del 1o. Por lo tanto bastaría
considerar los desplazamientos debidos únicamente al primero y ahorrarse la
combinación modal.

INGENIERÍA SISMORRESISTENTE
REFERENCIAS 39
REFERENCIAS
1. Biggs, J.M., Dynamic Analysis of One-Degree Systems, en Notas del
curso Fundamentals of Earthquake Engineering for Buildings .
Massachusetts Institute of Technology. Cambridge, Massachusetts. 1972
2. Röesset, J.M. Structural Dynamics. Notas de clase. Massachusetts
Institute of Technology. Cambridge, Mass. 1974.
3. Biggs, J.M., Introduction to Structural Dynamics. McGraw-Hill. New
York. 1964
4. Craig Jr., R.R., Structural Dynamics. John Wiley & Sons. New York.
1981
5. Clough, R.W. & Penzien, J. Dynamics of Structures. McGraw-Hill. New
York. 1975
6. Okamoto, S. Introduction to Earthquake Engineering. Halsted Press.
John Wiley & Sons. New York. 1973
7. Piqué, J., Echarry, A. A modal combination for Dynamic Analysis of
Reinforced Concrete Frames. 9a. Conferencia Mundial de Ingeniería
Antisísmica. Tokyo-Kyoto. Japón. 1988
8. O.I.N., Ministerio de Vivienda y Construcción. Normas Peruanas de
Diseño Sismo Resistente. Lima, 1977
9. Wilson, E.L., Der Kiureghian, A., Bayo, E.P. A Replacemente for the
SRSS Method in Seismic Analysis Journal of Earthquake Engineering
and Structural Dynamics, Vol. 9, 187'194. 1981
10. SENCICO, Ministerio Transportes Comunicaciones Vivienda y
Construcción. Norma Técnica de Edificación, NTE E-030 Diseño
Sismorresistente. Lima, 1997
11. Piqué, J., Scaletti, H., Análisis Sísmico de Edificios, Ediciones Capítulo
de Ingeniería Civil. Lima, Perú. 1991
12. Bazán, E., Meli, R. Diseño Sísmico de Edificios, Editorial Limusa.
Balderas, México. 2002
13. Scaletti, H., "Análisis Sísmico Seudo Tridimensional" en Notas del curso
de Análisis Estructural II. Universidad Nacional de Ingeniería. Lima,
Perú. 2003

Dr. JAVIER PIQUÉ DEL POZO