clear all
clc

Db=368;       %Prùmìr bubnu v mm
sb=24;          %Tlouš�ka stìny bubnu [mm]
dl=14;          %Prùmìr lana [mm]
t=16;              %Rozteè závitù bubnu [mm]
zb=60;          %Poèet závitù bubnu [1]
nb=4;           %Poèet nyvinutých vrstev [1]
Fl=40400;   %Síla v lanì [N]

b=t*zb;

delta=sqrt(dl^2-(t/2)^2);

l=0;
for i=1:nb
D(i)=Db+2*delta*(i-1);
l=l+pi*zb*D(i);
end

l=(l-3*pi*Db)/1000 %Pracovní délka lana [m] 

Mk=D(nb)/2*Fl;

Rc=Db/2+delta*(nb-1)+dl/2;
Dmax=Rc*2

D1=Db-dl;

Dstr=(D1+Dmax)/2

D0=D1-2*sb;

Wo=pi*D1^3/32*(1-(D0/D1)^4);
Wk=pi*D1^3/16*(1-(D0/D1)^4);

sigmaO=(Fl/2*(b+80)/2)/Wo; %napìtí v ohybu
tauk=Mk/Wk;                               %napìtí v krutu
sigmatl=Fl/(sb*t);                       %sevøení lanem

sigmared=sqrt(sigmaO^2+sigmatl^2-sigmaO*sigmatl+3*tauk^2)%Redukované napìtí [Nmm-2]