clc
clear all
%% Vstupní data materiálu 42CrMo4
Rm=1283; % Mez pevnosti v tahu [MPa] 
Re=900; % Mez kluzu [MPa]
sigmac_TT=495; % Mez únavy v tahu/tlaku [MPa]
sigmac_OHYB=525; % Mez únavy v tahu/tlaku [MPa]
% Ostatní vstupní data
sigma_VMA=392.7; % Maximální napětí von Mises, zátěžový stav A [MPa] 
sigma_VMB=74.9; % Maximální napětí von Mises, zátěžový stav B [MPa]
sigma_1A=434.5; % První hlavní napětí, zátěžový stav A [MPa] 
sigma_1B=53.2; % První hlavní napětí, zátěžový stav B [MPa] 
sigma_3A=26.5; % Třetí hlavní napětí, zátěžový stav A [MPa] 
sigma_3B=0.9; % Třetí hlavní napětí, zátěžový stav B [MPa] 

sigma_eX=392.7; % Maximální hodnota lokálního napětí [MPa] 
sigma_eX1=164.8; % Napětí ve vzdálenosti xx1 pod povrchem [MPa]
xx1=2.26; % Vzdálenost dvou krajních uzlů elementu [mm]
d_vzorek=7.5; % Průměr zkušebního vzorku [mm]
Do=66; %Vnější průměr ojničního čepu [mm]
eta_sigma=1.189*Do^-0.097; % Součinitel vlivu velikosti [-]
ny_sigma=0.75; % Vliv pravděpodobnosti přežití [-]
%% Tuhost pro uložení klikového hřídele
Fpmax=99349; % Maximální síla od tlaku plynů [N]
cl=0.09; % Ložisková vůle [mm]
ks=(0.5*Fpmax)/(0.8*cl); % Tuhost pro stlačení ve středu hlavního ložiska [Nmm-1]
%% Stanovení bezpečnosti
kappaR=(1/sigma_eX)*((sigma_eX-sigma_eX1)/xx1); % Poměrný gradiet [mm-1]
fG=1+((((sigmac_OHYB/sigmac_TT)-1)/(2/d_vzorek))*kappaR); % Korekční součinitel [mm-1]
beta_alfa=1+(sqrt(kappaR)*10^-(0.35+(Re/810))); % Poměr součinitele vrubu a koncentrace napětí [-]

sigma_emax=sigma_VMA; % Maximální ekvivalentní napětí [MPa]
sigma_emin=sigma_VMB; % Minimální ekvivalentní napětí [MPa]
sigma_ea=(sigma_emax-sigma_emin)/2; % Amplituda ekvivalentního napětí [MPa]
sigma_em=(sigma_emax+sigma_emin)/2; % Střední hodnota ekvivalentního napětí [MPa]

kz=1/((beta_alfa*(sigma_ea/(sigmac_OHYB*eta_sigma*ny_sigma*fG)))+(sigma_em/Rm)); % Bezpečnost prostého zatěžování [-]
disp(['Bezpečnost prostého zatěžování = ',num2str(kz),])
k_kaleno=1.3*kz; % Bezpečnost po zakalení přechodových rádiusů [-]
disp(['Bezpečnost po zakalení rádiusů = ',num2str(k_kaleno),])