I.

INTRODUCCIÓN

Los diseños de ingeniería cada vez requieren más una adecuada respuesta
dinámica. Esto puede deberse a que las cargas realmente se apliquen de forma
muy rápida, como a asignar una mayor importancia a aspectos como el
mantenimiento de la funcionalidad, la resistencia, y el confort ante las vibraciones.
Estas condiciones de diseño a menudo se añaden a las puramente estáticas, de
estabilidad y resistencia en la posición de equilibrio.
El sismo o terremoto es producido por la liberación repentina de una acumulación
tensional, originada por los movimientos de placas de la corteza terrestre. Esta
“explosión” de tensiones se transmite de forma ONDULATORIA a través de la
superficie (como las olas en la superficie del mar). Es fundamental la comprensión
del sismo como fenómeno ondulatorio, y no como simple oscilación horizontal,
sobre todo en estructuras grandes donde no podemos despreciar la componente
vertical del movimiento sísmico (se produce una oscilación horizontal y otra
vertical)
Diremos que en estas estructuras despreciamos la interacción suelo-estructura, de
tremenda importancia en otros casos (puentes, presas, rascacielos...), en las que
hay que considerar que el suelo se deforma con el sismo. Sin embargo, antes de
enfrentarnos al fenómeno sísmico (cargas dinámicas) hemos de analizar como
vibran las estructuras, pues dependiendo de cómo lo hagan, responderán de una
forma u otra al sismo.
Recurriremos al análisis modal de la estructura previo a cualquier cálculo.
Posteriormente podremos ampliar el análisis al análisis espectral, donde ya
consideramos la respuesta de la estructura al sismo.
Para los análisis modal y espectral, el método matemático más extendido de
solución es la discretización por elementos finitos (para los análisis en los que se
considera la interacción suelo-estructura, el método más ajustado a la realidad es
el modelado por elementos de contorno).
El análisis modal es una técnica utilizada para determinar las características
vibratorias de una estructura con comportamiento elástico y lineal, nos suministra
las frecuencias y formas naturales con las que vibra una estructura. Es el tipo de
análisis dinámico más importante, y debe de hacerse siempre antes de cualquier
otro análisis, pues las características vibratorias de una estructura determinan
como responde ante las cargas dinámicas (fundamental para evitar la resonancia).
Si sometemos el sistema a una excitación, éste responderá vibrando de una
determinada manera, propia y exclusiva del sistema concreto. Esta respuesta se
podrá expresar en un gráfico desplazamiento - frecuencia, (puede no ser
armónico, pero que puede descomponerse por la transformada de Fourier en
varias funciones armónicas).
A esta vibración la llamaremos MODO de vibración del sistema, y a su frecuencia,
FRECUENCIA NATURAL del sistema. Cada modo tendrá una forma propia de
vibración Si realizamos esta misma operación con todos los posibles sistemas de
un grado de libertad obtendremos todos los modos de vibración y sus formas, así
como sus frecuencias naturales.
Un ESPECTRO DE RESPUESTA es una curva que representa la respuesta de un
sistema idealizado ante una excitación (la respuesta puede ser la fuerza, el
desplazamiento, la velocidad o la aceleración).Cojamos el acelerograma de un
terremoto y un sistema de un grado de libertad, con una masa y rigidez,
sometemos ese sistema al acelerograma, y de todas las aceleraciones que sufre
el sistema nos quedamos con la máxima. Cojamos otro sistema de 1 G.D.L
distinto, con otra masa, y rigidez, y lo volvemos a someter al acelerograma y nos
quedamos con la máxima respuesta de aceleración... si continuamos así, al final
obtenemos una curva de máximas aceleraciones de todos los sistemas de un
grado de libertad según sus frecuencias naturales.
 II. MARCO TEÓRICO

Análisis Modal Espectral

 Ecuaciones modales para sistemas no amortiguados

Ecuación de equilibrio dinámico para un sistema lineal no amortiguado

El vector de desplazamientos se puede expresar en términos de las contribuciones

modales

Al sustituir la ecuación anterior en la ecuación de equilibrio dinámico

Premultiplicando se obtiene

Por el principio de ortogonalidad, todos los términos en cada una de las

sumatorias desaparecen excepto en donde r=n, por lo que se obtiene

Donde:
 Análisis de la historia de respuesta

mu  g  cu  ku  0
  1u
 n  2 nn qn  n  nu
2
q g

nT m1
n 
Mn
Γn=coeficiente de participación del modo n

Contribución del modo n al desplazamiento del nodo u

