Análisis Sísmico – 2° Semestre 2021

Espectros de respuesta
Espectros de respuesta 1 de 15

ESPECTROS DE RESPUESTA
 Espectros de respuesta 2 de 15

La definición de espectro de respuesta se le atribuye a Housner.

Corresponde a una definición que sirve para caracterizar movimiento del
terreno y con aquello determinar el efecto en las estructuras.
El método para calcularlo corresponde a determinar la respuesta de un
oscilador de 1 GDL. Para un amortiguamiento fijo, se determina la respuesta
máxima para el oscilador para distintas frecuencias (o periodos).

Se pueden determinar espectros de respuesta de:

-Aceleración
-Velocidad
-Deformación
Espectro de deformación
proporciona toda la información
necesaria
-> Deformación máxima
-> Fuerzas internas
 Espectros de respuesta 3 de 15

Espectro de deformación

Fuente: Dinámica de estructuras - Chopra

 Espectros de respuesta 4 de 15

Actividad:
Utilizando MATLAB, calcule el espectro de desplazamiento en base al registros
compatible con DS61 del terremoto 27F para la locación de Constitución.
Considerar Tn hasta los 5 segundos.
Utilice el “método exacto” o “método de diferencia centrada” para la
determinación de la respuesta del oscilador de 1 GDL.
Solución:
Recuerdo dinámica aplicada – Método exacto

Fuente: Dinámica de estructuras - Chopra

 Recuerdo Dinámica

Método exacto:
• Se realiza una interpolación de la excitación en cada instante de tiempo y luego resuelve la
ecuación de movimiento con metodologías convencionales.

La ecuación a resolver es:

Para cada intervalo la respuesta u(t) está compuesta de:
1) respuesta producto de la vibración debida a las condiciones iniciales de desplazamiento y
velocidad u y u respectivamente.
2) Respuesta producto de la fuerza bajo condiciones iniciales nulas
3) Respuesta de carga bajo condiciones iniciales nulas
Recuerdo Dinámica
Recuerdo Dinámica
 Recuerdo Dinámica

Método diferencia centrada:

• Método basado en diferencias finitas
• Si se considera un paso de tiempo constante Dt:

Reemplazando en ecuaciones de mov:

 Espectros de respuesta 5 de 15

clear function [Sde] =Espect_delt_DS61(T,xi)

close all % Función que genera el espectro elástico de
desplazamientos en base al DS61
clc % Para tipo de suelo B, I=1 y zona 3
%________________________________ % Carlos Aguilar
% Rev A 27-08-14
___________ % Los valores que se deben definir son
% Generación de un vector de % T : periodo
T0=0.3; % s Para suelo tipo B
periodos hasta 5 seg p=1.5; %Para suelo tipo B
Bi=0.02; I=1;
G=981; %cm/s2
T=0:0.02:5; %vector de periodos n=length(T);
A0=0.4*G;
%________________________________ for k=1:n
___________
if T(k)<=0.47
% Generación de vector de Cd=1;
desplazamientos alf=((1+4.5*((T(k)/T0)^p))/(1+((T(k)/T0)^3)));
Sde(k)=(0.05/xi)^0.4*((T(k)^2)*alf*A0*Cd)/(4*pi^2);
% y gráfico Periodo v/s elseif T(k)<=2.02
Desplazamiento Cd=(0.95*T(k))+(0.55);
alf=((1+4.5*((T(k)/T0)^p))/(1+((T(k)/T0)^3)));
% en base al DS 61 Sde(k)=(0.05/xi)^0.4*((T(k)^2)*alf*A0*Cd)/(4*pi^2);
[Sde]= Espect_delt_DS61(T,Bi); elseif T(k)<=5
Cd=(0.065*T(k)^2)-(0.75*T(k))+(3.72);
alf=((1+4.5*((T(k)/T0)^p))/(1+((T(k)/T0)^3)));
Sde(k)=(0.05/xi)^0.4*((T(k)^2)*alf*A0*Cd)/(4*pi^2);
end

end
Espectros de respuesta 6 de 15
Espectros de respuesta 7 de 15

