ANÁLISIS MODAL ESPECTRAL

U NI VE R SI DA D PE D A GÓGIC A Y TE C N OL ÓGI C A D E C O LO MB IA
E S PE CI A LI Z AC I ÓN EN ES T RU C TURAS
I NG . I VO NN E A LE J AND RA M ON RO Y M OL AN O
Adaptado del capítulo 15 del libro Dinámica Estructural Aplicada al Diseño Sísmico
 FORMULACIÓN DEL ANÁLISIS
MODAL ESPECTRAL
De acuerdo con lo presentado en las secciones anteriores, sabemos que las ecuaciones de movimiento para
un sistema sometido a una excitación en la base son de la forma:

𝑀 𝑈̈ + 𝐾 𝑈̇ = − 𝑀 𝛾 𝑋̈ (1)

En donde,
𝑀 : Matriz de masa de la estructura
𝐾 : Matriz de rigidez de la estructura
𝛾 : Matriz que se ve afectada por el número de componentes de acelerograma (pueden existir hasta tres
componentes del acelerograma). Puede tener dimensiones nx1, nx2 ó nx3 respectivamente, siendo n el
número de grados de libertad
𝑋̈ : Vector de aceleraciones en la base
Posteriormente por medio de la solución del problema de valores y vectores propios (eigen values), podemos
obtener los modos de vibración y las frecuencias de la estructura, Φ 𝑦 𝜔 , con base en sus propiedades de
vibración libre representadas por:

𝑀 𝑈̈ + 𝐾 𝑈̇ = 0

La solución del sistema de ecuaciones diferenciales simultáneas se obtiene desacoplando el sistema por medio
de las siguiente transformación de coordenadas:
𝑈 = Φ 𝜂 (2)
Y derivando dos veces contra el tiempo sabemos que:
𝑈̈ = Φ 𝜂̈ (3)

Reemplazando (2) y (3) en (1):

Φ 𝑀 Φ 𝜂̈ + Φ 𝐾 Φ 𝜂 =−Φ 𝑀 𝛾 𝑋̈

𝐼 𝜔 𝛼

Tanto 𝐼 como 𝜔 , son matrices diagonales, y por esto el sistema se desacopla, lo cual implica que tienen n
ecuaciones independientes de un grado de libertad del tipo:
𝜂 ̈ + 𝜔 𝜂 = −𝛼 𝑋 ̈
Y si se aplica amortiguamiento modal:
𝜂 ̈ + 2𝜉𝜔 𝜂̇ + 𝜔 𝜂 = −𝛼 𝑋 ̈

De acuerdo con la definición de espectro de respuesta, sabemos que el máximo valor que puede tener el
desplazamiento relativo u, entre la base y la masa de un sistema de un grado de libertad sometido a un
acelerograma en su base 𝑋̈ (𝑡), es precisamente el valor que se lee del espectro de desplazamiento Sd(T,ξ),
calculado para el mismo acelerograma, utilizando los valores del periodo T, y el amortiguamiento ξ, del mismo
sistema de un grado de libertada.

Por lo tanto el máximo valor que puede tomar 𝜂 corresponde al valor leído del espectro de desplazamientos
de la excitación amplificado por el coeficiente de participación 𝛼 .

Entonces:
(𝜂 ) = 𝛼 ∗ 𝑆 (𝑇 , 𝜉 )

Recordemos que:
2𝜋
𝑇 =
𝜔
Adicionalmente
𝜉 : Corresponde al valor del amortiguamiento modal
Cuando se dispone del espectro de aceleraciones puede utilizarse la transformación determinada a partir de la
relación entre Sa, Sv y Sd:
1 𝑇
(𝜂 ) = 𝛼 ∗ ∗ 𝑆 (𝑇 , 𝜉 ) = 𝛼 ∗ ∗ 𝑆 (𝑇 , 𝜉 )
𝜔 4𝜋

