FACULTAD DE INGENIERÍA

CARRERA DE INGENIERÍA CIVIL

Curso: INGENIERÍA SISMORRESISTENTE

Mg. Ing. Omar Olivos

Lima - Perú

UNIDAD 3
RESPUESTA ELÁSTICA DE ESTRUCTURAS PARA SISTEMAS DE VGDL
Sesión 5: Análisis Sísmico Dinámico Modal Espectral – Seudo Tridimensional

UNIDAD 3 RESPUESTA ELÁSTICA DE

ESTRUCTURAS PARA
SEMANA 04 Y 05 SISTEMAS DE VGDL –
MÉTODO DINÁMICO
LOGRO ESPECÍFICO
DE APRENDIZAJE

Al finalizar la unidad, el estudiante estructura con criterio

sismorresistente sistemas lineales, haciendo uso de los conceptos
dinámicos, análisis modal, análisis espectral, en función de la
normatividad vigente.

1
 CONTENIDO DE LA UNIDAD 3

I. RESPUESTA Y ANÁLISIS SÍSMICO DE SISTEMAS DE VGDL

II. SISTEMAS GENERALIZADOS
III. ANÁLISIS SÍSMICO MODAL ESPECTRAL
IV. ANÁLISIS ESTÁTICO
V. ANÁLISIS SÍSMICO MODAL ESPECTRAL
SEUDOTRIDIMENSIONAL

LOGRO DE LA SESIÓN

Al finalizar la sesión los alumnos aprenderán la

respuesta dinámica de estructuras de varios
niveles mediante el análisis dinámico MODAL
pseudo tridimensional considerando tres
grados de libertad por nivel o piso.

2
TEMA
ANÁLISIS SÍSMICO MODAL SEUDO TRIDIMENSIONAL

CONTENIDO DEL TEMA

1. Introducción
2. Hipótesis de análisis
3. Procedimiento
4. Análisis sísmico seudo tridimensional
5. Análisis sísmico modal seudo tridimensional
6. Ejemplos

3
1. INTRODUCCIÓN

El análisis debe tener en cuenta todos los GDL necesarios para

representar completamente todos los modos posibles de deformación y
las fuerzas de inercia significativas.
Las ecuaciones de movimiento tienen la misma forma matricial que en
los casos básicos anteriores.
Las matrices de masas y de rigideces tienen usualmente más elementos
y pueden no ser necesariamente diagonales.
No parece ser necesario refinar demasiado un modelo elástico que
representa a una estructura que se espera incursione en
comportamiento no lineal para la intensidad del sismo de diseño.
Por ello, conviene usar tantos GDL como sean necesarios para
representar las deformaciones y fuerzas relevantes, salvo estructuras
especiales o de comportamiento complejo.

2. HIPÓTESIS DE ANÁLISIS

Una estructura espacial puede ser modelada como un ensamble de

pórticos planos, con propiedades de rigidez solamente en sus planos
respectivos, admitiendo que las rigideces ortogonales a sus planos son
bastante menores y pueden no ser consideradas.
La hipótesis fundamental es la relativa a las losas de piso, las cuales son
consideradas como cuerpos rígidos que conectan a los pórticos. Para
fines del análisis sísmico, los grados de libertad para las losas de piso son
tres: dos traslaciones horizontales y una rotación torsional en planta.
De este modo, tampoco se toman en cuenta las deformaciones axiales
en las columnas.

4
2. HIPÓTESIS DE ANÁLISIS

Tres GDL por piso o nivel: dos traslaciones y una rotacional en planta

3. PROCEDIMIENTO

3.1 DETERMINACIÓN DE LA MATRIZ DE MASAS

En las ecuaciones de movimiento, las fuerzas de inercia deben incluir los

productos de los momentos de inercia de las masas respecto a un eje
vertical con las aceleraciones rotacionales correspondientes.
La matriz de masas adopta la forma:

5
3. PROCEDIMIENTO

3.2 DETERMINACIÓN DE MATRIZ DE RIGIDEZ DE LA ESTRUCTURA

La matriz de rigidez de la estructura puede formularse mediante el
método de la rigidez directa, implementado como sigue:

Paso 1. DETERMINACIÓN DE LA RIGIDEZ LATERAL DE CADA PÓRTICO PLANO

Cada pórtico está en unas coordenadas locales, en base a los cuales
están referidos los grados de libertad considerados para el pórtico.

donde:
uij : desplazamiento del nivel j del pórtico i
i : indicador del pórtico
j : indicador del piso

3. PROCEDIMIENTO

Paso 1. DETERMINACIÓN DE LA RIGIDEZ LATERAL DE CADA PÓRTICO PLANO