Espectros de Diseño
 Los espectros de sismos reales tienen forma irregular y presentan variaciones
bruscas de la respuesta máxima en función del periodo natural. Por lo anterior,
para fines de diseño, los reglamentos de construcción indican espectros
suavizados en los que se ensanchan los picos y se eliminan los valles
Análisis Modal Espectral
 Estructuras sin torsión:
Los desplazamientos de las masas de los sistemas de varios grados de libertad
pueden expresarse como la suma de los desplazamientos debidos a la
participación de cada uno de los modos naturales. Para el desplazamiento de la j-
ésima masa:
N N
u j (t )   u jn (t )   n jn Dn (t )
n 1 n 1

Si la excitación sísmica se define por un acelerograma, la función Dn(t) es posible

calcularla con cualquiera de los métodos numéricos vistos anteriormente. Con la
expresión anterior es posible determinar los desplazamientos de la estructura en
cualquier instante; conocidos estos es posible el cálculo de otras cantidades de
interés para el diseño, como es el caso de las fuerzas cortantes o momentos
flexionantes.
El análisis espectral se basa en que interesa la respuesta máxima que producirá
un sismo futuro. Debido a que no es posible predecir con exactitud los
acelerogramas de eventos sísmicos venideros, los reglamentos normalmente
prescriben la intensidad sísmica de diseño mediante espectros suavizados que
suministran la seudoaceleración máxima An, para cada periodo Tn. An es igual al
desplazamiento máximo por la frecuencia ω al cuadrado; por lo tanto, el espectro
de diseño proporciona el valor máximo de Dn(t) dado por An/ωn2.

u jn máx  An n jn / n
2

Nos brinda la contribución máxima del modo n al desplazamiento de la masa j

Combinación de respuestas modales

La expresión anterior permite conocer cualquier respuesta (fuerza cortante,
deformación de entrepiso, momento de volteo, etc.) máxima rn de la estructura
debido al modo n. sin embargo, para fines de diseño interesa determinar la
respuesta máxima R de la estructura por la participación de todos los modos.
Debido a que las respuestas máximas de cada uno de los modos no coinciden
entre ellos, no se puede hacer una suma directa de los máximos y se han
propuesto distintos métodos para combinar los máximos de las respuestas
modales; uno de los métodos más sencillos es el que considera:
 R r
n
n
2

El cual es recomendable utilizar cuando los periodos naturales difieren al menos

en un 10% entre sí.
La regla de combinación de modos es la siguiente:
 Se calcula la respuesta de interés para cada modo
 Se combinan los resultados de cada modo mediante la raíz cuadrada de la
suma de los cuadrados de las contribuciones de cada modo.
EJERCICIO DE ANALISIS MODAL ESPECTRAL
Método de superposición modal espectral

Se tiene la siguiente estructura

Figura No. 1 Vista en planta de la estructura

Figura No. 2 vista frontal de la estructura.

Las vigas de la estructura tienen la siguiente forma
35 × 40 [𝑐𝑚] − − − − − − − − − 𝑒𝑗𝑒 𝐵
30 × 40 [𝑐𝑚] − − − − − − − − − 𝑒𝑗𝑒 𝐴 𝑦 𝐶
Edificio de 4 pisos con los siguientes parámetros de espectro
𝐴𝑎 = 0.25
𝐴𝑣 = 0.25
𝐹𝑎 = 1.15
𝐹𝑣 = 1.55
𝐼 = 1.0

Dadas los pesos por piso registrados en la figura anterior se tiene que la matriz de
masa [Ton] de la estructura es:

12.2324 0 0 0
0 12.2324 0 0
[𝑀] = ( ) (𝑡𝑜𝑛)
0 0 12.2324 0
0 0 0 7.1356

Se realizara un análisis de movimiento en sentido E-W, por lo tanto es importante

considerar que para este caso particular de estructura, las vigas del eje B son
distintas a las vigas de los ejes A y C
Con ayuda del programa de Matlab MRIG_KL, ofrecido por el profesor, fue posible
determinar las matrices laterales de rigidez de las ejes Ay C (la misma) y la del eje
B. Estas matrices están dadas en unidades de [T/m]
Eje B
5.6731 −4.3807 1.1788 −0.1630
−4.3807 7.2395 −4.6078 0.9686 𝑡𝑜𝑛
[𝐾𝐿 𝐵 ] = 1000 ∗ ( )( )
1.1788 −4.6078 6.6844 −3.1183 𝑚
−0.1630 0.9686 −3.1183 2.2934
Eje A-C
5.6097 −4.3507 1.2288 −0.1750
−4.3507 7.1009 −4.5653 0.9952 𝑡𝑜𝑛
[𝐾𝐿,𝐸𝐽𝐸 𝐴−𝐶 ] = 1000 ∗ ( )( )
1.2288 −4.5653 6.4822 −2.9956 𝑚
−0.1750 0.9952 −2.9956 2.1538

