DINAMICA DE ESTRUCTURAS Profesor: JOHN P.

ARAGON BROUSSET

5.5. Espectros de respuesta

Si observamos la integral de Duhamel (definida para sistemas de 1 GDL), ésta depende,

exclusivamente, del acelerograma elegido (𝑢̈ 𝑔) y de las características vibratorias del
sistema expresadas a través de su frecuencia de vibración “w” o de su periodo de
2𝜋
vibración “T” (𝑤 = √𝐾⁄𝑚 y 𝑇 = ). Luego, los periodos fundamentales de vibración de
𝑤
los edificios comúnmente construidos están, aproximadamente, entre 0 y 4 seg. Si
calculamos la integral de la ecuación (5.26) para este rango de periodos y graficamos
los valores; cada vez que queramos conocer el valor de la ecuación (5.26) ya no
tendremos que realizar la integración, sino que acudiremos al gráfico, entrando a él
únicamente con el periodo fundamental de la estructura. A este grafico se le conoce
como “espectro de respuesta de seudo-velocidad” y a partir de él se puede construir los
espectros de respuesta de seudo-aceleración y seudo-desplazamientos.

Por lo tanto los espectros representan, gráficamente, la respuesta de sistemas de 1 GDL

(respuesta medida en términos de desplazamientos, velocidades y aceleraciones)
caracterizados por sus periodos de vibración T; respecto de un sismo determinado y
representado por su acelererograma respectivo.

A continuación, en los gráficos N° 5.15, N° 5.16 y N° 5.17; se muestran espectros típicos

de seudo-velocidades, seudo-desplazamientos y seudo-aceleraciones
respectivamente.
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

Sv

Figura N° 5.15 T

Sd

Figura N° 5.16
T

Sa

Figura N° 5.17

T
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

5.6. Espectros de Fourier

Las Series de Fourier se utilizan para descomponer la aceleración sísmica en

armónicos. De esta manera se tendrá a un sismo representado por un conjunto de
armónicos, cada uno de ellos definidos por su frecuencia, amplitud y ángulo de fase. El
Espectro de Fourier es la representación gráfica de la importancia de cada armónico
dentro del conjunto.

Por lo tanto, la relevancia de un Espectro de Fourier, en el análisis sísmico, está en que

ellos permiten conocer el contenido de frecuencias del sismo en cuestión.

5.7. Respuesta de un sistema continúo a un movimiento sísmico

A continuación determinaremos la ecuación de movimiento de un sistema como el

mostrado en la figura:

X
F (x,y)

m (x)
K (x)
C (x)

y
S
Figura N° 5.18
 DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

El sistema tiene masa, rigidez y amortiguamiento distribuido es su longitud; y está sujeto,

simultáneamente, a un desplazamiento acelerado de su base y a una fuerza dinámica
lateral distribuida en su longitud.

La solución de la ecuación de movimiento nos permitirá conocer la respuesta vibratoria

del sistema a través de los desplazamientos, en el tiempo, de cualquier punto del
sistema: 𝑌(𝑥,𝑡) .

Si asumimos una función de forma 𝜓(𝑥) , invariable en el tiempo, la respuesta del mismo
sistema puede ser expresada por la siguiente función:

𝑌(𝑥,𝑡) = 𝜓(𝑥) 𝑍(𝑡)

Donde Z(t) corresponde al desplazamiento de un punto específico del sistema. Por lo

tanto, lo que se ha hecho es convertir el sistema continuo en un sistema de 1 GDL.

Consideremos un instante de la respuesta del sistema, en el que el desplazamiento del

suelo es S. Por lo tanto el desplazamiento absoluto de cualquier punto del sistema será:

𝑌𝑎𝑏𝑠 (𝑥,𝑡) = 𝑆 + 𝜓(𝑥) 𝑍(𝑡)

Y por lo tanto, la aceleración absoluta de cualquier punto del sistema será:

̈
𝑌𝑎𝑏𝑠(𝑥,𝑡) = 𝑆̈ + 𝜓(𝑥) 𝑍̈(𝑡)

Si tomamos un diferencial de la longitud del sistema, y en él definimos las fuerzas

generadas, se tiene que:

δw = δu (5.29)

FEXT
F (x,t) dx FI FD

FE dx

Figura N° 5.19
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

Fuerzas externas = 𝐹(𝑥,𝑡) 𝑑𝑥