Para el i-ésimo pórtico se determina mediante el siguiente
procedimiento:
a) Defina el grado de libertad para i-ésimo pórtico: los desplazamientos laterales
al nivel de los pisos ui = [ ui1 ui2 …… uiN ]T y los desplazamientos verticales y
rotacionales para cada nodo.
b) Obtenga la matriz de rigidez completa para el i-ésimo pórtico con referencia al
grado de libertad del pórtico.
c) Condense estáticamente todos los grados rotacionales y verticales para
obtener la matriz de rigidez lateral del i-ésimo pórtico con tamaño N x N, en
indicada mediante kxi o kyi si el pórtico está orientado en la dirección X o Y,
respectivamente.
uiN(t)

ui3(t)

GDL del
ui2(t)
Pórtico i
ui1(t)

6
3. PROCEDIMIENTO

3.2 DETERMINACIÓN DE MATRIZ DE RIGIDEZ

PASO 2. DETERMINACIÓN DE LA MATRIZ DE TRANSFORMACIÓN DEL DESPLAZAMIENTO
La matriz de transformación del desplazamiento que relaciona el grado
de libertad ui, definido en el paso 1(a) para el i-ésimo pórtico con el
grado de libertad global del edificio. Esta matriz es de N x 3N

uij = uoj cos αi + voj sen αi + θojri

ri = (xi - xo) sen αi - (yi - yo) cos αi

3. PROCEDIMIENTO

3.2 DETERMINACIÓN DE MATRIZ DE RIGIDEZ

PASO 2. DETERMINACIÓN DE LA MATRIZ DE TRANSFORMACIÓN DEL DESPLAZAMIENTO

Se determina la matriz de transformación que relaciona el

desplazamiento lateral
uij = gij . uoj
Sea el vector ui , que define los desplazamientos locales de todos los
pórticos:

ui = Gi uo

Gi es la matriz de transformación y uo la matriz de desplazamientos

(del centro de masas da cada piso) de toda la estructura.

7
3. PROCEDIMIENTO

3.2 DETERMINACIÓN DE MATRIZ DE RIGIDEZ

PASO 3. DETERMINACIÓN DE LA MATRIZ DE RIGIDEZ DE PÓRTICO TRASFORMADA A
GDL DEL EDIFICIO
La matriz Ki de 3N x 3N es la contribución del i-ésimo pórtico a la matriz
de rigidez del edificio

Ki = GiT . kLi . Gi

PASO 4. MATRIZ DE RIGIDEZ DEL EDIFICIO

Sumar las matrices de rigidez para todos los pórticos a fin de obtener
la matriz de rigidez del edifício:

K = Σ Ki

3. PROCEDIMIENTO

3.2 DETERMINACIÓN DE MATRIZ DE RIGIDEZ

PASO 4. MATRIZ DE RIGIDEZ DEL EDIFICIO

Sumar las matrices de rigidez para todos los pórticos a fin de obtener
la matriz de rigidez del edifício:

K = Σ ki =

8
3. PROCEDIMIENTO

3.3 PLANTEAMIENTO DE LA ECUACIÓN DE EQUILBIRIO

Para el piso j:

El vector de cargas para el sistema será:

F=
…

4. ANÁLISIS SÍSMICO SEUDO TRIDIMENSIONAL

4.1 ECUACIÓN DEL MOVIMIENTO

La ecuación de equilibrio dinámico tiene la siguiente forma:

+ + =− . . !

Cambio de variable u= Y: "# + $# + # = −". %. #!

En las ecuaciones de movimiento, mientras las aceleraciones de los

pisos ocurren en las dos direcciones horizontales y tienen un
componente rotacional, las aceleraciones del terreno ocurren
solamente en la dirección de análisis.

Por ello, en el segundo miembro de la ecuación, el vector R contiene

unos (1) en los lugares correspondientes a los GDL orientados en la
dirección de análisis y ceros (0) en los demás lugares.

De este modo, el movimiento del suelo aparece en las ecuaciones

correspondientes a la dirección del movimiento de la base.

9
4. ANÁLISIS SÍSMICO SEUDO TRIDIMENSIONAL

4.1 ECUACIÓN DEL MOVIMIENTO

Al trabajar por separado en cada dirección se tiene para:

Fuerzas solo en “ x ” : M# + $# + # = −". % . #!

Fuerzas sólo en “ y” "# + $# + # = −". % . #!

RX Ry

4. ANÁLISIS SÍSMICO SEUDO TRIDIMENSIONAL

4.2 ECUACIÓN DEL MOVIMIENTO

Ejemplo en el edificio de 2 pisos

' 1 0
' 1 0
(' 0 1
& = (' % = 0 % = 1
)' 0 0
)' 0 0

10
5. ANÁLISIS MODAL SEUDO TRIDIMENSIONAL

5.1 FACTOR DE PARTICIPACIÓN MODAL