Es importante considerar una matriz lateral general para todo el edificio que no es
más que la suma de las dos matrices laterales obtenidas anteriormente,
considerando que la matriz del eje A y C debe multiplicarse por 2.
 𝐾 = 2 ∗ 𝐾𝐿 𝐴 𝑦 𝐶 + 𝐾𝐿𝐵
1.6892 −1.3082 0.3636 −0.0513
[𝐾𝑡𝑜𝑡𝑎𝑙 ] = 10000 (−1.3082 2.1441 −1.3738 0.2959
)
0.3636 −1.3738 1.9649 −0.9110
−0.0513 0.2959 −0.9110 0.6601
Al resolver el problema de valores propios de la ecuación de movimiento, con
ayuda del software MODOS EDIFICIO, y empleando la matriz de rigidez K general
de la estructura, se obtiene

Modos y periodos del edificio en la dirección este-oeste.

MODO 𝑇(𝑠)
1 0.9747
2 0.2946
3 0.1572
4 0.1052
Tabla No.1 Modos y Periodos de la estructura

La matriz [Φ]4x4 contiene los modos normalizados con respecto a la masa:

−0.0762 0.1741 0.1802 −0.1147
−0.1378 0.1476 −0.0835 0.1844
[Φ] = ( )
−0.1812 −0.0423 −0.1408 −0.1652
−0.2034 −0.2185 0.1963 0.1118
Espectro de aceleraciones.

Figura No. 1 Espectro de aceleraciones.

Como primera medida se realiza el cálculo de los límites de periodo en los que se
encuentra la estructura, es decir:
𝐴𝑣 𝐹𝑣 0.25 ∗ 1.55
𝑇0 = 0.1 ( ) = 0.1 ( ) = 𝟎. 𝟏𝟑𝟒𝟖 𝒔
𝐴𝑎 𝐹𝑎 0.25 ∗ 1.15
𝐴𝑣 𝐹𝑣 0.25 ∗ 1.55
𝑇𝑐 = 0.48 ( ) = 0.48 ∗ ( ) = 𝟎. 𝟔𝟒𝟕𝟎 𝒔
𝐴𝑎 𝐹𝑎 0.25 ∗ 1.15
𝑇𝐿 = 2.4𝐹𝑣 = 2.4 ∗ 1.55 = 𝟑. 𝟕𝟐𝟎𝟎 𝒔

Una vez calculados los periodos limite, se observa en la figura No.1 y se efectúa el
cálculo de las aceleraciones (𝑆𝑎 ).
Análisis Modal
Modo 1. 𝑇 = 0.9747 𝑠, 𝑇 > 𝑇0 > 𝑇𝑐 𝑦 𝑇 < 𝑇𝐿
𝐴𝑣 𝐹𝑣 𝐼 0.25 ∗ 1.55 ∗ 1
𝑆𝑎 = 1.2 ( ) = 1.2 ∗ ( ) = 0.4771
𝑇 0.9747
Modo 2. 𝑇 = 0.2946 𝑠, 𝑇0 < 𝑇 < 𝑇𝑐
𝑆𝑎 = 2.5𝐴𝑎 𝐹𝑎 𝐼 = 2.5 ∗ 0.25 ∗ 1.15 ∗ 1 = 0.7188
Modo 3. 𝑇 = 0.1572 𝑠, 𝑇0 < 𝑇 < 𝑇𝑐
𝑆𝑎 = 2.5𝐴𝑎 𝐹𝑎 𝐼 = 2.5 ∗ 0.25 ∗ 1.15 ∗ 1 = 0.7188
Modo 4. 𝑇 = 0.1052 𝑠, 𝑇 < 𝑇0
0.6𝑇 0.6∗0.1052
𝑆𝑎 = 2.5 ∗ 𝐴𝑎 𝐹𝑎 𝐼 (0.4 + ) = 2.5 ∗ 0.25 ∗ 1.15 ∗ 1 ∗ (0.4 + ) = 0.6240
𝑇0 0.1348

A continuación se presenta el cálculo del factor de participacion modal, {Γ}4𝑥1 =

[Φ]𝑇 4𝑥4 [𝑀]4𝑥4 {𝛾}4𝑥1