̈
Fuerza de inercia = 𝑚(𝑥) 𝑑𝑥 𝑌𝑎𝑏𝑠(𝑥,𝑡) = 𝑚(𝑥) (𝜓(𝑥) 𝑍̈(𝑡) + 𝑆̈)𝑑𝑥
̇
Fuerza de amortiguamiento = 𝐶(𝑥) 𝑑𝑥 𝑌(𝑥,𝑡) = 𝐶(𝑥) 𝜓(𝑥) 𝑍̇(𝑡) 𝑑𝑥
Fuerza restitutiva = 𝐾(𝑥) 𝑑𝑥 𝑌(𝑥,𝑡) = 𝐾(𝑥) 𝜓(𝑥) 𝑍(𝑡)

Se sabe que en cualquier instante de tiempo, se cumple:

W=U

Si se da un desplazamiento virtual δy y calculamos la variación tanto del trabajo externo

como de la energía de deformación interna; se tiene que:

δw = δu

Si 𝛿𝑦 = 𝜓(𝑥) 𝛿𝑧, entonces:

El trabajo virtual de las fuerzas de inercia es:

𝛿𝑤𝑝 = ∫ 𝐹(𝑥,𝑡) 𝑑𝑥𝜓(𝑥) 𝛿𝑧 = 𝛿𝑧 ∫ 𝐹(𝑥,𝑡) 𝜓(𝑥) 𝑑𝑥 (5.30)

El trabajo virtual de las fuerzas externas es:

𝛿𝑤𝑖 = ∫ 𝑓𝐼 𝜓(𝑥) 𝛿𝑧 = ∫ 𝑚(𝑥) (𝜓(𝑥) 𝑍̈(𝑡) + 𝑆̈)𝑑𝑥 𝜓(𝑥) 𝛿𝑧

2
= 𝛿𝑧 (∫ 𝑚(𝑥) 𝜓(𝑥) 𝑑𝑥) 𝑍̈(𝑡) + 𝛿𝑧 (∫ 𝑚(𝑥) 𝜓(𝑥) 𝑑𝑥) 𝑆̈(𝑡) (5.31)

Donde:

2
Masa generalizada 𝑀∗ = ∫ 𝑚(𝑥) 𝜓(𝑥) 𝑑𝑥

Masa participante 𝐿∗ = ∫ 𝑚(𝑥) 𝜓(𝑥) 𝑑𝑥

Por lo tanto, la ecuación (5.31) puede escribirse así:

𝛿𝑤𝑖 = 𝛿𝑧𝑀∗ 𝑍̈(𝑡) + 𝛿𝑧𝐿∗ 𝑆̈ (5.32)

El trabajo virtual de las fuerzas de amortiguamiento es:

2
𝛿𝑤𝑐 = ∫ 𝑓𝑐 𝜓(𝑥) 𝛿𝑧 = ∫ 𝐶(𝑥) 𝜓(𝑥) 𝑍̇(𝑡) 𝑑𝑥𝜓(𝑥) 𝛿𝑧 = 𝛿𝑧(∫ 𝐶(𝑥) 𝜓(𝑥) 𝑑𝑥)𝑍̇(𝑡) (5.33)
Donde:
2
Amortiguamiento generalizado 𝐶 ∗ = ∫ 𝐶(𝑥) 𝜓(𝑥) 𝑑𝑥
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

Por lo tanto, la ecuación (5.33) puede escribirse así:

𝛿𝑤𝑐 = 𝛿𝑧𝐶 ∗ 𝑍̇(𝑡) (5.34)

El trabajo virtual de las fuerzas restitutivas es:

2
𝛿𝑤𝑘 = ∫ 𝑓𝑒 𝜓(𝑥) 𝛿𝑧 = ∫ 𝐾(𝑥) 𝑑𝑥𝜓(𝑥) 𝑍(𝑡) 𝜓(𝑥) 𝛿𝑧 = 𝛿𝑧(∫ 𝐾(𝑥) 𝜓(𝑥) 𝑑𝑥)𝑍 (5.35)

Donde:

2
Rigidez externa generalizada 𝐾𝑒∗ = ∫ 𝐾(𝑥) 𝜓(𝑥) 𝑑𝑥

Por lo tanto, la ecuación (5.35) puede escribirse así:

𝛿𝑤𝑘 = 𝛿𝑧 𝐾𝑒∗ 𝑍 (5.36)

Las fuerzas internas de una estructura las medimos en términos de: momentos flectores,
fuerzas cortantes, fuerzas axiales y momentos torsores. En este caso en particular
consideraremos que solo existen momentos flectores. Por lo tanto el trabajo de estos
momentos sobre la rotación causada por el δy nos dará la energía interna de
deformación δu:

1 𝑑 2 (𝛿𝑦(𝑥,𝑡) ) 𝑑 2 (𝜓(𝑥) 𝛿𝑧(𝑡) ) ′′

= = = 𝜓(𝑥) 𝛿𝑧 (5.37)
𝜌 𝑑𝑋 2 𝑑𝑋 2

