Análisis espectral
El análisis espectral es determinar las acciones sísmicas que pueden actuar sobre una estructura de acuerdo
al espectro de diseño de un ciudad. El espectro de diseño de una localidad obedece a un estudio estadístico
que toma en cuenta el historial sísmico, las fallas geológicas con potencial sísmico, el tipo de suelo, entre
otros, para determinar un espectro en función del periodo de la estructura y aceleración del suelo.
figura1
La figura 1 muestra las caracteristicas del espectro de diseño del Ecuador de acuerdo a la NEC-11. En el
eje de las x se tiene el periodo y en el eje de lasy las aceleraciones en un suelo tipo roca (PGA). El espectro
se determina a travéz de 3 funciones, una función lineal, una función constante y una función no lineal. Los
parámetros que definen el espectro están en función de la region a la que pertenece la localidad en estudio y al
tipo de suelo.
1
2
Tipos de suelo
3
clear
load modelo
load KM_Phi_Omega
T=2*pi()./omega;% Periodos de la estructura
Modelo
figure
plotnodes(Nodes,'Numbering','off')
hold on
plotelem(Nodes,Elements,Types,'Numbering','off')
4
Definición del espectro para Cuenca
z=0.25;%[g]
eta=2.48;%[-]
En Cuenca es común el tipo de suelo C
fa=1.3;%[-]
fd=1.28;%[-]
fs=0.94;%[-]
r=1;%[-]
[Sa_,T_]=ecuspec(fa,fd,fs,eta,r,z); % peridoos y aceleraciones espectrales
%
Sa=interp1(T_,Sa_,T); % Se desea interpolar para conocer las acelraciones que afectan a la estr
figure
plot(T_,Sa_,T,Sa,'ro')
% plot(T_,Sa_)
xlabel('Periodo [s]')
ylabel('Sa [g]')
5
El espectro de que se muestra representa la acelaración que el sismo de diseño le puede causar a un
estructura particular localizada en el lugar que define al espectro. Es importante notar que el espectro define
acelaraciones para un rango de periodos de 0 a 10 s.
Análisis estático
DLoads=accel([0 0 9.81],Elements,Types,Sections,Materials);
W=elemloads(DLoads,Nodes,Elements,Types,DOF);
Ust=K\W;
% figure
% plotdisp(Nodes,Elements,Types,DOF,Ust,'DispScal',1000)
Análisis espectral
Influencia del sismo sobre la estructura
En la direccion y, se seleccionan los grados de libertad que correspondan a y(2)
Egy=sum(selectdof(DOF,0.02),1).'; % Se selecionan los grados de libertad en la dirección y
En la dirección x se selecciona 1
Egx=sum(selectdof(DOF,0.01),1).'; % Se selecionan los grados de libertad en la dirección x
Participación durante la acción sísmica
6
Distribucion de masa en la estructura
sy=M*Egy;% en la dirección y
sx=M*Egx;% en la dirección x
Transformación a fuerzas modales que dependen del factor de distribución Γ. Se obtiene del vector de
fuerza en función de la importancia de cada modo de vibrar
Gammay=phi.'*sy;% Gamma y
Gammax=phi.'*sx;% Gamma x
Masa participante modal (por modo) en las direcciones x y y
Masty=(Gammay.*Gammay);
Mastx=(Gammax.*Gammax);
% Mpar=Mast/Mtot*100;
Vector de fuerza efectiva (el que actúa como representante del sismo sobre la estrcutura)
sjy=M*phi*diag(Gammay);
sjx=M*phi*diag(Gammax);
Desplazamientos causados por sj de manera estática
usjy=K\sjy;
usjx=K\sjx;
Desplazamientos causados por el sismo (dinámica) y que estan en funcion de sj
Udy=usjy*Sa*9.81;
Uy=Udy+Ust;
Udx=usjx*Sa*9.81;
Ux=Udx+Ust;
Deformaciones en X y en Y
%plotdisp(Nodes,Elements,Types,DOF,Ux,'DispScal',100)
%plotdisp(Nodes,Elements,Types,DOF,Uy)
Determinación de las derivas para asesorar la capacidad sísmica de la edificación.
Las derivas es la diferencia entre la deformación dos puntos. En estructuras interesa conocer la deriva entre
plantas que se define por:
en donde es la deformación máxima del piso i+1, es la deformación máxima del piso i y H es la altura
de entrepiso. Por ejemplo determinar la deriva entre la planta baja y la primera planta del modelo mostrado en
7
las direcciones x y y, determinar la dirección crítica del edificio y calcular la deriva en la dirección crítica entre la
primera planata ata y el novel donde empieza la cubierta.
• Se sabe que las deformaciones de todos los nodos en las direcciones x y y son 0 en la planta baja por la
que ya se conoce , tambien conocemos H
• Se determina el nodo que tiene la mayor deformación absoluta en el entrepiso en la primera planta alta
Deriva en la dirección x nivel 0 hasta 2.7 m
[Nodsel,Selem]=selectelem(Nodes,Elements,Types,[0 0 2.7 100 100, 2.70]);% Seleción de nodos y e
plotnodes(Nodsel,'Numbering','off')
plotdisp(Nodsel,Elements(Selem,:),Types,DOF,Ux,'DispScal',100)
8
De nodos seleccionados se requieren unicacmente los grados de libertad en x
seldofx=Nodsel(:,1)+0.01;% se selecionan solo los grados de libertad en x
seldofx=selectdof(DOF,seldofx);% se determinan la posición de esos grados de libertad
Ux_s=seldofx*Ux;% se multiplica la matriz de grados de libertad seleccionados con vector de def
[d_i1x,~]=max(abs(Ux_s));% deformación máxima del entrepiso i+1 (primera planta alta)
disp(num2str(d_i1x,5))
0.004788
Se define la deriva
d_rx1=(d_i1x-0)/2.70;
disp(num2str(d_rx1,5))
0.0017733
Deriva en la dirección y
figure
plotdisp(Nodsel,Elements(Selem,:),Types,DOF,Uy,'DispScal',100)
9
seldofy=Nodsel(:,1)+0.02;% se selecionan solo los grados de libertad en x
seldofy=selectdof(DOF,seldofy);% se determinan la posición de esos grados de libertad
Uy_s=seldofx*Uy;% se multiplica la matriz de grados de libertad seleccionados con vector de def
[d_i1y,~]=max(abs(Uy_s));% deformación del entrepiso i+1
disp(num2str(d_i1y,5))
0.00023014
Se define la deriva
d_ry=(d_i1y-0)/2.70;
disp(num2str(d_ry,5))
8.5238e-05
Determinar la deriva entre la primera planta y el nivel de cubierta.
La dirección crítica es la x
[Nodsel,Selem]=selectelem(Nodes,Elements,Types,[0 0 5.4 100 100 5.4]);
seldofx=Nodsel(:,1)+0.01;% se selecionan solo los grados de libertad en x
seldofx_=selectdof(DOF,seldofx);% se determinan la posición de esos grados de libertad
Ux_s2=seldofx_*Ux;% se multiplica la matriz de grados de libertad seleccionados con vector de d
[d_i2x,pos]=max(abs(Ux_s2));% deformación en nivel de cubierta
disp(num2str(d_i2x,5))
0.028872
10
El nodo que más se deforma es
seldofx(pos)
ans =
499.01
figure
plotnodes(Nodsel)
Se define la deriva
d_rx2=(d_i2x-d_i1x)/2.70;
disp(num2str(d_rx2,5))
0.0089201
Asesoramiento: Si la deriva en la dirección x o y es mayor a 0.006 el edificio tiene peligro de colapso. Si la
deriva es mayor a 0.003 pero menor a 0.006 el edificio puede presentar daños y si la deriva es menor 0.003
el edifcio no sufre daños considerables. En este caso la edificación muestra peligro de colpaso en la dirección
x en los niveles primera planta a cubierta, pues la deriva es mayor a 0.006 (6 por mil). En el análisis modal se
pudo ver que la pared ubicada en...... presentaba un movimiento mayor con respecto al resto de paredes. Se
esperaba que esta pared sea afectada durante la acción sismica.
De las deformaciones totales en x y en y se puede concluir que el área más peligrosa es...... y la menos
peligrosa es... se recomienda que durante una acción sísmica si no se puede salir la familia se ubique
alrededor de....
11
12