Donde {𝛾}4𝑥1 es un vector unitario de 4 filas

𝑇
−0.0762 0.1741 0.1802 −0.1147
{Γ} = (−0.1378 0.1476 −0.0835 0.1844
)
−0.1812 −0.0423 −0.1408 −0.1652
−0.2034 −0.2185 0.1963 0.1118
12.2324 0 0 0 1
0 12.2324 0 0 1
∗( )∗( )
0 0 12.2324 0 1
0 0 0 7.1356 1
−6.2849
{Γ} = ( 1.8582 )
0.8620
−0.3699
Ahora se procede a calcular los de los desplazamientos máximos a partir del
𝑇2
espectro de aceleraciones, 𝑍𝑖 𝑚𝑎𝑥 = |Γ𝑖 | 4𝜋𝑖 2 ∗ 𝑆𝑎 (𝑇𝑖 , 𝜉𝑖 )

0.97472
𝑍1 𝑚á𝑥 = 9.81 ∗ 6.2849 ∗ ∗ 0.4771 = 0.7078
4𝜋 2
0.29462
𝑍2 𝑚á𝑥 = 9.81 ∗ 1.8582 ∗ ∗ 0.7188 = 0.0288
4𝜋 2
0.15722
𝑍3 𝑚á𝑥 = 9.81 ∗ 0.8620 ∗ ∗ 0.7188 = 0.0038
4𝜋 2
0.10522
𝑍4 𝑚á𝑥 = 9.81 ∗ 0.3699 ∗ ∗ 0.6240 = 6.3420 ∗ 10−4
4𝜋 2
A continuación se presenta una tabla con el resumen de los pasos anteriores.
Modo 𝑇 (𝑠) 𝑆𝑎 [𝑔] Γ𝑖 𝑍𝑖 𝑚á𝑥
1 0.9747 0.4771 -6.2849 0.7078
2 0.2946 0.7188 1.8582 0.0288
3 0.1572 0.7188 0.8620 0.0038
4 0.1052 0.6240 -0.3699 6.3420 ∗ 10−4
Tabla No.2 Resumen de cálculo

Desplazamientos modales
Para calcular los desplazamientos modales {𝑢𝑖 }, se puede realizar de dos formas:
la primera trabajando con los vectores fila de la matriz de modos normalizados (Φ)
{Φi }|Γ𝑖 |𝑆𝑎 (𝑇𝑖 ,𝜉𝑖 )
es decir: {𝑢𝑖 } = donde el subindice 𝑖 indica el numero del modo. El
𝜔𝑖 2
otro medio es directamente de forma matricial, un proceso mucho más rápido, se
tiene [𝑈] = [Φ][Zm ] donde la matriz [Zm ] es una matriz diagonal con los valores
máximos 𝑍𝑖 𝑚á𝑥 . Asi se tiene que:
−0.0762 0.1741 0.1802 −0.1147
−0.1378 0.1476 −0.0835 0.1844
[𝑈] = ( )
−0.1812 −0.0423 −0.1408 −0.1652
−0.2034 −0.2185 0.1963 0.1118
0.7078 0 0 0
0 0.0288 0 0
∗( )
0 0 0.0038 0
0 0 0 6.3420 ∗ 10−4
−0.0539 0.0050 0.0007 −0.0001
[𝑈] = (−0.0975 0.0043 −0.0003 0.0001
)
−0.1283 −0.0012 −0.0005 −0.0001
−0.1440 −0.0063 0.0007 0.0001
Con estos resultados es posible obtener los máximos factibles de desplazamiento
para cada piso, con la siguiente expresión:

𝑢𝑗𝑚á𝑥 = √∑(𝑢𝑗𝑖 )2
𝑖

Donde el subíndice 𝑗 indica el piso.

𝑢4𝑚á𝑥 = √0.05391 + 0.00502 + 0.00072 + 0.00012 = 0.0541 𝑚

𝑢3𝑚á𝑥 = √0.09752 + 0.00432 + 0.00032 + 0.00012 = 0.0976 𝑚

𝑢2𝑚á𝑥 = √0.12832 + 0.00122 + 0.00052 + 0.00012 = 0.1283 𝑚

𝑢1𝑚á𝑥 = √0.14402 + 0.00632 + 0.00072 + 0.00012 = 0.1441 𝑚

Y el vector de los desplazamientos máximos factibles es:

0.0541
{𝑈𝑚á𝑥 } = (0.0976) (𝑚)
0.1283
0.1441