En las ecuaciones de movimiento, mientras las aceleraciones de los
pisos ocurren en las dos direcciones horizontales y tienen un
componente rotacional, las aceleraciones del terreno ocurren
solamente en la dirección de análisis.
En un análisis dinámico modal, el factor de participación estática
Representará la dirección del sismo.

Para el análisis seudo tridimensional se tiene:

.,/ " % .,/ " %

Γ,, = Γ,, =
.,/ " ., .,/ " .,
Donde:
Rx, Ry son los vectores de localización en las direcciones X y Y,
respectivamente

5. ANÁLISIS MODAL SEUDO TRIDIMENSIONAL

5.2 MATRIZ DE MASAS

Matriz de masas

Jo es la matriz de momentos polares de inercia de la masa.

11
5. ANÁLISIS MODAL SEUDO TRIDIMENSIONAL

5.3 MATRIZ DE RIGIDECES

K=

6. EJEMPLO 1

ANÁLISIS ESTÁTICO SEUDO TRIDIMENSIONAL

Se desarrollará el análisis de una estructura aporticada de concreto armado
(E=2x106t/m2) de un piso, considerando que el diafragma de piso es
suficientemente rígido para realizar el análisis con un modelo Seudo
Tridimensional.

12
6. EJEMPLO 1

ANÁLISIS ESTÁTICO SEUDO TRIDIMENSIONAL

13
6. EJEMPLO 1

6. EJEMPLO 1

14
6. EJEMPLO 1

6. EJEMPLO 1

15
6. EJEMPLO 1

6. EJEMPLO 1

Gi
KL i
Póritco i (t/cm) cos αi sen αi ri (cm)
1 10.325 1 0 450
2 10.325 1 0 -150
3 6.848 1 0 -750
A 10.325 0 1 -450
B 10.325 0 1 150
C 6.848 0 1 750

T
Ki = Gi . KL i . Gi

16
6. EJEMPLO 1

Rigidez global del pórtico 1

1
K1 = 0 10.325 1 0 450
450

10.325
K1 = 0 1 0 450
4646.25 3x1 1X3

10.325 0 4646.25
K1 = 0 0 0
4646.25 0 2090812.5
3x3

6. EJEMPLO 1

Rigidez global del pórtico 2

1
K2 = 0 10.325 1 0 -150
-150

10.325
K2 = 0 1 0 -150
-1548.75 3x1 1X3

10.325 0 -1548.75
K2 = 0 0 0
-1548.75 0 232312.5
3x3

17
6. EJEMPLO 1

Rigidez global del pórtico 3

1
K3 = 0 6.848 1 0 -750
-750

6.848
K3 = 0 1 0 -750
-5136 3x1 1X3

6.848 0 -5136
K3 = 0 0 0
-5136 0 3852000
3x3

6. EJEMPLO 1

Rigidez global del pórtico A

0
KA = 1 10.325 0 1 -450
-450

0
KA = 10.325 0 1 -450
-4646.25 3x1 1X3

0 0 0
KA = 0 10.325 -4646.25
0 -4646.25 2090812.5
3x3

18
6. EJEMPLO 1

Rigidez global del pórtico B

0
KB = 1 10.325 0 1 150
150

0
KB = 10.325 0 1 150
1548.75 3x1 1X3

0 0 0
Kb = 0 10.325 1548.75
0 1548.75 232312.5
3x3

6. EJEMPLO 1

Rigidez global del pórtico C

0
KC = 1 6.848 0 1 750
750

0
KC = 6.848 0 1 750
5136 3x1 1X3

0 0 0
KC = 0 6.848 5136
0 5136 3852000
3x3

19
6. EJEMPLO 1

K edificio = K1 + K2 + K3 + KA + KB + KC

27.5 0.0 -2038.5

K edificio = 0.0 27.5 2038.5
-2038.5 2038.5 12350250.0

6. EJEMPLO 1

Vector de fuerzas y vector de desplazamientos

10 uo
F= 0 Uo = vo
0 θo

c) Determinación de desplazamientos de pisos

-1
[ F ] = [ K ] [ Uo ] → [ Uo ] = [ K ] [ F ]

uo 0.03682243 -0.00045615 6.1531E-06 10 0.36822 cm

Uo = vo = -0.00045615 0.03682243 -6.1531E-06 0 = -0.00456 cm
θo 6.1531E-06 -6.1531E-06 8.30013E-08 0 0.000062 rad

20
6. EJEMPLO 1

6. EJEMPLO 1

21
6. EJEMPLO 1

6. EJEMPLO 1

22
TEMA
ANALISIS SISMICO DINAMICO SEUDO-
TRIDIMENSIONAL

45

6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

Para el edificio de un piso mostrado en la figura, se pide realizar el análisis sísmico

matricial seudo tridimensional usando el método de las rigideces. Aplicar el
método de Análisis sísmico dinámico (modal espectral) conforme a la Norma E.030

