clear all; clc;
%Parametry  
m=2325;     %pohotovostní hmotnost vozu v kg
S_n=265;    %nominální šířka pneumatiky v mm
AR=35;      %poměrná výška profilu pneumatiky
D_r=21;     %průměr ráfku kola v palcích
g=9.81;     %gravitační zrychlení v m/s^-2

D=S_n*AR/50+D_r*25.4;
R=D/2;
G=m*g/4; 


%Výpočet závislosti stahovací síly na poloměru zvedacích válců
r1=20:1:35;  %průměr zvedacích válců v mm
h1=60;       %výška uchycení zvedacích válců v mm
alfa01=asin((R-h1)./(R+r1));   
F_max1=(G/2)./(tan(alfa01));
L_max1=2.*(R+r1).*cos(alfa01);

figure('Name','Závislost velikosti maximální stahovací síly na poloměru zvedacích válců','NumberTitle','off');
plot(r1,F_max1, 'Linewidth', 1.5,'Color',[0.5255, 0.2549, 0.6784]);
title({'Závislost velikosti maximální stahovací síly' 'na poloměru zvedacích válců'});
xlabel('Poloměr zvedacích válců [mm]'); ylabel('Stahovací síla [N]');
set(gca,'FontName','Technika Book','FontSize', 12,'XGrid','on','YGrid','on');
print -djpeg Zavislost1

figure('Name','Závislost maximální rozteče zvedacích válců na jejich poloměru','NumberTitle','off');
plot(r1,L_max1,'Linewidth', 1.5,'Color',[0.3294, 0.6275, 0.0471]);
title({'Závislost maximální rozteče zvedacích válců' 'na jejich poloměru'});
xlabel('Poloměr zvedacích válců [mm]'); ylabel('Maximální rozteč zvedacích válců [mm]');
set(gca,'FontName','Technika Book','FontSize', 12,'XGrid','on','YGrid','on');
print -djpeg Zavislost2


%Výpočet závislosti stahovací síly na výšce uchycení zvedacích válců
r2=25;       %průměr zvedacích válců v mm
h2=30:1:60;  %výška uchycení zvedacích válců v mm
alfa02=asin((R-h2)./(R+r2));   
F_max2=(G/2)./(tan(alfa02));
L_max2=2.*(R+r2).*cos(alfa02);

figure('Name','Závislost velikosti maximální stahovací síly na výšce uchycení zvedacích válců','NumberTitle','off');
plot(h2,F_max2, 'Linewidth', 1.5,'Color',[0.5255, 0.2549, 0.6784]);
title({'Závislost velikosti maximální stahovací síly' 'na výšce uchycení zvedacích válců'});
xlabel('Výška uchycení zvedacích válců [mm]'); ylabel('Stahovací síla [N]');
set(gca,'FontName','Technika Book','FontSize', 12,'XGrid','on','YGrid','on');
print -djpeg Zavislost3

figure('Name','Závislost maximální rozteče zvedacích válců na výšce jejich uchycení','NumberTitle','off');
plot(h2,L_max2,'Linewidth', 1.5,'Color',[0.3294, 0.6275, 0.0471]);
title({'Závislost maximální rozteče zvedacích válců'  'na výšce jejich uchycení'});
xlabel('Výška uchycení zvedacích válců [mm]'); ylabel('Maximální rozteč zvedacích válců [mm]');
set(gca,'FontName','Technika Book','FontSize', 12,'XGrid','on','YGrid','on');
print -djpeg Zavislost4


%Výpočet průběhu stahovací síly
r=25;        %průměr zvedacích válců v mm
h=40;        %výška uchycení zvedacích válců v mm
alfa0=asin((R-h)/(R+r));  

alfa=linspace(alfa0,pi/2,50);
L=2*(R+r).*cos(alfa);
F=(G/2)./(tan(alfa));
y=h+r.*sin(alfa)-R.*(1-sin(alfa));

figure('Name','Průběhy stahovací síly a zdvihu kola při uvažování absolutně tuhých těles','NumberTitle','off');
yyaxis left
plot(L,F, 'Linewidth', 1.5);
title({'Průběhy stahovací síly a zdvihu kola' 'při uvažování absolutně tuhých těles'});
xlabel('Rozteč zvedacích válců [mm]'); ylabel('Stahovací síla [N]');
yyaxis right
plot(L,y, 'Linewidth', 1.5);
ylabel('Zdvih kola [mm]');
set(gca,'FontName','Technika Book','FontSize', 12,'XGrid','on','YGrid','on');
print -djpeg PrubehAbsTuha

Fmax=max(F)
Lmax=max(L)
