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

% dati geometrici e inerziali
AB=0.4;
BC=1.2;
IA=0.2667;
ID=0.2;     % momento d'inerzia del volano 
mC=12;

qpunto0=6;      % velocità iniziale
n=300;          % suddivisione dell'angolo giro
q=0;            % angolo di manovella

% calcolo di A0
% quando la manovella è orizzontale il rapporto di trasmissione è nullo
A0=IA;
delta_q=2*pi/n;       % incrementi dell'angolo di manovella
for i=1:n
    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_q(i)=q;                         % salavataggio deigli angoli di manovella nel vettore
    vet_qpunto(i)=sqrt(A0/A)*qpunto0;   % calcolo velocità angolare di manovella e memorizzazione
    
    % introduzione del volano
    A=IA+ID+mC*tau^2;                      % inerizia ridotta
    vet_qpunto_volano(i)=sqrt((A0+ID)/A)*qpunto0;   % calcolo velocità angolare di manovella e memorizzazione

    
    q=q+delta_q;                        % incremento dell'angolo di manovella
end

figure(1)
hold on
plot(vet_q,vet_qpunto,'r');             % velocità angolari senza volano
plot(vet_q,vet_qpunto_volano,'g');             % velocità angolari senza volano

grid
axis equal