% sintesi del volano
% visualizzazione dell'andamento della velocità

clear all
clf

% dati geometrici e inerziali
AB=0.4;
BC=1.2;
IA=0.2667;
ID=0.2;
mC=12;

qpunto0=6;      % velocità iniziale
n=3000;          % suddivisione dell'angolo giro
q=0;            % angolo di manovella
q_finale=2*pi;  % angolo finale


% coppie di disturbo e motrice
Cmax=7;        % massimo della coppia di disturbo che ha andamento Cr=Cmax*sin(q/2)[Nm]
Cm=2*Cmax/pi;   % coppia motrice costante durante un periodo di rotazione

% calcolo di A0
% quando la manovella è orizzontale il rapporto di trasmissione è nullo
A0=IA;
delta_q=q_finale/n;       % incrementi dell'angolo di manovella

% Energia cinetica iniziale
E0=1/2*A0*qpunto0^2;

E=E0;           %inizializzazione dell'energia
for i=1:n
    Cr=Cmax*sin(q/2);               % coppia di disturbo
    E=E+(Cm-Cr)*delta_q;            % aggiornamento dell'energia

    phi2=-asin((AB/BC)*sin(q));         % calcolo dell'angolo di biella
    tau=(cos(q)*tan(phi2)-sin(q))*AB;   % rapporto di trasmissione
    A=IA+mC*tau^2;                      % inerizia ridotta

    vet_qpunto(i)=sqrt(2*E/A);          % velocità di manovella
    vet_q(i)=q;                         % memorizzazione dell'angolo di manovella
    q=q+delta_q;                        % incremento dell'angolo di manovella
end

figure(1)
hold on
plot(vet_q,vet_qpunto,'r');             % velocità angolari senza volano
grid
axis equal