Antisismia Trabajo Informe - Matlab

UNIVERSIDAD DE HUANUCO
FACULTAD DE INGENIERÍA
ESCUELA ACADEMICA PROFESIONAL:
INGENIERÍA CIVIL
DOCENTE: ING. ERICK REM LOVON DAVILA
CURSO: INGENIERIA ANTISISMICA
SECCION: “B”
CICLO: IX
INTEGRANTES:
 MATEO VERASTEGUI STIVEN
HUANUCO- PERU
2018
ANTISISMICA
Código 141509021 - INGENIERÍA ANTISÍSMICA

Trabajo 2 - Parte A: Espectro Elástico y Matriz de Rigidez
Fecha de Entrega Parte A y B: 26 de noviembre del 2018 a las 9:00 pm
Ejercicio 1
Utilizar el registro de Viña del Mar Centro del terremoto del Maule de 2010
proporcionado en clases. Considere la componente NS del registro.
a. Grafique la aceleración versus tiempo e indique la aceleración máxima del
suelo (PGA).
b. Integre el registro para obtener la velocidad del suelo. Grafique la velocidad
versus tiempo del suelo e indique la velocidad máxima del suelo (PGV). Utilice
las condiciones iniciales que se indican en el registro.
c. Integre el registro para obtener el desplazamiento máximo del suelo. Grafique
el desplazamiento versus tiempo del suelo e indique el desplazamiento máximo
del suelo (PGD). Utilice las condiciones iniciales que se indican en el registro.
d. Utilizando el método de integración exacto (función del taller 1), calcule el
desplazamiento de un sistema de un grado de libertad. Considere un periodo Tn
= 1 seg y un amortiguamiento de 5%.
SOLUCIÓN EN MATLAB
SISMO 1
clc;
clear;
disp('***********ING. ANTISISMICA********')
disp('******* MATEO VERASTEGUI STIVEN **********')
disp('***GRUPO: “B”****')
%numero de datos
n=12502;
DEFINIMOS LA FUNCIÓN
a=[1.694 ……….. 7.122];
INTDIFCEN INDICANDO EL
tt=0.01; IMPUT Y OUTPUT Y ADEMAS
%************* DEFINIMOS LOS PARAMETROS
t=zeros(1,n);
t(1,1)=0;
d=zeros(1,n);
v=zeros(1,n);
d(1,1)=input('Ingrese Posición Inicial: -0.009 ');
v(1,1)=input('Ingrese Velosidad Inicial: -0.053 ');
for i=1:n-1
t(1,i+1)=t(1,i)+tt;
v(1,i+1)=(a(1,i)*tt/2)+(a(1,i+1)*tt/2)+v(1,i);
d(1,i+1)=(a(1,i)*(tt^2)/3)+(a(1,i+1)*(tt^2)/6)+(v(1,i)
*tt)+d(1,i);
ANTISISMICA
end
disp(t);
disp(d);
disp(v);
vm=0;
for i=1:n
if abs(v(1,i))>vm;
vm=abs(v(1,i));
end
end
dm=0;
for i=1:n
if abs(d(1,i))>dm;
dm=abs(d(1,i));
end
end
disp('Velosidad Máxima');
disp(vm);
disp('Desplazamiento Máxima');
disp(dm);
subplot(3,1,1);
plot(t,a,'b');
title('Grafica de Tiempo vs. Aceleración, En el
Suelo')
xlabel('Tiempo (Seg)')
ylabel('Aceleración (cm/s2)')
subplot(3,1,2);
plot(t,v,'r');
title('Grafica de Tiempo vs. Velosidad, En el Suelo')
ylabel('Velosidad (cm/s)')
subplot(3,1,3);
plot(t,d,'m');
title('Grafica de Tiempo vs. Desplazamiento, En el
Suelo')
ANTISISMICA
ylabel('Dezplazamiento (cm)')
OBTENEMOS LOS VALORES DE LOS PARAMETROS Y LA SOLUCION
ANTISISMICA
METODO EXACTO CON INTERPOLACIÓN
GRAFICO RESULTANTE
Luego de haber definido todos los parámetros esta es la curva que se obtuvo
SISMO 2
n=12502;
pto=[1.694 1.779 3.04 4.48…………
%intervalo de tiempo
tt=0.01;
%termino de edicion
%*****************************************************
*********************
m=input('Ingrese la masa: ');
e=input('Ingrese razón de amortiguamiento: ');

t=input('Ingrese el periodo: ');
wn=2*pi/t;
ANTISISMICA
to=input('Ingrese Tiempo Inicial: ');

do=input('Ingrese Desplazamiento Inicial:-0.009 ');
vo=input('Ingrese Velosidad Inicial:-0.053 ');
u=zeros(1,n);
uu=zeros(1,n);
u(1,1)=do;
uu(1,1)=vo;
ttt=zeros(1,n);
ttt(1,1)=to;
k=m*(wn^2);
wd=wn*(1-e^2)^0.5;
x=sin(wd*tt);
y=cos(wd*tt);
A1=(wn*e)/wd;
ee=exp(-e*wn*tt);
C1=(2*e)/(wn*tt);
C2=(1-2*(e^2))/(wd*tt);
C3=wn/((1-e^2)^0.5);
A=ee*(A1*x+y);
B=ee*x/wd;
C=(C2-A1)*x-(1+C1)*y;
C=(1/k)*(C1+ee*C);
D=(-C2*x+C1*y);
D=(1/k)*(1-C1+ee*D);
AA=-ee*C3*x;
BB=ee*(y-A1*x);
CC=((C3+(A1/tt))*x)+(y/tt);
CC=(1/k)*(-(1/tt)+ee*CC);
DD=A1*x+y;
DD=(1/(k*tt))*(1-ee*DD);
disp('Coeficiente A sera:');
disp(A);
disp('Coeficiente B sera:');
disp(B);
disp('Coeficiente C sera:');
disp(C);
disp('Coeficiente D sera:');
disp(D);
disp('Coeficiente A'' sera:');
disp(AA);
disp('Coeficiente B'' sera:');
disp(BB);
ANTISISMICA
disp('Coeficiente C'' sera:');

disp(CC);
disp('Coeficiente D'' sera:');
disp(DD);
for i=1:n-1
u(1,i+1)=A*u(1,i)+B*uu(1,i)+C*pto(1,i)+D*pto(1,i+1);
uu(1,i+1)=AA*u(1,i)+BB*uu(1,i)+CC*pto(1,i)+DD*pto(1,i+
1);
ttt(1,i+1)=ttt(1,i)+tt;
end
sd=0;
for i=1:n
if abs(u(1,i))>sd;
sd=abs(u(1,i));
end
end
sv=wn*sd;
sa=(wn^2)*sd;
disp('Tiempo:');
disp(ttt);
disp('Desplazamiento:');
disp(u);
disp('Velosidad:');
disp(uu);
disp('Periodo');
disp(t);
disp('Espectro Desplazamiento Máximo')
disp(sd);
disp('Espectro Pseudo-Velosidad')
disp(sv);
disp('Espectro Pseudo-Aceleración')
disp(sa);
%subplot(1,2,1);
plot(ttt,u,'r');
ANTISISMICA
title('Grafica Dezaplamiento vs Tiempo, En la

Estructura');
xlabel('Tiempo (Seg.)')
ylabel('Dezplamiento (Cm.)')
%subplot(1,2,2);
%plot(ttt,uu,'b');
%title('Grafica Velosidad vs Tiempo, En la
Estructura');
%xlabel('Tiempo (Seg.)')
%ylabel('Velosidad (Cm/Seg)')
clear all; clc;

disp('**********PREGUNTA N 1**********')
disp('CURSO: INGENIERIA ANTISISMICA')
disp('DOCENTE: Ing. Eric Lovon Davila')
disp('ALUMNO: MARCO CONDEZO GONZALES')
disp('grupo: A')
%% DEFINICION DE LAS PROPIEDADES DEL SISTEMA
m=1; %kg Masa

Tn=1; %seg Periodo no Amortiguado
wn=2*pi/Tn; %rad/s Frecuencia Angular
k=wn^2*m; % Rigidez del Sistema
DEFINIMOS LAS
po=10; %N PROPIEDADESDEL
td=0.5; %seg Periodo Amortiguado SISTEMA,
r=0.05; DESPLAZAMIENTO
c=r*2*m*wn; ESTATICO Y
wd=wn*sqrt(1-r^2);
% Definicion del Desplazamiento Estatico VECTOR TIEMPO
ANTISISMICA
ust=po/k
% Definicion del Vector Tiempo
t=0:0.1:1;
DEFINIMOS EL
% Definicion del Desplazamiento DESPLAZAMIENTO
PARA LOS VALORES Ti
Y PROCEDEMOS A
for i=1:length(t) GRAFICAR
ti=t(i);
ut(i)=ust/(2*r)*(exp(-r*wn*ti)*(cos(wd*ti)+r/sqrt(1-
r^2)*sin(wd*ti))-cos(wn*ti));
end
t2=0:0.001:1
for i=1:length(t2)
ti=t2(i);
ut2(i)=ust/(2*r)*(exp(-r*wn*ti)*(cos(wd*ti)+r/sqrt(1-
r^2)*sin(wd*ti))-cos(wn*ti));
end
%Graficando los resultados
figure
plot(t,ut,'-o');hold on
plot(t2,ut2,'r');hold on
METODO EXACTO
CON
INTERPOLACIÓN
GRAFICO
RESULTANTE
ANTISISMICA
EJERCICIO Nº 3
function[Sd,Sv,Sa,Svr]=
CalculoEspectro(ag,dtag,nu,T,uo,vo)
%UNTITLED3 Summary of this function goes here
% Detailed explanation goes here
m=1;
n=length(ag);
t=[0:dtag:(n-1)*dtag]';
for i=1:length(T)
if T(i)<0.2
dt(i)=T(i)/20;
else
dt(i)=0.01;
end
t2=[0:dt:t(end)]';
%interpolacion de rspuestas
ag2=interpl(t,ag,t2);
wn(i)=2*pi/T(i); % frecuencia angular
k(i)=wn(i)^2 * m;
c(i)=2*m*nu*wn(i);
pt2=m*ag2/100; % metros
%Integracion de la respuesta
[]=IntExacto()
ANTISISMICA
end
end

Grafica 1
ANTISISMICA
EJERCICIO Nº 4
function Kc = RigidezColumna(L,E,I,A)
P=12*E*I/L^3;
Q=6*E*I/L^2;
R=4*E*I/L;
S=A*E/L;
Kc=[P 0 -Q -P 0 -Q
0 S 0 0 -S 0
-Q 0 R Q 0 R/2
-P 0 Q P 0 Q
0 -S 0 0 S 0
-Q 0 R/2 Q 0 R];
end
function Kv = RigidezViga(L,E,I,A)
P=12*E*I/L^3;
Q=6*E*I/L^2;
R=4*E*I/L;
Kv=[P Q -P Q
Q R -Q R/2
-P -Q P -Q
Q R/2 -Q R];
end
function K = ColocarK(K,kelem,tabla,indice)
tablai=tabla(indice,:); %lee los grados de
libertad
i=find(tablai~=0); %indices rigidos del
elemento
tablai=tablai(tablai~=0); %
K(tablai,tablai)=K(tablai,tablai)+kelem(i,i);
ANTISISMICA
EJERCICIO Nº 5
close all; clear all; clc;
disp('******ING. ANTISISMICA******')
disp('******MATEO VERASTEGUI STIVEN ******')
disp('******GRUPO:"B" ******')
disp('******ELASTICIDAD Y RIGIDEZ******')
E=20000*10^6;
Lc=2.8;
Ic=(0.5*0.5^3)/12;
Ac=0.5*0.5;
ANTISISMICA
%consideraciones para la viga
Lv=7;
Iv=(0.3*0.6^3)/12;
Kc=RigidezColumna(Lc,E,Ic,Ac);
Kv=RigidezViga(Lv,E,Iv);
TablaC=[0 0 0 25 1 2
0 0 0 25 3 4
25 1 2 26 5 6
25 3 4 26 7 8
26 5 6 27 9 10
26 7 8 27 11 12
27 9 10 28 13 14
27 11 12 28 15 16
28 13 14 29 17 18
28 15 16 29 19 20
29 17 18 30 21 22
29 19 20 30 23 24];
TablaV=[1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16
17 18 19 20
21 22 23 24];
K=zeros(30,30);
for i=1:12
K=ColocarK(K,Kc,TablaC,i);
end
for i=1:6
K=ColocarK(K,Kv,TablaV,i);
end
%aplicar el vector de fuerzas externas del sistema;
F=zeros(30,1);
F(30,1)=10000; %vector de fuerzas en Newton (N)
u=K\F; %desplazamiento del sistema
%Desplazamiento lateral
ANTISISMICA
DespLat=u(25:30);
Fcol1=Kc*[0;0;0;u(25);u(1);u(2)];
Fcol2=Kc*[0;0;0;u(25);u(1);u(2)];
K11=K(1:24,1:24);
K12=K(1:24,25:30);
K21=K12';
K22=K(25:25:30);
Kd=inv(K11)*K12;
Kcon=K22-K21*Kd

GRAFICO
ANTISISMICA

Antisismia Trabajo Informe - Matlab

Cargado por

Copyright:

Formatos disponibles

Antisismia Trabajo Informe - Matlab

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Antisismia Trabajo Informe - Matlab

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD DE HUANUCO

DOCENTE: ING. ERICK REM LOVON DAVILA

CURSO: INGENIERIA ANTISISMICA

 MATEO VERASTEGUI STIVEN

Código 141509021 - INGENIERÍA ANTISÍSMICA

OBTENEMOS LOS VALORES DE LOS PARAMETROS Y LA SOLUCION

METODO EXACTO CON INTERPOLACIÓN

e=input('Ingrese razón de amortiguamiento: ');

to=input('Ingrese Tiempo Inicial: ');

disp('Coeficiente C'' sera:');

title('Grafica Dezaplamiento vs Tiempo, En la

clear all; clc;

%% DEFINICION DE LAS PROPIEDADES DEL SISTEMA

m=1; %kg Masa

% Definicion del Vector Tiempo

%Graficando los resultados

METODO EXACTO CON INTERPOLACIÓN

close all; clear all; clc;

%consideraciones para la viga

%aplicar el vector de fuerzas externas del sistema;

METODO EXACTO CON INTERPOLACIÓN

También podría gustarte