% Geislingeruv hltic
format longG

% 1) Zvolen moment setrvacnosti vzhledem k rozmerum -> cim vetsi, tim lepsi
% Moment setrvacnosti disku
I_vnitrni=0.00025509          % [kgm2]
I_vnejsi=0.0145303          % [kgm2]
I_celk=I_vnitrni+I_vnejsi   % [kgm2]

% 2) Zvolena frekvence s pripojenou setrvacnosti hltice, na kterou bude hltic naladen
omega=641   % frekvence s pripojenou setrvacnosti hltice [rad/s]

% 3) Vypoctena tuhost hltice
k_vypoctena=I_vnejsi*omega^2 % tuhost Nm/rad

% Nyni byly v grafu prenosove funkce v souboru zjednoduseny_nahradni_model_smykadla.slx
% v zalozce "Apps/Model linearizer/Bode" urceny dva pracovni body (FP1 a FP2)
% v pruseciku puvodni prenosove funkce a prenosove funkce s hlticem.
% Temito body musi prochazet tlumena i netlumena funkce.
% Tuhost je nastavena tak, aby pracovni body mely stejnou aplitudu.
% Pokud ma FP1 vyssi amplitudu nez FP2, je treba snizit
% tuhost. Pokud ma FP1 nizsi amplitudu nez FP2, je treba tuhost zvysit. Pokud je
% amplituda obou pracovnich bodu prilis vysoka, je mozne ji zvysenim momentu
% setrvacnosti snizit.
k_zvolena=5900               % [Nm/rad]

% 4) Odhad tlumeni podle prenosove funkce tak, aby krivka mezi FP1 a FP2 byla plocha
b_geis_zvolene1=2   % Zvolene tlumeni podle grafu [Nms/rad]

% Konstrukce sestav listovych pruzin
n_sestav=4        % pocet sestav pruzin
E_planzeta=206e9  % Yanguv modul pruznosti pruzinove oceli TiAl6V4 (titanova slitina) [Pa]
b=0.00515         % sirka planzety [m]
a=0.001           % tloustka planzety [m]
l_planzety=0.028          % delka planzety [m]
A_planzety=l_planzety*b   % plocha planzety [m2]
R=0.040                   % prumer pod listem pruziny [m]
k_geis=(((R+l_planzety)^2)*E_planzeta*a*n_sestav*b^3)/(4*l_planzety^3) % skutecna tuhost hltice [Nm/rad]

% Parametry silikonu v hltici
ro_olej=971  % hustota silikonu kg/m3
ny_olej=0.100  % kinematicka viskozita m2/s (1 cSt = 1 mm2/s) - odpovida 108 000 cSt (centistokes)
eta=ny_olej*ro_olej % dynamicka viskozita Pas

% Urceni skutecneho tlumeni
R_drazka=0.044       % polomer stredu drazky od osy otaceni hridele
a_drazka=0.0015      % sirka drazky [m]
b_drazka=0.005      % vyska drazky [m]
l_drazka=0.008    % delka tunelu otevreneho do dutiny [m]
b_geis=(12*A_planzety^2*l_drazka*eta*R_drazka^2)/(a_drazka*b_drazka^3) % skutecne tlumeni hltice [Nms/rad]

% Namahani teplem
ksi=b_geis/(2*sqrt(k_geis*I_vnejsi))
M_vibrace=8.207   % Moment vibraci pohlceny hlticem - odecteno z grafu [Nm]
n_smyk=436         % otacky smykadla [ot/min]
f_buzeni=58.0244   % frekvence budicich kmitu [Hz]
i=f_buzeni/(n_smyk/60) % frekvence budicich kmitu na jednu otacku smykadla [Hz]
P_ztratovy=(pi/60)*(ksi/(1+ksi^2))*((M_vibrace^2)*i*n_smyk)/k_geis % tepelny vykon hltice [W]
