clc
%wszystko w m i N
Q=17658;
Lp=1.6*sqrt(Q);
Lk=50;
Lw2=2*(Lp+Lk)/1000;
Dzn1=0.053;
sigma_kryt1=274000000; %% zeby dr2 wyszla w m
dr2=sqrt((4*10*Q+pi*sigma_kryt1*Dzn1^2)/(pi*sigma_kryt1)) %% 61mm

Imin=(pi/64)*(0.064^4-Dzn1^4);
Ar2=(pi/4)*(dr2^2-Dzn1^2);
i_min=sqrt(Imin/Ar2);
lambda2=Lw2/i_min;
sigma_plast=280;
E=2.1*100000;
lambda_graniczna=sqrt(((pi^2)*E)/(0.9*sigma_plast));

b2=(280-252)/lambda_graniczna;
a2=280;
sigma_kryt2=a2-b2*lambda2; %273 MPa

%sigma zastepcza

Ms=0.5*Q*0.07*tan(8.5*pi/180); %Nm
Wo=(pi/16)*((0.064^4-Dzn1^4)/0.064);
tau=Ms/Wo;
sigma_c=Q/Ar2;
sigma_z=sqrt((sigma_c^2)+3*(Ms/Wo)^2) %% Pa, 28 MPa

xw2=sigma_kryt2/(sigma_z/1000000)

%nakretka
P=10;
D1=0.065; %% tr 75
d=0.075; %% tr 75
d3=0.064;
En=100000;
D4=0.076;
Az2=(pi/4)*(d^2-D1^2);
Pdopr2=15000000*(31/70)
n2=Q/(Az2*Pdopr2)
ncalk=n2+1.5; %%4
Hn2=ncalk*P; %%40


%rownosc odksztalcen - to jest lepsze

Dzn2=sqrt(((d3^2)*E+(D4^2)*En)/En) % 120
Dw0=D4+0.002



% to drugie gowno
Dzn2_2=sqrt(4*Q/(40*pi)+78^2) % 89
%