Dado que, por medio de cualquiera de los dos procedimientos, se dispone de unos valores máximos de los
desplazamientos de los grados de libertad de la estructura 𝑈 . Desafortunadamente, este procedimiento es
errado debido a que los valores máximos de los desplazamientos, o aceleraciones, que se coleccionan en el
espectro de respuesta no ocurren en el mismo instante en el tiempo.

Adicionalmente, debe anotarse que el signo, positivo o negativo, que indica la dirección de la respuesta
también se pierde, debido a que al espectro se llevan únicamente los valores absolutos de las respuestas.

En principio, las respuestas modales individuales son correctas y corresponden a los máximos valores que
tendrían cada una de ellas, simplemente hay que tener en cuenta que pueden ser tanto positivas como
negativas. La dificultad radica en determinar una manera apropiada de combinarlas para obtener una
respuesta modal apropiada. Esto se logra por medio de lo que se conoce como Métodos de Combinación
Modal Espectral.
Los desplazamientos dinámicos modales máximos que se presentan en la estructura, correspondientes a cada
modo individual, modo i, pueden obtenerse por medio de:
𝑼𝒎𝒐𝒅 (𝒊) = 𝝓(𝒊) ∗ 𝜼𝒊 𝒎𝒂𝒙 = 𝝓 𝒊 ∗ 𝜶𝒊 ∗ 𝑺𝒅 𝑻𝒊 , 𝝃𝒊

En la ecuación anterior debe tenerse en cuenta que el resultado multiplicado por -1 también es factible, dado
que se trata de un movimiento alternante derivado de un fenómeno ondulatorio. Esta posibilidad de cambio
de signo se manifiesta en todas las diferentes formas de la respuesta modal.

Para cada modo individual (i), las fuerzas dinámicas inerciales modales máximas que se presentan en la
estructura pueden obtenerse multiplicando los desplazamientos modales máximos por la matriz de rigidez de
la estructura:
𝑭𝒎𝒐𝒅 (𝒊) = 𝑲 𝑼𝒎𝒐𝒅 (𝒊) = 𝑲 𝝓(𝒊) (𝜼𝒊 )𝒎𝒂𝒙 = 𝑲 𝝓(𝒊) 𝜶𝒊 ∗ 𝑺𝒅 (𝑻𝒊 , 𝝃𝒊 )

Cada una de estas fuerzas modales máximas pueden utilizarse como un conjunto de fuerzas estáticas y con
ellas, independientemente, por medio de un análisis estático convencional, llegar a encontrar las fuerzas
internas causadas por el modo (i) en cada uno de los elementos de la estructura. Análogamente, estas fuerzas
internas modales máximas pueden obtenerse, también, utilizando los desplazamientos modales máximos
obtenidos. Las dos alternativas conducen a resultados idénticos.
En este punto se tendría la respuesta máxima, individual por modo, de los diferente parámetros relevantes
causados por unas fuerzas inerciales aplicadas a la estructura como si fueran fuerzas estáticas externas. Estos
parámetros comprenden las fuerzas internas en los elementos de la estructura, las derivas de piso, el corte
basal y el momento de vuelco, entre otros. Habrá tantos conjuntos independientes de parámetros como
modos tenga la estructura. Tan solo bastaría combinarlos.
EJEMPLO

Mg

kN/m
Por la tanto las ecuaciones del movimiento son:

𝑀 𝑈̈ + 𝐾 𝑈̇ = − 𝑀 𝛾 𝑋̈

𝑈 ̈ 𝑈
24 0 0 0 0 0 23161.854 0 0 −28169.3 0 0
0 24 0 0 0 0 𝑈 ̈ 0 24371.97 0 0 −28952 0 𝑈
0 0 122 0 0 0 𝑈 ̈ 0 0 364109.3 0 0 −436626 𝑈
+
0 0 0 24 0 0 𝑈 ̈ −28169.348 0 0 62166.82 0 0 𝑈
0 0 0 0 24 0 𝑈 ̈ 0 −28952 0 0 63127.05 0 𝑈
0 0 0 0 0 122 0 0 −436626 0 0 956686 𝑈
𝑈 ̈
24 0 0 0 0 0 1
0 24 0 0 0 0 1
=− 0 0 122 0 0 0 1 𝑋̈
0 0 0 24 0 0 1
0 0 0 0 24 0 1
0 0 0 0 0 122 1
Resolviendo el problema de valores propios tenemos que:

