%Função para cálculo dos parâmetros do material viscoelástico
function [G, Ga] = viscoelastico(omega, omegaa)

%Declarando propriedades do material viscoelástico
%***************************
%Material: Borracha Butílica
%***************************
T0 = 273; %[K]
theta1 = 15.1; 
theta2 = 171; 
G0 = 1.53e6; %[N/m2]
Goo = 1.11e8; %[N/m2]
beta = 0.396; 
b1 = 1.34e-2; 

%Informar aqui a temperatura de trabalho do material viscoelástico
T = 25 + 273;

%Cálculo de alpha(T)
alpha = 10^((-theta1*(T-T0))/(theta2+T-T0));

%omega = 0:0.5:2*2*pi*10000;

for iomeg=1:length(omega)

    G(iomeg) = (G0 + Goo*b1*(i*(omega(iomeg))*alpha)^(beta))/(1 + b1*((i*(omega(iomeg))*alpha)^(beta)));
    
end

    Ga = (G0 + Goo*b1*((i*(omegaa))*alpha)^(beta))/(1 + b1*((i*(omegaa)*alpha)^(beta)));

% figure(2)
% subplot(2,1,1)
% loglog((omega*alpha), (real(G)));
% subplot(2,1,2)
% loglog((omega*alpha), (imag(G)./real(G)));