23
6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

Vigas: 0.30 x 0.60 m

Columnas: 0.30 x 0.30 m

DATOS GENERALES: DATOS DE DISEÑO:

- Ubicación : Lima, Pueblo Libre f’c = 210 kg/cm2
- Uso : Vivienda fy = 4200 kg/cm2
- Sistema estructural : Aporticado Ec = 15100√ 210 = 217,371 tn/cm2
- Suelo de cimentación : Tipo S1(suelo rígido)

6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

1) Cálculo de matriz de rigideces laterales locales de los pórticos

Por Wilbur: (Base Empotrada)

24
6. EJEMPLO 2

ANÁLISIS SÍSMICO SEUDO TRIDIMENSIONAL

1) Cálculo de matriz de rigideces laterales locales de los pórticos

6. EJEMPLO 2

ANÁLISIS SÍSMICO SEUDO TRIDIMENSIONAL

1) Cálculo de matriz de rigideces laterales locales de los pórticos

25
6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

2) Centro de masas
VALOR APRÓXIMADO: Centro geométrico del área en planta.

CM (Xo,Yo)

6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

2) Centro de Rigidez (CR)

26
6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

4) Excentricidad (e)

6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

4) Excentricidad (e)

27
6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

5) Momentos Torsores: Mt

6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

6) MATRIZ DE MASAS
Para “n” pisos:
Piso n

Piso 3

Para el piso “i”:

28
6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

6) MATRIZ DE MASAS
Para “n” pisos: Para el piso “i”:

6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

6) MATRIZ DE MASAS
Para el cálculo del momento polar de masas, es frecuente aproximarlo
considerando que las masas de cada piso están aproximadamente
uniformemente distribuidas, y por consiguiente los momentos polares de
inercia pueden determinarse a partir de las dimensiones del área de losa en
planta, suponiendo que el radio de giro de las masas es el mismo que el de las
áreas.

29
6. EJEMPLO 2

ANÁLISIS SÍSMICO DINÁMICO SEUDO TRIDIMENSIONAL

6) MATRIZ DE MASAS

6. EJEMPLO 2
6) Matriz de masas: M

30
6. EJEMPLO 2
7) Cálculo de la matriz de rigidez de la estructura: K

La rigidez en coordenadas globales del pórtico “i”

, = 0,/ · 2, · 0,
Resolviendo para el pórtico “i”:

La matriz de rigidez en coordenadas globales del edificio será la suma de los Ki

4 4

=3 , = 3 0,/ · 2, · 0,
,5 ,5

6. EJEMPLO 2
7) Cálculo de la matriz de rigidez de la estructura: K

31
6. EJEMPLO 2
7) Cálculo de la matriz de rigidez de la estructura: K

6. EJEMPLO 2
7) Cálculo de la matriz de rigidez de la estructura: K

32
6. EJEMPLO 2
7) Cálculo de la matriz de rigidez de la estructura: K

6. EJEMPLO 2
8) Cálculo de los periodos y formas de modo de vibración:

33
6. EJEMPLO 2
8) Cálculo de los periodos y formas de modo de vibración:
Para nuestro caso, como nuestro edificio es de un solo nivel los factores de
participación modal serán:
.,/ " %
Γ,, =
.,/ " .,

Las formas de modo φi han sido normalizados de manera que el producto:

6. EJEMPLO 2
8) Cálculo de los periodos y formas de modo de vibración:

34
6. EJEMPLO 2
9) Cálculo de la aceleración espectral: Sai

6. EJEMPLO 2
9) Cálculo de la aceleración espectral: Sai

35
6. EJEMPLO 2
10) Cálculo del vector de desplazamientos u

6. EJEMPLO 2
10) Cálculo del vector de desplazamientos u

36
6. EJEMPLO 2
10) Cálculo del vector de desplazamientos u

6. EJEMPLO 2
10) Cálculo del vector de desplazamientos u

37
6. EJEMPLO 2
10) Cálculo del vector de desplazamientos u

6. EJEMPLO 2
10) Cálculo del vector de desplazamientos u

38
6. EJEMPLO 2
10) Cálculo del vector de desplazamientos u

¿ PREGUNTAS ?

39
 CONCLUSIONES

¿Qué hemos aprendido el día de hoy?

La respuesta en una estructura debido a fuerzas sísmicas depende de sus

propiedades dinámicas (periodos, formas modales y amortiguamiento). En los
análisis es común obtener esta respuesta por medio del método de
"Superposición Modal", el cual dice que la respuesta total en la estructura es
la suma de las respuestas individuales de cada modo.

A continuación se muestra un video que muestra la superposición modal de

forma conceptual, donde se observa que algunos modos se excitan más que
otros, esto es debido a que el movimiento sísmico no contiene frecuencias
que puedan excitar tales modos.

40