INGENIERÍA SISMORRESISTENTE

MÉTODOS PARA ANÁLISIS MODAL

Ms. Ing. Gonzalo H. Díaz García
 ASIGNACION DEL TRABAJO DE CAMPO

[Link]án conformar grupos de 6 estudiantes uno de los cuales será el coordinador.

[Link] trabajo de campo consistirá de un análisis sísmico de una edificación usando el ETABS y
sus soluciones correspondientes. Los entregables se detallan a continuación:

a) Segundo entregable (T2): Análisis Modal espectral modelo 3D. Presentación de modelos
en software
RESPUESTA SÍSMICA DE SISTEMAS DE VARIOS GRADOS DE
LIBERTAD Y CRITERIOS ESTRUCTURALES SISMORRESISTENTES

Al finalizar la unidad, el estudiante estructura con criterio sismorresistente sistemas lineales

y de varios grados de libertad, haciendo uso de los conceptos dinámicos, análisis modal,
análisis espectral, en función de la normatividad vigente
 LOGRO:

• Al termino de la sesión el estudiante, de forma colaborativa, determinar

las fuerzas laterales y desplazamientos de cada nivel usando el
espectro de respuesta; según el caso de estudio asignado por su
docente, mostrando en el proceso criterios de orden y precisión.
 CONTENIDO

1.Métodos para análisis modal de sistemas de “n” gdl.

2.Métodos manuales.
3.Métodos matriciales.
[Link] de los modos de vibración.
[Link]ón de los modos de vibración.
INTRODUCCIÓN

Resulta complejo elegir entre el análisis dinámico plano o tridimensional, éste

último representado por sus dos componentes horizontales (cargas
reversibles), lo cual no es posible en el plano. El análisis dinámico
tridimensional, requerirá la evaluación de la estructura con varios grados de
libertad por medio de métodos sofisticados como el de los elementos finitos,
que ayudaría a resolver las ecuaciones diferenciales de movimiento existentes
por cada grado de libertad, es una herramienta poderosa, sin embargo su
modelación e interpretación de resultados no es sencilla y sólo se justificaría
su uso en obras de magnitud.
ECUACIÓN DE MOVIMIENTO

• Una estructura de varios niveles mostrada en la Figura 1, se puede idealizar como un

pórtico de varios niveles con diafragma de cuerpo rígido asumiendo que la masa está
concentrada en cada nivel, las columnas se suponen axialmente inextensibles pero
lateralmente flexibles. La respuesta dinámica del sistema está representada por el
desplazamiento lateral de las masas con el número de grados de libertad dinámica o
n modos de vibración que son iguales al número de masas. La vibración resultante
del sistema esta dada por la superposición de las vibraciones de cada masa. Cada
modo individual de vibración tiene su propio periodo y puede ser representado por un
sistema simple del mismo periodo.

Pórtico Sistema Modo Modo Modo Modo Modo

Equivalente 1 2 3 4 5
ECUACIÓN DE MOVIMIENTO

Las fuerzas de inercia en la ecuación 10.1 son simplemente:

f Ia  m a  ua
f Ib  mb  ub (10
f Ic  m c  uc

En forma matricial:

 f Ia  m a 0 0  ua 
    
 f Ib    0 mb 0   ub  (10
f   0 mc   
 Ic   0 uc 
O más generalmente:

FI   M  U
ECUACIÓN DE MOVIMIENTO

En forma matricial:

 f Sa  k aa k ab k ac  u a 
    
 f Sb   k ba k bb k bc   u b  (10.6)
 f  k k cc   
 Sc   ca k cb u c 
O más generalmente:
Fs   K  U  (10.7)

Donde {Fs} es el vector de fuerzas elásticas, [K] es la matriz de rigidez y {U} es el vector de desplazamientos.

Por analogía, las fuerzas de amortiguamiento en la ecuación 10.1 pueden expresarse como:

FD   C U 

ECUACIÓN DE MOVIMIENTO


Donde {FD} es el vector de fuerzas de amortiguamiento, [C] es la matriz de amortiguamiento y U es el vector
de velocidades. En general no es práctico determinar c y el amortiguamiento es expresado en términos del
coeficiente de amortiguamiento ().

Aplicando las ecuaciones 10.4, 10.7 y 10.8 las ecuaciones de equilibrio dinámico (10.1) pueden escribirse
generalmente como:
{FI}+{FD}+{FS}={p(t)} (10.9)
Lo cual es equivalente a:
Lo  
M cual  
U  es
Cequivalente
U  K U   a:
p(t ) (10.10)
 
M  U  C U  K U   p(t )  
RESPUESTA DINÁMICA: ANÁLISIS MODAL

• Para determinar la respuesta dinámica de una estructura de varios grados de libertad se puede utilizar el
procedimiento de análisis modal. Se obtiene la respuesta máxima por separado para cada modo, modelando cada
uno de ellos como un sistema de simple grado de libertad. Debido a que los valores máximos no pueden ocurrir
simultáneamente, estos valores son combinados estadísticamente para obtener la respuesta total.

• El método de la raíz cuadrada de la suma de los cuadrados, SRSS, es aplicable para estructuras bidimensionales
cuando la relación entre los periodos de cualquier modo alto con cualquier modo bajo es 0.75 o menor, y la relación
de amortiguamiento no excede el 5%. El análisis modal puede ser enfocado mediante métodos matriciales,
numéricos o métodos iterativos.

Alan Williams, pp 21 [ref. 4]

MÉTODO MATRICIAL

• Como la respuesta dinámica de una estructura depende de la frecuencia o periodo

de vibración y de la forma desplazada (forma modal), el primer paso en un análisis
de un sistema de varios grados de libertad es encontrar las frecuencias y las formas
modales de vibración libre. En este caso no existen fuerzas externas y el
amortiguamiento es considerado cero.

• Cada grado de libertad dinámico provee una ecuación de equilibrio dinámico, la

vibración resultante del sistema consiste de n de éstas ecuaciones, y puede ser
expresado en forma matricial para vibración libre no amortiguada como:

M  U K  U   0
• La vibración libre descrita gráficamente por las Figuras mostradas de un sistema no
amortiguado en uno de sus modos de vibración natural puede describirse
matemáticamente por:
MÉTODO MATRICIAL

Vibración libre de un sistema no amortiguado en su primer modo natural de vibración: (a) Pórtico de dos
niveles; (b) Forma de la deformada en los instantes de tiempo a, b, c, d y e; (c) Coordenada modal q 1(t); (d)
Histograma de desplazamiento
MÉTODO MATRICIAL

Vibración libre de un sistema no amortiguado en su segundo modo natural de vibración: (a) Pórtico de dos
niveles; (b) Forma de la deformada en los instantes de tiempo a, b, c, d y e; (c) Coordenada modal q 2(t); (d)
Histograma de desplazamiento
MATRIZ MODAL Y ESPECTRAL

• Los N eigenvalores y los N modos pueden ser acoplados en forma matricial. El modo natural o eigenvector n
correspondiente a la frecuencia natural n tiene elementos jn, donde j indica el DOF. De este modo los N
eigenvectores pueden presentarse o disponerse en una matriz cuadrada, de la cual cada columna es un modo:

 11 12  1 N 
  22   2 N 
   jn   
 
21
   
 
 N1  N 2   NN 
• Donde [] es llamada matriz modal. Los N eigenvalores n2 pueden ser acoplados en una ,matriz diagonal 2, la
cual es conocida como matriz espectral.

 12 
 
  2


 22



 
  N 
2
MATRIZ MODAL Y ESPECTRAL

• Cada eigenvalor y eigenvector satisfacen la ecuación 3 la cual puede ser reescrita como:

K   n   M   n  n2
• Utilizando la matriz modal y espectral es posible representar esta ecuación en una ecuación matricial simple:

K    M    2 
• Esta ecuación presenta en forma compacta las ecuaciones relacionando todos los eigenvalores y eigenvectores.
Ortogonalidad de los modos

• Los modos naturales correspondientes a diferentes frecuencias naturales se muestran a continuación para
satisfacer la siguiente condición de ortogonalidad. Cuando nr (entiéndase que r también es una frecuencia
natural).

 nT  k   r  0  nT  m   r  0
• La demostración de esta propiedad es la siguiente: la enésima frecuencia natural y el modo que satisfacen la
ecuación 3 multiplicados por rT, la transpuesta de r, da:

 rT  k   n   n2   rT  m   n
• Análogamente se realiza lo mismo con la erésima frecuencia natural y el modo que satisface la ecuación 3; de esta
manera k·r = r2·m·r multiplicando por nT da:

 nT  k   r   r2   nT  m   r
• Donde se ha utilizado la propiedad de simetría de la matriz de masa y rigidez. se tiene:

( n2   r2 )   nT  m   r  0
ORTOGONALIDAD DE LOS MODOS

• Debido a que m y k son definidos positivos, los elementos de la diagonal de K y M son positivos, y están
relacionados por:

kn  n
2
 mn
MÉTODO ITERATIVO

• Para edificios de pocos niveles, que no excedan a cinco plantas, el análisis modal
puede limitarse al modo fundamental. El sistema estructural puede ser modelado
como un pórtico con losas de entre piso rígidas. Los desplazamientos laterales de los
nudos son entonces el resultado de la flexión de las columnas sin incluir rotación en
los nudos. La rigidez de un nivel en particular esta dada por:
• La masa en cada nivel se asume concentrada en las losas de entre piso como se
muestra en la Figura 10.4. utilizando estos supuestos se han desarrollado técnicas
iterativas basadas en métodos propuestos por Rayleigh, Stodola y Holzer. A
continuación se presenta una adaptación del método de Holzer. El modelo dinámico
describe que: cuando un nudo alcanza su desplazamiento lateral máximo ui, la
velocidad es cero y la fuerza de inercia en el nudo está dada por:

I FI  mi  ui
k i  12  E 
h3 FI  mi   2  u i

Alan Williams, pp 35-40 [ref. 4]

MÉTODO ITERATIVO

Análisis modal de una estructura resistente a fuerza lateral

m 1 2 x1
m 3 2 x3 m 2 2 x2

m3 x3
m 3 2 x3
k3 k3  3
m2 x2
m 2 2 x2
k2 k2 2
m1 x1
m 1 2 x1
k1 k1 1
 1  2 3

Sistema Modelo Desplazamiento Fuerzas Esfuerzos

estructural dinámico del nudo laterales en de corte
equivalente equivalente y deriva los nudos de piso
MÉTODO ITERATIVO

• La fuerza cortante en cualquier nivel es igual al producto de la rigidez del nivel por el
desplazamiento del mismo. El incremento en la fuerza de corte en el nudo es
producido por la fuerza de inercia en ese nivel. El incremento de la fuerza cortante
esta dado por:
FS  k i   i  k i 1   i 1
• Donde ki·i = fuerza cortante total en el nivel i.
• Igualando la fuerza de inercia y el incremento de la fuerza cortante se tiene:

FI  FS
mi   2  u i  k i   i  k i 1   i 1
• La solución de esta ecuación se puede obtener asumiendo una forma modal inicial
con un desplazamiento unitario en el nivel superior; a partir del cual se calcula la
fuerza de inercia o el incremento de fuerza cortante en términos de la frecuencia
natural, en cada nivel. Sumando el incremento de fuerza cortante a partir del nivel
superior hacia abajo se tiene la fuerza cortante total en cada piso. Dividiendo este
valor por la rigidez apropiada de cada nivel se obtiene el desplazamiento (deriva) de
cada piso. Dividiendo estos desplazamientos por el desplazamiento en la parte
superior de la estructura se obtiene la forma modal corregida. Esta forma modal
corregida puede ser usada como una nueva forma modal inicial en el proceso de
iteración hasta que coincidan la forma modal corregida con la inicial.
APLICACIÓN 01:

• Un pórtico de dos niveles tiene las propiedades mostradas en la Figura A y una

relación de amortiguamiento de 5%, está localizado sobre un sitio rocoso cerca del
origen del sismo de Loma Prieta. Determinar las fuerzas laterales y desplazamientos
de cada nivel usando el espectro de respuesta.

1.5
w2 = 70 [t]
w2 = 70 [t]
4 m. k2 = 3 [t/cm] k2 = 3 [t/cm]
Roca, derca de la fuente w1 = 70 [t]
1.0
w1 = 70 [t]
Aceleración espectral, g

4 m. k1 = 7 [t/cm] k1 = 7 [t/cm]


Roca blanda, sur de San

Francisco

0.5

Roca suave, Oackland

Roca, San Francisco

0.0 1.0 2.0 3.0 4.0

Periodo natural de vibración T n, s
APLICACIÓN 01:

• El desplazamiento de cada nudo es:

 k11 k12  7  3  3  10  3
K        3 3 
 21
k k 22    3 3   
• la matriz diagonal de masa es:
w  m g  m 
w 70
 
 0.071 t  s 2 / cm  M   
0.071 0 
0.071
g 981  0
• la ecuación de eigenvalores es:

K   2
n 
 M    n   0   10  3

 3 3    2 0.071

n 
0   1n 
      0
   0 0.071   2n 

• el determinante a resolver es

10  0.071  n2 3 30  0.71  n2  0.213   n2  0.005   n4  9  0

0
3 3  0.071  n2
0.005   n4  0.923   n2  21  0
APLICACIÓN 01:

• resolviendo el polinomio se obtiene las frecuencias naturales correspondientes a los

modos de vibración:
1  5.15 rad / s
 2  12.57 rad / s
• el periodo natural correspondiente es: T1  1.22 s
T2  0.5 s 
• A partir del espectro de respuesta, la aceleración espectral es:
 
A1  0.23  g  225.63 cm / s 2
A2  0.83  g  814.23 cm / s 
2

• Sustituyendo estos valores en la ecuación de eigenvalores y estableciendo el primer

componente de cada modo igual a la unidad se obtiene cada uno de los
eigenvectores: Para w1=5.15

10  0.71(5.15) 2 3  11  0  8.1129  11  3   21  0

      
 3 3  0.071(5.15) 2   21  0  3  11  1.1169   21  0

 1 
11  1   21  2.706  
2.706
APLICACIÓN 01:

• Para w2=12.57 se resuelve análogamente y se obtiene el eigenvector:

 1 
 
 0.406
• De este modo se obtiene la matriz de eigenvectores:
 1 1 
2.706  0.406
 

• Los componentes de la matriz modal normalizada están dados por:

u jn
 jn 
m jj  u 2jn  1
2

• Para el modo 1:

m jj  u 2j1   0.071  1  2.706 
1
2 2
1
2
 0.768

• Para el modo 2:


m jj  u 2j 2   0.071  1  0.406 
1
2 2
1
2
 0.288
APLICACIÓN 01:

 1 1 
•  0.288   1.302 3.477 
Entonces la matriz modal normalizada es:
   02..768
706  0.406  3.520  1.413
 
 0.768 0.288 

• El vector de coeficiente de participación esta definido por:


P  TT  M  
1
  M  
• como para la matriz modal normalizada, entonces [P] se reduce a: P  T  M  1
P  
1.302 3.520  0.071 0  1 0.342
    1  0.147
3.477  1.413  0 0.071   
• Asumiendo que la estructura se comporta elásticamente, la matriz de desplazamiento
esta dada por:
U    P D   P A  2 
5.15 2 0 
U   
1.302 3.477  0.342 0  225.63 0 

    
3.520  1.413  0 0.147  0 814.23  0 12.57 2 

U   
3.78 2.616 

10.219  1.063
APLICACIÓN 01:

• Los desplazamientos máximos resultantes de cada nudo se obtienen a través de la raíz cuadrada de la suma de
los cuadrados, SRSS, de la fila respectiva a cada nudo y esta dado por:

U c   
4.60 

10.27
• La matriz de fuerzas laterales en cada nudo esta dado por:
Fs   K  U 
Fs   
7.143 29.349 
• 19.317 la
Utilizado la raíz cuadrada de la suma de los cuadrados .037 lateral en cada nudo es:
11fuerza

Fsc   
30.20
• El vector de cortante basal es: 
22.25

• V   F T
 1 basal
Utilizado la raíz cuadrada de la suma de los cuadrados la cortante Vbes:  26 .46
T
18 .312 

Vb  
V112  V122  1
2
 32.178 t 
ACTIVIDAD 2:

• Realizar el análisis del caso en estudio; determinando las fuerzas laterales y

desplazamientos de cada nivel usando el espectro de respuesta:
GRACIAS…