nT=length(T);
% Generación de vector de w=zeros(nT,1);
desplazamientos D=zeros(nT,1);
% en base a registro compatible DS61
global S DT N %function
Consti_2010 %Unidades [u,v,t]=Exacta(bi,wi,dt,N,uo,vo,p)
cm/seg^2
DT = dT; S = ux; N = length(S); for i=2:nT
w(i)=2*pi()/T(i);
tt=N*DT; %100;
w2=w(i)^2;
dt=DT; [qi,qdi,t2] =
uo=0; Exacta(Bi,w(i),DT,N,0,0,S');
vo=0; D(i)=max(abs(qi));
uu = zeros(N,1); V(i)=w(i)*D(i);
vv = zeros(N,1); A(i)=w(i)*V(i);
acc = zeros(N,1);
%___________________________________ end

Fuente: Dinámica de estructuras - Chopra

Espectros de respuesta 8 de 15

function
[u,v,t]=Exacta(bi,wi,dt,N,uo,vo
,p) Ap=-f1*(wi/f3*sn);
k=wi^2; Bp=f1*(cs-f2*sn);
WD=wi*sqrt(1-bi^2); Cp=(1/k)*(-
1/dt+f1*((wi/f3+f2/dt)*sn+1/dt*cs));
f1=exp(-1*bi*wi*dt); Dp=1/(k*dt)*(1-f1*(f2*sn+cs));
f2=bi/sqrt(1-bi^2);
f3=sqrt(1-bi^2); %Definición de los vectores del sistema
f4=WD*dt; t = zeros(N,1);
f5=wi*dt; u = zeros(N,1); %Solución ecuación de movimiento
v = zeros(N,1); armonico método integracion
sn=sin(WD*dt); exacta
cs=cos(WD*dt); u(1) = uo; %Variables
v(1) = vo; %bi=factor de amortiguamiento
A=f1*(f2*sn+cs); critico
B=(f1/WD)*sn; for i =1:1:(N-1) %wi=frecuencia del modo i
C=(1/k)*(2*bi/f5+f1*(((1- ui=u(i); %k=coeficientes de rigidez
2*bi^2)/f4-f2)*sn- vi=v(i); %dt=intervalo de tiempo del sismo
(1+2*bi/f5)*cs)); pi=p(i); %ug=aceleración sísmica
D=(1/k)*(1-2*bi/f5+f1*((2*bi^2- pi1=p(i+1); %
1)/f4*sn+2*bi/f5*cs)); %A,B,C,D coeficientes del metodo
u(i+1)=A*ui+B*vi+C*pi+D*pi1; %Ap, Bp, Cp, Dp, coeficientes del
v(i+1)=Ap*ui+Bp*vi+Cp*pi+Dp*pi1; metodo

t(i+1) = t(i) + dt;

end
Espectros de respuesta 9 de 15
 Espectros de respuesta 10 de 15

Espectro de pseudo-velocidad

Para un sistema de 1 GDL con

frecuencia natural ω𝑛 , si se 2𝜋
V= ω𝑛 𝐷 = 𝐷
tiene que D es la deformación 𝑇𝑛
máxima, se puede considerar V:

La energía de deformación
𝑚𝑉 2
almacenada en el sistema 𝐸𝑆𝑜 =
durante el sismo es: 2

V: Pseudo – velocidad 𝑢ሶ ≠ 𝑉
m: Masa estructural

Espectro de Pseudo velocidad se genera a

partir del espectro de desplazamiento
 Espectros de respuesta 11 de 15

Espectro de pseudo-velocidad 2𝜋
V= ω𝑛 𝐷 = 𝐷
𝑇𝑛 Tn = 0.5 s

D = 2.67 in

2  in
V =  D = 33.55
Tn s

Tn = 1 s

D = 5.97 in

2  in
V =  D = 37.51
Tn s

Tn = 2 s

D = 7.47 in

2  in
V =  D = 23.47
Tn s
 Espectros de respuesta 12 de 15

Espectro de pseudo-velocidad
Tn = 0.5 s

D = 2.67 in

2  in
V =  D = 33.55
Tn s

Tn = 1 s

D = 5.97 in

2  in
V =  D = 37.51
Tn s

Tn = 2 s

D = 7.47 in

2  in
V =  D = 23.47
Tn s
 Espectros de respuesta 13 de 15

Espectro de pseudo-aceleración

Para un sistema de 1 GDL con

frecuencia natural ω𝑛 , si se 2𝜋
A= ω𝑛 2 𝐷 = ( 𝑇 )2 𝐷
tiene que D es la deformación 𝑛
máxima, se puede considerar A:

A posee unidades de
𝑉𝑏𝑜 = 𝑓𝑠𝑜 = 𝑚A 𝑢ሷ ≠ 𝐴
aceleración y se relaciona con el
corte basal como:

El corte basal máximo se puede 𝐴

escribir como: 𝑉𝑏𝑜 = 𝑤 𝑤: 𝑃𝑒𝑠𝑜 𝑑𝑒 𝑙𝑎 𝑒𝑠𝑡𝑟𝑢𝑐𝑡𝑢𝑟𝑎
𝑔
𝐴
:Coeficiente sísmico
𝑔

Espectro de Pseudo aceleración se genera a

partir del espectro de desplazamiento
 Espectros de respuesta 14 de 15

Espectro de pseudo-aceleración

Tn = 0.5 s

D = 2.67 in
2
2   D
A =   = 1.09
Tn  g
 

Tn = 1 s

D = 5.97 in
2
A =  2    D = 0.61
T  g
 n

Tn = 2 s

D = 7.47 in
2
A = 
2   D
 = 0.191
Tn  g
 
 Espectros de respuesta 15 de 15

Espectro combinado (Newmark)

Espectro combinado El Centro (2%) - Chopra