1 𝑀 ′′
Sabemos que: 𝜃=𝜌 𝑥 y 𝐸𝐼
= 𝜓(𝑥)

′′
Entonces: 𝜃 = (𝜓(𝑥) 𝛿𝑧) 𝑍

Por lo tanto la variación de la energía interna de deformación será:

′′
𝛿𝑈 = ∫ 𝑀(𝑥) 𝜓(𝑥) 𝛿𝑧𝑍𝑑𝑥

′′ ′′ ′′2
= ∫ 𝐸𝐼(𝑥) 𝜓(𝑥) 𝑍(𝑡) 𝑑𝑥𝜓(𝑥) 𝛿𝑧 = 𝛿𝑧(∫ 𝐸𝐼(𝑥) 𝜓(𝑥) 𝑑𝑥)𝑍(𝑡) (5.38)

Donde:
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

′′2
Rigidez generalizada 𝐾𝑢∗ = ∫ 𝐸𝐼(𝑥) 𝜓(𝑥) 𝑑𝑥

Por lo tanto, la ecuación (5.38) puede escribirse así:

𝛿𝑈 = 𝛿𝑍 𝐾𝑈∗ 𝑍(𝑡) (5.39)

Reemplazando las ecuaciones (5.32), (5.34), (5.36) y (5.39) en (5.29), se tiene:

𝑀∗ 𝑍̈(𝑡) + 𝐶 ∗ 𝑍̇(𝑡) + 𝐾 ∗ 𝑍(𝑡) = 𝑃∗ − 𝐿∗ 𝑆̈ (5.40)

La expresión (5.40) corresponde a la ecuación de movimiento de un sistema de 1 GDL,

de masa generalizada M*, amortiguamiento generalizado C*, rigidez generalizada K*
(𝐾𝑒∗ + 𝐾𝑢∗ ) y carga externa generalizada 𝑃∗ − 𝐿∗ 𝑆̈.

Para el caso general en que se tenga una estructura con propiedades distribuidas y
discretas:

2
Carga generalizada 𝑃∗ = ∫ 𝐹(𝑥,𝑡) 𝜓(𝑥) 𝑑𝑥 + ∑ 𝐹𝑖 𝜓(𝑥𝑖)
2 2
Masa generalizada 𝑀∗ = ∫ 𝑚(𝑥) 𝜓(𝑥) 𝑑𝑥 + ∑ 𝑀𝑖 𝜓(𝑥𝑖)
2 2
Amortiguamiento generalizado 𝐶 ∗ = ∫ 𝐶(𝑥) 𝜓(𝑥) 𝑑𝑥 + ∑ 𝐶𝑖 𝜓(𝑥𝑖)
′′ 2 2 2
Rigidez generalizada 𝐾 ∗ = ∫ 𝐸𝐼(𝑥) 𝜓(𝑥) 𝑑𝑥 + ∫ 𝐾(𝑥) 𝜓(𝑥) ∑ 𝐾𝑖 𝜓(𝑥𝑖)
Masa participante 𝐿∗ = ∫ 𝑚(𝑥) 𝜓(𝑥) 𝑑𝑥 + ∑ 𝑀𝑖 𝜓(𝑥𝑖)
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

Cap VI. Análisis Sísmico de edificios

Introducción

Hasta el momento, lo que se tiene es:

1.- Un sismo de diseño expresado a través de una aceleración máxima del suelo, que
tiene una probabilidad de excedencia de 10% en 50 años de exposición.
2.- Un modelo de análisis en el que está representada la masa, rigidez y
amortiguamiento del edificio real; y cuya cinemática se define a través de 3 GDL por
piso (siempre y cuando los techos se comporten como diafragmas rígidos).

Lo que queda ahora es TRADUCIR la aceleración del suelo en fuerzas de inercia

horizontales actuando sobre la masa discreta del modelo; y así poder determinar las
demandas de desplazamientos y de resistencia generada por el sismo de diseño. A esta
metodología se le conoce como “diseño basado en fuerzas”, y es la que desarrollaremos
en el curso.

P3(t)

P2(t)

P1(t)

Figura N° 6.1
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

6.1. Respuesta de un sistema discreto de varios GDL sujeto a un

movimiento horizontal acelerado en su base

Consideremos un modelo plano de masas concentradas y rigideces laterales de

entrepiso Ki sujeto a un sistema de cargas laterales variables en el tiempo:

M3
𝑃3𝑒𝑞(𝑡)
K3

M2
𝑃2𝑒𝑞(𝑡)
K2

M1
𝑃1𝑒𝑞(𝑡)

K1

Figura N° 6.2

