%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% RUN SDOF OSCILLATOR
%
% created by Michalis Fragiadakis, Dec 2013
% mfrag@mail.ntua.gr
%
% please report any bugs
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
includeFile

% RUN A NONLINEAR SDOF
Tn = 1; %sec
xi = 0.05;
dt = acc.dtacc;
hardening = 0.01; % post yield over elastic modulus
fy = 0.01;
nsteps = acc.nsteps;
niter  = 10;
sfactor = 1; % record scaling factor
iplot = 0; % draw hysteretic plot

[umax,vmax,amaxRel,amaxTot,fsmax,u,fs]=...
    sdfNL(Tn,xi,dt,hardening,fy,acc.rec,nsteps,niter,sfactor,iplot);


% plot hysteretisis
figure()
hold on; grid on; box on;
plot(u,fs,'k.-')

xlabel('displacement','FontSize',18)
ylabel('force','FontSize',18)


% plot displacement over time
figure()
hold on; grid on; box on;
plot((1:(nsteps-1))*dt,u,'k-')

xlabel('time (sec)','FontSize',18)
ylabel('displacement','FontSize',18)

% plot force over time
figure()
hold on; grid on; box on;
plot((1:(nsteps-1))*dt,fs,'k-')

xlabel('time (sec)','FontSize',18)
ylabel('force','FontSize',18)