Cálculo de las derivas de piso

Con los resultados de desplazamientos por piso para cada modo, es posible
calcular las derivas modales del piso, 𝐷𝑗 = 𝑢𝑗 − 𝑢𝑗−1 y estos valores agruparlos en
una matriz de derivas modaldes [𝐷].
[𝐷]
−0.0539 − (−0.0974) 0.0050 − (0.0043) 0.0007 − (−0.0003) −0.0001 − (0.0001)
−0.0974 − (−0.1283) 0.0043 − (−0.0012) −0.0003 − (−0.0005) 0.0001 − (−0.0001)
=( )
−0.1283 − (−0.1440) −0.0012 − (−0.0063) −0.0005 − (0.0007) −0.0001 − (0.0001)
−0.1440 − (0) −0.0063 − (0) 0.0007 − (0) 0.0001 − (0)
0.0436 0.0007 0.001 −0.0002
[𝐷] = ( 0.0308 0.0055 0.0002 0.0002
)
0.0157 0.0051 −0.0012 −0.0002
−0.1440 −0.0063 0.0007 0.0001
Como ya se tienen la matriz de derivas de la estructura es posible determinar sus
valores máximos factibles, con la siguiente expresión (bajo el mismo concepto
aplicado en el cálculo de máximos factibles de desplazamiento):

𝐷𝑗𝑚á𝑥 = √∑(𝐷𝑗𝑖 )2
𝑖

𝐷4𝑚á𝑥 = √0.04362 + 0.00072 + 0.0012 + 0.00022 = 0.0436 𝑚

𝐷3𝑚á𝑥 = √0.03082 + 0.00552 + 0.00022 + 0.00022 = 0.0313 𝑚

𝐷2𝑚á𝑥 = √0.01572 + 0.00512 + 0.00122 + 0.00022 = 0.0166 𝑚

𝐷1𝑚á𝑥 = √0.14402 + 0.00632 + 0.00072 + 0.00012 = 0.1441 𝑚

El vector de derivas de piso máximas factibles es:
0.0436
{𝐷𝑚á𝑥 } = ( 0.0313 ) (𝑚)
0.1283
0.1441
Fuerzas Modales
Las fuerzas modales se obtienen con la siguiente expresión, [𝐹] = [𝐾][𝑈]
1.6892 −1.3082 0.3636 −0.0513
[𝐹] = 10000 ∗ ( −1.3082 2.1441 −1.3738 0.2959
)
0.3636 −1.3738 1.9649 −0.9110
−0.0513 0.2959 −0.9110 0.6601
−0.0539 0.0050 0.0007 −0.0001
−0.0975 0.0043 −0.0003 0.0001
∗( )
−0.1283 −0.0012 −0.0005 −0.0001
−0.1440 −0.0063 0.0007 0.0001
−27.4018 27.9095 13.3976 −3.1762
[𝐹] = (−49.5726 23.6530 −6.2054 5.1087
) (𝑡𝑜𝑛)
−65.2002 −6.7858 −10.4675 −4.5759
−42.6837 −20.4296 8.5139 1.8058
Se procede a calcular los valores máximos factibles de la siguiente forma:

𝐹𝑗𝑚á𝑥 = √∑(𝑓𝑗𝑖 )2
𝑖

𝐹4𝑚á𝑥 = √27.40182 + 27.90952 + 13.39762 + 3.17622 = 41.4654 𝑡𝑜𝑛

𝐹3𝑚á𝑥 = √49.57262 + 23.65302 + 6.20542 + 5.10872 = 55.5114 𝑡𝑜𝑛

𝐹2𝑚á𝑥 = √65.20022 + 6.78582 + 10.46752 + 4.57592 = 66.5404 𝑡𝑜𝑛

𝐹1𝑚á𝑥 = √42.68372 + 20.42962 + 8.51392 + 1.80582 = 48.1146 𝑡𝑜𝑛

El vector de fuerzas modales máximas factibles es:

41.4654
{𝐹𝑚á𝑥 } = ( 55.5114 ) (𝑡𝑜𝑛)
66.5404
48.1146
Cortantes Modales de Piso
Finalmente se calculan los valores de cortante modal de piso, el cual se define
como:
𝑚

𝑉𝑗𝑖 = ∑ 𝐹𝑘𝑖
𝑗

−27.4018 27.9095 13.3976 −3.1762

−76.9744 51.5625 7.1922 1.9325
𝑉=( )
−142.1746 44.7767 −3.2753 −2.6434
−184.8583 24.3471 5.2386 −0.9376