Para plantear las ecuaciones de equilibrio dinámico, asumiremos una configuración

deformada 𝜓(𝑥𝑖) , así:

Ψ3=X3

Ψ2=X2

Ψ1=X1

Figura N° 6.3
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

Ecuación de equilibrio sobre el GDL 1:

𝐾2 (𝑋2 − 𝑋1)

𝑀1 𝑋̈1 𝑃1𝑒𝑞

𝐾1 𝑋1

𝑀1 𝑋̈1 + 𝐾1 𝑋1 − 𝐾2 (𝑋2 − 𝑋1 ) = 𝑃1𝑒𝑞

𝑀1 𝑋̈1 + 𝐾1 𝑋1 − 𝐾2 𝑋2 + 𝐾2 𝑋1 = 𝑃1𝑒𝑞
𝑀1 𝑋̈1 + 𝑋1 (𝐾1 + 𝐾2 ) − 𝑋2 𝐾2 = 𝑃1𝑒𝑞

Ecuación de equilibrio sobre el GDL 2:

𝐾3 (𝑋3 − 𝑋2)

𝑀2 𝑋̈2
𝑃2𝑒𝑞

𝐾2 (𝑋2 − 𝑋1)

𝑀2 𝑋̈2 + 𝐾2 (𝑋2 − 𝑋1 ) − 𝐾3 (𝑋3 − 𝑋2 ) = 𝑃2𝑒𝑞

𝑀2 𝑋̈2 + 𝐾2 𝑋2 − 𝐾2 𝑋1 − 𝐾3 𝑋3 + 𝐾3 𝑋2 = 𝑃2𝑒𝑞
𝑀2 𝑋̈2 − 𝐾2 𝑋1 + 𝑋2 (𝐾2 + 𝐾3 ) − 𝐾3 𝑋3 = 𝑃2𝑒𝑞

Ecuación de equilibrio sobre el GDL 3:

𝑀3 𝑋̈3
𝑃3𝑒𝑞

𝐾3 (𝑋3 − 𝑋2)

𝑀3 𝑋̈3 + 𝐾3 (𝑋3 − 𝑋2 ) = 𝑃3𝑒𝑞

𝑀3 𝑋̈3 − 𝐾3 𝑋2 − 𝐾2 𝑋1 + 𝐾3 𝑋3 = 𝑃3𝑒𝑞
𝑀1 𝑋̈1 + 𝑋1 (𝐾1 + 𝐾2 ) − 𝑋2 𝐾2 = 𝑃1𝑒𝑞
𝑀2 𝑋̈2 − 𝐾2 𝑋1 + 𝑋2 (𝐾2 + 𝐾3 ) − 𝐾3 𝑋3 = 𝑃2𝑒𝑞
𝑀3 𝑋̈3 − 𝐾3 𝑋2 − 𝐾2 𝑋1 + 𝐾3 𝑋3 = 𝑃3𝑒𝑞
DINAMICA DE ESTRUCTURAS Profesor: JOHN P. ARAGON BROUSSET

Expresando matricialmente este sistema de ecuaciones, se tiene:

[𝑀]{𝑋̈} + [𝐾]{𝑋} = {𝑃𝑒𝑞 }

Dónde:

𝑀1 0 0
[𝑀] = | 0 𝑀2 0|
0 0 𝑀3

𝐾1 + 𝐾2 −𝐾2 0
[𝐾] = | −𝐾2 𝐾2 + 𝐾3 −𝐾3 |
0 −𝐾3 𝐾3

𝑋̈1 𝑋1 𝑃1𝑒𝑞
̈
{𝑋} = {𝑋̈2 } {𝑋} = {𝑋2 } {𝑃𝑒𝑞 } = {𝑃2𝑒𝑞 }
𝑋̈3 𝑋3 𝑃3𝑒𝑞

Sabemos que para un desplazamiento horizontal acelerado Xs en la base, las fuerzas

Peq son iguales a:

𝑃𝑒𝑞 = −𝑀𝑋̈𝑠

Por lo tanto:

−𝑀1 𝑋̈𝑆 −𝑀1

{𝑃𝑒𝑞 } = {−𝑀2 𝑋̈𝑆 } 𝑂 {𝑃𝑒𝑞 } = {−𝑀2 } 𝑋̈𝑆
−𝑀3 𝑋̈𝑆 −𝑀3

Luego, el sistema de ecuaciones puede expresarse así:

[𝑀]{𝑋̈} + [𝐾]{𝑋} = −{𝑀}𝑋̈𝑆 (6.1)

Observando la matriz de masas y la matriz de rigideces, vemos que el sistema de

ecuaciones esta desacoplada dinámicamente, pero esta acoplado estáticamente.