% Sistema com 1 gdl sob desbalanceamento

clc;
clearvars;
close all;

% dados do sistema
m=120; % massa, kg
c=1000; % coeficiente de amortecimento, kg/s
k=8e5; % constante de rigidez, N/m

% parâmetros modais
omn=sqrt(k/m); % frequência natural, rad/s
zeta=c/(2*m*omn); % razão de amortecimento

% rotação de operação
rop=3000; % rpm
omo=rop*2*pi/60; % rad/s

% condições de desbalanceamento
F=374; % amplitude da força, N
me=0.01*m; % massa desbalanceada, kg

% item (a)

% excentricidade
ec=F/(me*omo^2); % m

% visualização de parâmetros modais,
% rotação de operação e excentricidade
disp(' ')
disp(['A frequência natural circular é ',...
      'igual a ',num2str(omn),' rad/s.'])
disp(' ')
disp(['A razão de amortecimento é igual a ',...
       num2str(zeta),'.'])
disp(' ')
disp(['A rotação de operação é igual a ',...
       num2str(omo),' rad/s.'])
disp(' ')
disp(['A excentricidade é igual a ',...
       num2str(ec),' m.'])
disp(' ')

% % item (b)
% 
% % amplitude de vibração na rotação de operação
% eop=omo/omn;
% Hop=1/sqrt((1-eop^2)^2+(2*zeta*eop)^2);
% Yop=(me*ec*eop^2*Hop)/m;
% 
% % amplitude de vibração na rotação crítica
% ecr=1;
% Hcr=1/sqrt((1-ecr^2)^2+(2*zeta*ecr)^2);
% Ycr=(me*ec*ecr^2*Hcr)/m;
% 
% % visualização de amplitudes de vibração,
% % nas rotações de operação e crítica
% disp(' ')
% disp(['A amplitude de vibração na rotação ',...
%       'de operação é ',num2str(Yop),' m.'])
% disp(' ')
% disp(['A amplitude de vibração na rotação ',...
%       'crítica é ',num2str(Ycr),' m.'])
% disp(' ')