Modo ω2 (rad/s)2 ω (rad/s) f (Hertz) T (s)

1 9738.21 98.682 15.706 0.064
2 3274.49 57.223 9.107 0.110
3 3205.25 32.985 9.011 0.111
4 1087.98 32.985 5.250 0.190
5 371.30 19.269 3.067 0.326
6 350.11 18.711 2.978 0.336
Normalizando los modos con respecto a la masa tenemos que:

Verificando la normalización de los modos tenemos:

𝐼 𝜔
Graficando los modos normalizados:
Posteriormente se calculan los coeficientes de participación:

Con esto ya podemos determinar las ecuaciones desacopladas::

𝜂 + 2𝜉 𝜔 𝜂 ̇ + 𝜔̈ 𝜂 = −4.588𝑋̈

𝜂 + 2𝜉 𝜔 𝜂 ̇ + 𝜔̈ 𝜂 = 2.014𝑋̈

𝜂 + 2𝜉 𝜔 𝜂 ̇ + 𝜔̈ 𝜂 = −2.066𝑋̈

𝜂 + 2𝜉 𝜔 𝜂 ̇ + 𝜔 ̈ 𝜂 = −14.932𝑋̈
̈
𝜂 + 2𝜉 𝜔 𝑣𝜂 ̇ + 𝜔 𝜂 = 6.629𝑋̈

𝜂 + 2𝜉 𝜔 𝜂 ̇ + 𝜔̈ 𝜂 = −6.613𝑋̈
La respuesta para cada una de estas ecuaciones desacopladas, se obtiene utilizando el espectro de
desplazamientos. Para el caso particular utilizaremos el espectro del NSR-10 para la ciudad de Tunja con un
perfil de suelo tipo D y un grupo de uso I.

Espectro Elástico de Desplazamientos de Diseño

0.700

0.600

0.500

0.400
Sd (m)

0.300

0.200

0.100

0.000
0 1 2 3 4 5 6 7
T (s)
0.62*Aa*Fa*I*T2*(0.4+0.6*T/T0) 0.62*Aa*Fa*I*T2 0.3*Av*Fv*I*T 0.3*Av*Fv*I*TL
Por lo tanto tenemos que:

MODO T (s) αi Sd(Ti,ξi) (m) (ηi)max=αixSd(Ti,ξi)

1 0.064 -4.588 0.000498 -0.00228
2 0.110 2.014 0.00181 0.00364
3 0.111 -2.066 0.00185 -0.00382
4 0.190 14.932 0.00627 -0.09362
5 0.326 6.629 0.01845 0.12231
6 0.336 -6.613 0.01860 -0.1230
Determinación de los desplazamientos por modo:

𝑈( ) 𝑈( ) 𝑈( ) 𝑈( ) 𝑈( ) 𝑈( )
Graficando los desplazamientos por modo:
Derivas de piso máximas:

Utilizando los desplazamientos anteriores es posible calcular para cada modo la deriva máxima que puede
tener cada piso de la estructura como el desplazamiento relativo entre piso y piso. Es costumbre expresar esta
deriva como porcentaje de la altura de cada piso.

Fuerzas inerciales máximas

Para determinar las fuerzas inerciales máximas por modo que impone el sismo sobre la estructura, se multiplica
la matriz de rigidez de toda la estructura por los desplazamientos máximos correspondientes a cada modo:
Cortante máximo modal de piso

El cortante máximo modal de piso se define como:

CORTANTE BASAL