No BSD License  

Highlights from
Calculating blending moment of the shaft

image thumbnail
from Calculating blending moment of the shaft by Riccardo
It calculate blending moment of the shaft when you insert a generic number of forces.

SottoGrafico (reazioni,forze)
function SottoGrafico (reazioni,forze)
i=0;
lunghezza=evalin('base','lunghezza');
position=[0 lunghezza];

% Disegna le reazioni vincolari sul grafico
for i=1:2
    if reazioni(1,i)>0
        react=line([position(i) position(i)], [0 -reazioni(1,i)],'marker','.','Color','r');        
        f_pallucca(position(i),0,reazioni(1,i)/10,'r');
    else
        react=line([position(i) position(i)], [0 reazioni(1,i)],'marker','.','Color','r');
        f_pallucca(position(i),reazioni(1,i),reazioni(1,i)/10,'r');        
    end
end

% Disegna le forze agenti sull'albero
i=0;
nFor=evalin('base','numeroforze');
for i=1:nFor
    if forze(i,2)>0
        forza=line([forze(i,3) forze(i,3)], [0 forze(i,2)],'marker','.','Color','b');
        f_pallucca(forze(i,3),forze(i,2),forze(i,2)/10,'b');
    else
        forza=line([forze(i,3) forze(i,3)], [0 -forze(i,2)],'marker','.','Color','b');
        f_pallucca(forze(i,3),0,forze(i,2)/10,'b');
    end
    text(forze(i,3),abs(forze(i,2)/2),[num2str(forze(i,2)),'KN']);
    a(i,:)=[forze(i,2)];
end
reazmax=max(abs(reazioni));
forzamassima=max(abs(a));
ylabel('Forze in KN'); xlabel('Lunghezza dell''albero in mm');
axis ([-position(2)/10 1.1*position(2) -reazmax forzamassima]);

Contact us at files@mathworks.com