% Ejercicio No1: Vibracion libre
% Se pide encontrar la respuesta en el tiempo de la estructura indicada en
% la figura para cada una de sus direcciones principales.
% Se pide determinar tanto la respuesta amortiguada, como la no amortiguada
% Se debe determinar la respuesta para los siguientes amortiguamientos: 5%,
% 10%, y 20%.
% Las condiciones de borde son: u(0)=20 [cm]

clear all
clc

% Calculo de frecuencias naturales (Ver hoja de resolucion)
wn.x=32.36; %[1/s]
wn.y=20.21; %[1/s]

% Calculamos los periodos (Esto es importante para determinar el tiempo de la respuesta que debemos graficar) 
Tn.x=2*pi()/wn.x;
Tn.y=2*pi()/wn.y;

Tn

% Determinamos el tiempo de la respuesta
t=[0:.001:3];

% Determinamos las condiciones de borde del problema
xo=20; %[cm]
vo=0; %[cm/s]

% Calculo de la respuesta en vibracion libre:
    % Resolucion de la ecuacion para vibracion libre sin amortiguamiento
    %   Los datos requeridos son: [xt,vt,at] = SDOF_VL(t,wn,xo,vo)
    %   [t] → intervalo de tiempo
    %   wn → frecuencia natural del sistema
    %   xo → desplazamiento inicial
    %   vo → velocidad inicial
[x_x_na,v_x_na,a_x_na]=SDOF_VL(t,wn.x,xo,vo);
[x_y_na,v_y_na,a_y_na]=SDOF_VL(t,wn.y,xo,vo);

% Calculo de la respuesta en vibracion amortiguada
    % Resolucion de la ecuacion para vibracion libre con amortiguamiento
    %   Los datos requeridos son: [xt,vt,at] = SDOF_VLA(t,wn,xo,vo,xi)
    %   [t] → intervalo de tiempo
    %   wn → frecuencia natural del sistema
    %   xo → desplazamiento inicial
    %   vo → velocidad inicial
    %   xi → ratio de amortiguamiento

xi=[0.05 0.10 0.20];

for i=1:length(xi)
[temp] = SDOF_VLA(t,wn.x,xo,vo,xi(1,i));
x_x_a{i}=temp;
end

% Graficamos la respuesta
figure
plot(t,x_x_na)
xlabel('Tiempo [s]')
ylabel('Desplazamiento x(t)')
title('Respuesta en el tiempo')
hold on

for i=1:length(xi)
plot(t,x_x_a{i})
hold on
end
legend('Sin amortiguamiento','Xi=0.05','Xi=0.10','Xi=0.20')