% Znaménková konvence dle DP
%% Počáteční příkazy
clear all                                                   % smazání všech proměnných
clf                                                         % vyčistit figures
clc                                                         % vyčistit cmd line
set(0,'DefaultFigureWindowStyle','docked')                  % kouzlo na automatické dokování oken
format compact                                              % zápis
format shortEng                                             % zápis
%% Zadané hodnoty
p0=101325;                                                  % atmosférický tlak (v klikové skříni)
vrtani=76.5;                                                % [mm] vrtání
zdvih=86.9;                                                 % [mm] zdvih
D=vrtani/1000;                                              % [m]
H=zdvih/1000;                                               % [m]
l=0.138;                                                    % [m] délka ojnice
mpc=0.394;                                                  % hmotnost posuvných částí (pístní skupina + část ojnice)
m2=0.334;                                                   % hmotnost rotačních částí ojnice
PV=[0 540 180 360];                                         % pořadí zážehů 1-3-4-2
koeficient_treni_prumerny=0.05;                             % [-]
a_vzdalenost_os_valcu=0.082;                                % vzdálenost os válců
%% Vypočtené hodnoty
r=H/2;                                                      % [m] poloměr kliky
lambda=r/l;                                                 % klikový poměr
Spistu=(pi*D^2)/4;                                          % [m2] plocha pístu
%% Načtení reálných hodnot tlaků centrického KM
hodnoty_tlaku_centricky=('SKODA-indik-diagramy.txt');       % načtení datového souboru
HTC=importdata(hodnoty_tlaku_centricky, '\t',2);            % \t oddělovačem je tabulátor; počet řádků, které se přeskakují
alfa_deg=HTC.data(:,1);                                     % načtení tlaků pro jednotlivé otáčky
alfa=(alfa_deg/180)*pi;                                     % převod natočení klikového hřídele na radiány pro další výpočty
beta=asin(lambda*sin(alfa));                                % úhel odklonu ojnice
%% Výpočet základních kinematických veličin a sil v centrickém klikovém mechanismu
for i=2:1:18                                                % beru 17 sloupců hodnot tlaků pro celý rozsah otáček
tlak_bar(:,i)=HTC.data(:,i);                                % načtení hodnot tlaků do proměnné tlak_bar
    for k=1:length(alfa)                                    % spočítání kinematiky klikového mechanismu a sil působících v klikovém mechanismu
                                                            % v závislosti na natočení klikového hřídele a daných otáčkách 
        n(i)=(1000+(i*250))/60;                                 % přiřazení otáček jednotlivým sloupcům i
        omega(i)=2*pi*n(i);                                     % vypočtení úhlové rychlosti [rad/s] pro dané otáčky  
    
        s(k)=r*((1-cos((alfa(k)))+((lambda/4)*(1-cos(2*alfa(k))))));    % dráha pístu v závislosti na natočení klikového hřídele
        v(k,i)=r*omega(i)*(sin(alfa(k))+(lambda/2)*sin(2*alfa(k)));     % rychlost pístu v závislosti na natočení klikového hřídele
        a(k,i)=r*omega(i)^2*(cos(alfa(k))+lambda*cos(2*alfa(k)));       % zrychlení pístu v závislosti na natočení klikového hřídele
   
        Fp(k,i)=(tlak_bar(k,i)*1e5-p0)*Spistu;                  % síla od tlaků plynů působící na píst
        Fs(k,i)=-a(k,i)*mpc;                                    % setrvačná síla působící na píst
        Fvp(k,i)=Fp(k,i)+Fs(k,i);                               % výsledná síla působící na píst
   
        Fn(k,i)=Fvp(k,i).*tan(beta(k));                          % boční síla působící na píst -> síla mezi pláštěm pístu a válcem
        Fo(k,i)=Fvp(k,i)./(cos(beta(k)));                        % výsledná síla v ojnici
        Ft(k,i)=Fo(k,i).*(sin(alfa(k)+beta(k)));                 % tečná složka síly na ojničním čepu
        Mk(k,i)=Ft(k,i)*r;                                       % točivý moment na ojničním čepu KH
    
                Mkv1=circshift(Mk,[PV(1),0]);                   % posunutí průběhu točivého momentu 1. válce
                Mkv2=circshift(Mk,[PV(2),0]);                   % posunutí průběhu točivého momentu 2. válce
                Mkv3=circshift(Mk,[PV(3),0]);                   % posunutí průběhu točivého momentu 3. válce
                Mkv4=circshift(Mk,[PV(4),0]);                   % posunutí průběhu točivého momentu 4. válce
                Mkc=Mkv1+Mkv2+Mkv3+Mkv4;                        % celkový průběh točivého momentu
        
        Fro(k,i)=-m2*r*omega(i)^2;                              % odstředivá síla na ojniční čep
        Fr(k,i)=Fo(k,i).*(cos(alfa(k)+beta(k)));                % radiální síla na OČ
        Frc(k,i)=Fr(k,i)+Fro(k,i);                              % výsledná radiální síla na OČ
        Fvoc(k,i)=sqrt(Frc(k,i)^2+Ft(k,i)^2);                   % výsledná síla na OČ
    end
    
        
end

%%
figure(40)
hold on
plot(alfa_deg,Ft(:,11),'g-','LineWidth',2)
plot(alfa_deg,Fr(:,11),'b-','LineWidth',2)
plot(alfa_deg,Frc(:,11),'r-','LineWidth',2)
plot(alfa_deg,Fvoc(:,11),'k--','LineWidth',2)
title('Síly působící na ojniční čep')
                xlabel('[°]')                               % popis osy x
                ylabel('[N]')                               % popis osy y
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
                grid on
                legend('Tečná síla přenášená ojnicí','Radiální síla přenášená ojnici','Výsledná radiální síla na ojničním čepu','Výsledná zatěžující ojniční čep')
    hold off 

%% Vykreslování grafů pro CENTRICKÝ MECHANISMUS

figure(1)                                                   % obrázek 1
plot(alfa_deg,tlak_bar(:,2:end),'LineWidth',2)              % vykreslení průběhu tlaků pro jednotlivé otáčky 
title('Průběhy tlaku centrického klikového mechanismu pro otáčky 1500 až 5500')
                xlabel('[°]')                               % popis osy x
                ylabel('[bar]')                             % popis osy y
                xlim([0,720])                               % omezení rozsahu osy x
                set(gca,'XTick',0:180:720)                  % kroky osy x po 180°
                legend on
                grid on   

figure(2)                                                   % obrázek 2
plot(alfa_deg,s,'LineWidth',2)                              % vykreslení
title('Dráha pístu')
                xlabel('[°]')                               % popis osy x
                ylabel('[m.s^-1]')                          % popis osy y
                xlim([0,720])                               % omezení rozsahu osy x
                set(gca,'XTick',0:180:720)                  % kroky osy x po 180°
                grid on                                     % zapnutí mřížky               
                
figure(3)                                                   % obrázek 2
plot(alfa_deg,v(:,2:end),'LineWidth',2)                     % vykreslení
title('Rychlosti pístu pro otáčky 1500 až 5500')
                xlabel('[°]')                               % popis osy x
                ylabel('[m.s^-1]')                          % popis osy y
                xlim([0,720])                               % omezení rozsahu osy x
                set(gca,'XTick',0:180:720)                  % kroky osy x po 180°
                grid on                                     % zapnutí mřížky
               
figure(4)
plot(alfa_deg,a(:,2:end),'LineWidth',2)                         % vykreslení
title('Zrychlení pístu pro otáčky 1500 až 5500')
                xlabel('[°]')                               % popis osy x
                ylabel('[m.s^-2]')                          % popis osy y
                xlim([0,720])                               % omezení rozsahu osy x
                set(gca,'XTick',0:180:720)                  % kroky osy x po 180°
                grid on                                     % zapnutí mřížky
                               
figure(5)
    hold on
plot(alfa_deg,Fp(:,2:end),'g-','LineWidth',1)
plot(alfa_deg,Fs(:,2:end),'b-','LineWidth',1)
plot(alfa_deg,Fvp(:,2:end),'r-','LineWidth',1)
title('Síly působící na píst - od tlaků plynů(g), setrvačná(b), výsledná(r)')
                xlabel('[°]')                               % popis osy x
                ylabel('[N]')                               % popis osy y
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
                grid on
    hold off       

figure(6)
plot(alfa_deg,Fn(:,2:end),'LineWidth',1)
title('Boční síla působící na píst')
                xlabel('[°]')                               % popis osy x
                ylabel('[N]')                               % popis osy y
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
                grid on

figure(7)
plot(alfa_deg,Fo(:,2:end),'LineWidth',1)
title('Výsledná síla působící v ojnici')
                xlabel('[°]')                               % popis osy x
                ylabel('[N]')                               % popis osy y
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
                grid on

figure(8)
plot(alfa_deg,Ft(:,2:end),'LineWidth',1)
title('Tečná síla působící na ojniční čep')
                xlabel('[°]')                               % popis osy x
                ylabel('[N]')                               % popis osy y
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
                grid on

figure(9)
plot(alfa_deg,Mkc(:,2:end),'LineWidth',2)
                title('Průběh točivého momentu')
                xlabel('[°]')                               % popis osy x
                ylabel('[Nm]')                              % popis osy y
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
                grid on
                        
figure(10)
plot(alfa_deg,Mkc(:,11),'LineWidth',2)
                title('Průběh točivého momentu v otáčkách 3750min-1')
                xlabel('[°]')                               % popis osy x
                ylabel('[Nm]')                              % popis osy y
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
                grid on
                
%% Výkon a točivý moment CENTRICKÉHO MECHANISMU
Mt_stredni=mean(Mkc);                                       % průměrný (střední) točivý moment motoru
for t=1:18                                                  % pro všechny sloupce, tzn. otáčky od 1500 do 5500min-1, provedu výpočet výkonu z točivého momentu
P(:,t)=(mean(Mkc(:,t))*((2*pi*(1000+(t*250))/60)))/1000;
end

otacky=[0 1500 1750 2000 2250 2500 2750 3000 3250 3500 3750 4000 4250 4500 4750 5000 5250 5500];    %pomocná matice pro zobrazení do grafu

A=[otacky; Mt_stredni];                                     % pomocná matice pro točivý moment
B=[otacky; P];                                              % pomocná matice pro výkon

figure(11)
hold on
plot(otacky(2:end),A(2,2:end),'b-','LineWidth',2)
plot(otacky(2:end),B(2,2:end),'r-','LineWidth',2)
                title('Průběh točivého momentu a výkonu motoru v závislosti na otáčkách')
                xlabel('[min-1]')                           % popis osy x
                ylabel('[Nm]')                              % popis osy y
                xlim([1000,6000])                           % graf bude ukončen na x=720°
                ylim([10,180])                              % graf bude ukončen na x=720°
                legend('Průbeh točivého momentu','Průběh výkonu')
                grid on
hold off
%% CENTRICKÝ KONEC

%% KINEMATIKA EXCENTRICKÉHO MECHANISMU

%% DRÁHA - členy binomické věty a jejich vliv

for k=1:length(alfa) 
s_br(k)=r*(1-cos(alfa(k)))+l*(1-sqrt(1-((lambda^2)*(sin(alfa(k))*sin(alfa(k))))));

prvni_clen(k)=1;
druhy_clen(k)=0.5*(lambda^2)*(sin(alfa(k))^2);
treti_clen(k)=0.125*(lambda^4)*(sin(alfa(k))^4);

prvni_clen_cely(k)=r*(1-cos(alfa(k)))+l*(1-prvni_clen(k));
druhy_clen_cely(k)=r*(1-cos(alfa(k)))+l*(1-(prvni_clen(k)-druhy_clen(k)));
treti_clen_cely(k)=r*(1-cos(alfa(k)))+l*(1-(prvni_clen(k)-druhy_clen(k)-treti_clen(k)));

end 

prvni_clen_rozdil=(mean(prvni_clen_cely)/mean(prvni_clen_cely));
druhy_clen_rozdil=abs(mean(prvni_clen_cely)-mean(druhy_clen_cely))/mean(prvni_clen_cely);
treti_clen_rozdil=abs(mean(druhy_clen_cely)-mean(treti_clen_cely))/mean(prvni_clen_cely);

sloupec=[prvni_clen_rozdil,druhy_clen_rozdil,treti_clen_rozdil];

figure(21)
hold on
bar([1], sloupec(1),'FaceColor','b');
bar([2:3], sloupec(2:3),'FaceColor','g');
hold off
                title('Velikost členů binomické věty pro daný klikový poměr')
                xlabel('Pořadí členu')                               % popis osy x
                ylabel('Poměrná velikost')                  % popis osy y

figure(22)
hold on
plot(alfa_deg,s_br,'r--','LineWidth',4)
plot(alfa_deg,s,'b','LineWidth',2)
hold off
                title('Dráha dle rozkladu - Binomická rovnice')
                xlabel('[°]')                               % popis osy x
                ylabel('[Nm]')                              % popis osy y
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
                grid on
             
%% Výpočet skutečné dráhy
r_excentricky=r;
l_exentricky=l;
lambda_pro_excentricky=lambda;

for o=1:1:15
    e(o)=(o*2)/1000;
    lambda_e(o)=e(o)/l_exentricky;
    
    
    for k=1:length(alfa)

s_e(k,o)=sqrt(((l_exentricky+r_excentricky)^2)-e(o)^2)-r_excentricky*((1/lambda_pro_excentricky)+cos(alfa(k))+lambda_e(o)*sin(alfa(k))-(lambda/4)*(1-cos(2*alfa(k))));
s_e_presny(k,o)=sqrt(((l_exentricky+r_excentricky)^2)-e(o)^2)-r_excentricky*(+cos(alfa(k))+(1/lambda_pro_excentricky)*(sqrt((1-(lambda_pro_excentricky*sin(alfa(k))-lambda_e(o))^2))));

EXCENTRMIN(o)=min(s_e(:,o));
korekce_zdvih(o)=abs(EXCENTRMIN(o));

s_e_po_korekci(k,o)=s_e(k,o)+korekce_zdvih(o); 

    end
end

EXCENTRMAX=max(s_e);
korekce_EXCENTRMAX=max(s_e_po_korekci);
CENTRMIN=min(s);
CENTRMAX=max(s);
presnyMAX=max(s_e_presny);
presnyMIN=min(s_e_presny);

rozdil_v_mm=(s_e_po_korekci-s_e_presny)*1000;
figure(23)
    hold on
        plot(alfa_deg,s*1000,'k-','LineWidth',2)
        plot(alfa_deg,s_e(:,15)*1000,'r','LineWidth',2)
        plot(alfa_deg,s_e_po_korekci(:,15)*1000,'b-.','LineWidth',3)
        plot(alfa_deg,s_e_presny(:,15)*1000,'g:','LineWidth',3)
    hold off
                xlabel('Natočení klikové hřídele [°]')
                ylabel('Dráha pístu [mm]')
                xlim([0,720])
                set(gca,'XTick',0:5:720)
                legend('Dráha pístu centrického mechanismu','Dráha pístu excentrického mechanismu - 2 členy binomické věty','Dráha pístu excentrického mechanismu - 2 členy binomické věty - včetně korekce','Dráha pístu excentrického mechanismu - plný tvar')
                grid on      
                title('Srovnání dráhy pístu dle způsobu výpočtu')
    
figure(24)
hold on
plot(alfa_deg,s*1000,'k','LineWidth',2)
plot(alfa_deg,s_e_po_korekci(:,2)*1000,'m','LineWidth',1.5)
plot(alfa_deg,s_e_po_korekci(:,4)*1000,'b','LineWidth',1.5)
plot(alfa_deg,s_e_po_korekci(:,7)*1000,'r','LineWidth',1.5)
plot(alfa_deg,s_e_po_korekci(:,10)*1000,'g','LineWidth',1.5)
plot(alfa_deg,s_e_po_korekci(:,15)*1000,'c','LineWidth',1.5)
hold off
                title('Průběh dráhy pístu')
                xlabel('Natočení klikového hřídele [°] ')
                ylabel('Dráha pístu [mm]')
                xlim([0,360])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:90:360)                  % měřítko grafu jde po 180° skocích
                legend('e=0mm','e=4mm','e=8mm','e=14mm','e=20mm','e=30mm')
                grid on

%% Výpočet zdvihu a posunutí úvratí
                
for o=1:1:15
e(o)=(o*2)/1000;
e_mm(o)=e(o)*1000;
ZHU(o)=sqrt(((l_exentricky+r_excentricky)^2)-e(o)^2);
ZDU(o)=sqrt(((l_exentricky-r_excentricky)^2)-e(o)^2);
Zdvih_e=ZHU(o)-ZDU(o)
alfa_HU=asin(e/(l+r));
alfa_DU=asin(e/(l-r))+pi;
rozdil_alfa=(alfa_DU(o)-alfa_HU(o))*(180/pi);
end

alfa_HU_deg=alfa_HU*(180/pi);
alfa_DU_deg=alfa_DU*(180/pi)-180;
Zdvih_e_graf=ZHU-ZDU
rozdil_alfa_graf=(alfa_DU-alfa_HU)*(180/pi)

figure(25)
hold on
[ax,p1,p2] = plotyy(e_mm,Zdvih_e_graf*1000,e_mm,[rozdil_alfa_graf*(180/pi)],'plot','plot');
set(p1,'LineWidth',2)
set(p2,'LineWidth',2)
hold off
                title('Závislost zdvihu a rozdílu v natočení KH mezi úvratěmi na hodnotě excentricity')
                xlabel('Excentricita [mm]')
                ylabel(ax(1),'Zdvih [mm]')
                ylabel(ax(2),'Rozdíl mezi horní a dolní úvratí oproti CKM [°]')
                grid on
                xlim(ax(1),[2,30])
                xlim(ax(2),[2,30])
                
% Výpočet změny zdvihového objemu celého motoru
format long
zdvih_puvodni=86.9/1000;
zdvih_pokus=87.401/1000;
vrtani_cc=76.5/1000;

Vpuvodni=4*zdvih_puvodni*((pi/4)*(vrtani_cc^2))
Vnovy=4*zdvih_pokus*((pi/4)*(vrtani_cc^2))
e14=14;
r_mm=r*1000;
for nn=1:1:500
    l_v_o(nn)=nn;
    ZHU_v_o(nn)=sqrt(((l_v_o(nn)+r_mm)^2)-(e14^2));
    ZDU_v_o(nn)=sqrt(((l_v_o(nn)-r_mm)^2)-(e14^2));
    Z_v_o(nn)=ZHU_v_o(nn)- ZDU_v_o(nn);
    
end
Z_v_o(138:300)
plot(l_v_o(138:300),Z_v_o(138:300))
%% Výpočet rychlosti pístu excentrického klikového mechanismu pro 5500RPM

for o=1:1:15
        e(o)=(o*2)/1000;
        lambda_e(o)=e(o)/l_exentricky;
        n5500=5500;
        omega5500=2*pi*(n5500/60);
    for k=1:length(alfa)   
       v_e_5500(k,o)=r_excentricky*omega5500*(sin(alfa(k))-lambda_e(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
    end
end

figure(26)
hold on
plot(alfa_deg,v(:,18),'k--','LineWidth',1.5)
plot(alfa_deg,v_e_5500(:,2),'m:','LineWidth',1.5)
plot(alfa_deg,v_e_5500(:,4),'y:','LineWidth',1.5)
plot(alfa_deg,v_e_5500(:,6),'b:','LineWidth',1.5)
plot(alfa_deg,v_e_5500(:,8),'g','LineWidth',1.5)
plot(alfa_deg,v_e_5500(:,10),'r','LineWidth',1.5)
plot(alfa_deg,v_e_5500(:,15),'m','LineWidth',1.5)
hold off
                title('Průběh rychlosti pístu pístu při 5500min-1')
                xlabel('natočení klikového hřídele [°]')
                ylabel('rychlost pístu pístu [ms-1]')
                xlim([0,720])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:90:720)                  % měřítko grafu jde po 180° skocích
                legend('e=0mm','e=4mm','e=8mm','e=12mm','e=16mm','e=20mm','e=30mm')
                grid on
                
% výpočet středních rychlostí - centrický KM
str_rychlost_kladna = v(:,18);
Kladna_vetev = str_rychlost_kladna > 0;
str_rychlost_kladna(~Kladna_vetev) = 0;
vysledek_kladna = sum(str_rychlost_kladna) ./ sum(Kladna_vetev)

str_rychlost_zaporna = v(:,18);
Zaporna_vetev = str_rychlost_zaporna < 0;
str_rychlost_zaporna(~Zaporna_vetev) = 0;
vysledek_zaporna = sum(str_rychlost_zaporna) ./ sum(Zaporna_vetev)

% výpočet středních rychlostí - EXCENTRICKÝ KM
%e=20mm
str_rychlost_kladna_excentr = v_e_5500(:,10);
Kladna_vetev_excentr = str_rychlost_kladna_excentr > 0;
str_rychlost_kladna_excentr(~Kladna_vetev_excentr) = 0;
vysledek_str_v_kladna_excentr = sum(str_rychlost_kladna_excentr) ./ sum(Kladna_vetev_excentr)

str_rychlost_zaporna_excentr = v_e_5500(:,10);
Zaporna_vetev_excentr = str_rychlost_zaporna_excentr < 0;
str_rychlost_zaporna_excentr(~Zaporna_vetev_excentr) = 0;
vysledek_str_v_zaporna_excentr = sum(str_rychlost_zaporna_excentr) ./ sum(Zaporna_vetev_excentr)

%e=10mm
str_rychlost_kladna_excentr10 = v_e_5500(:,5);
Kladna_vetev_excentr10 = str_rychlost_kladna_excentr10 > 0;
str_rychlost_kladna_excentr10(~Kladna_vetev_excentr10) = 0;
vysledek_str_v_kladna_excentr10 = sum(str_rychlost_kladna_excentr10) ./ sum(Kladna_vetev_excentr10)

str_rychlost_zaporna_excentr10 = v_e_5500(:,5);
Zaporna_vetev_excentr10 = str_rychlost_zaporna_excentr10 < 0;
str_rychlost_zaporna_excentr10(~Zaporna_vetev_excentr) = 0;
vysledek_str_v_zaporna_excentr10 = sum(str_rychlost_zaporna_excentr10) ./ sum(Zaporna_vetev_excentr10)

%porovnání
stredni_rychlost_celkova_centricky=mean(abs(v(:,18)))   
stredni_rychlost_celkova_excentricky=mean(abs(v_e_5500(:,5)))


%% Načtení dat 5500RPM - pro graf vlivu excentricity na průběh spalovacích tlaků
nacteni_tlaky_excentricky_komplet_5500rpm=xlsread('tlaky_excentrickyKM_usporadane_5500rpm.xlsx');

tlaky_excentr_5500rpm_2az30=nacteni_tlaky_excentricky_komplet_5500rpm;
alfa_deg_exc=tlaky_excentr_5500rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;
ulozeni_tlaku_excentr_2az30mm_5500=tlaky_excentr_5500rpm_2az30(:,2:end);

figure(28)  
hold on
plot(alfa_deg_exc,tlak_bar(:,18),'k--','LineWidth',2)
% plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,5),'g','LineWidth',2)
% plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,10),'b','LineWidth',2)
% plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,15),'r','LineWidth',2)
% legend('e=0mm','e=10mm','e=20mm','e=30mm')

plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,2),'g','LineWidth',2)
plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,4),'b','LineWidth',2)
plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,6),'r','LineWidth',2)
plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,8),'m','LineWidth',2)
plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,10),'y','LineWidth',2)
plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,12),'b--','LineWidth',2)
plot(alfa_deg,ulozeni_tlaku_excentr_2az30mm_5500(:,15),'r--','LineWidth',2)
legend('e=0mm','e=4mm','e=8mm','e=12mm','e=16mm','e=20mm','e=24mm','e=30mm')

hold off
title('Průběhy tlaku excentrického klikového mechanismu')
                xlabel('Natočení klikového hřídele [°]')                               % popis osy x
                ylabel('Tlak ve spalovacím prostoru [bar]')                             % popis osy y
                xlim([0,720])                               % omezení rozsahu osy x
                set(gca,'XTick',0:180:720)                  % kroky osy x po 180°
                grid on 

%% Zrychlení pístu excentrického klikového mechanismu pro 5500 RPM + výpočet střední velikosti zrychlení do DP

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n5500=5500;
        omega5500=2*pi*(n5500/60);

        
    for k=1:length(alfa)   
        
        a_e_5500(k,o)=r_excentricky*(omega5500^2)*(cos(alfa(k))+lambda_e(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
        beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));
    end
end

figure(27)
hold on
plot(alfa_deg,a(:,18),'k:','LineWidth',1.5)
plot(alfa_deg,a_e_5500(:,2),'m','LineWidth',1.5)
plot(alfa_deg,a_e_5500(:,4),'y','LineWidth',1.5)
plot(alfa_deg,a_e_5500(:,6),'g','LineWidth',1.5)
plot(alfa_deg,a_e_5500(:,8),'k','LineWidth',1.5)
plot(alfa_deg,a_e_5500(:,10),'b','LineWidth',1.5)
plot(alfa_deg,a_e_5500(:,15),'r','LineWidth',1.5)
hold off
                title('Průběh zrychlení pístu pístu při 5500min-1')
                xlabel('natočení klikového hřídele [°]')
                ylabel('zrychlení pístu [ms-2]')
                xlim([0,360])                               % graf bude ukončen na x=720°
                set(gca,'XTick',0:90:360)                  % měřítko grafu jde po 180° skocích
                legend('e=0mm','e=4mm','e=8mm','e=12mm','e=16mm','e=20mm','e=30mm')
                grid on
                
str_a_HU_DU_centr=mean(abs(a(1:720,18)))
str_a_DU_HU_centr=mean(abs(a(361:720,18)))

str_a_HU_DU_excentr=mean(abs(a_e_5500(4:188,7)))
str_a_DU_HU_excentr=mean(abs(a_e_5500(189:364,7)))

abscelk_stredni_a_excentr=mean(abs(a_e_5500(:,7)))
abscelk_stredni_a_centr=mean(abs(a(:,18)));           
                
celk_stredni_a_excentr=mean(a_e_5500(:,7))
celk_stredni_a_centr=mean(a(:,18))


%% Ověření výpočtu bety
for k=1:length(alfa_exc)
lambda_e30mm=0.03/l_exentricky;    
beta_exc_30(k)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e30mm);
end               

figure(29)  
hold on
plot(alfa_deg_exc,beta_exc_rozsah(:,15),'k--','LineWidth',2)
plot(alfa_deg_exc,beta_exc_30,':','LineWidth',2)
xlim([0,720])
set(gca,'XTick',0:180:720)
grid on 
title('Porovnání vypočtených úhlu beta pro dva způsoby výpočtu')
xlabel('Natočení klikového hřídele [°]')
ylabel('Úhel odklonu ojnice [°]')

%% VÝPOČET ZTRÁT EXCENTRICKÉHO KM PRO CELÝ ROZSAH OTÁČEK

%% Ztrátový výkon centrický mechanismus, 1500-5500rpm + průměrná boční síla
for k=1:length(alfa)
PZCKM1500(k)=0.05*v(k,2)*Fn(k,2);
PZCKM1750(k)=0.05*v(k,3)*Fn(k,3);
PZCKM2000(k)=0.05*v(k,4)*Fn(k,4);
PZCKM2250(k)=0.05*v(k,5)*Fn(k,5);
PZCKM2500(k)=0.05*v(k,6)*Fn(k,6);
PZCKM2750(k)=0.05*v(k,7)*Fn(k,7);
PZCKM3000(k)=0.05*v(k,8)*Fn(k,8);
PZCKM3250(k)=0.05*v(k,9)*Fn(k,9);
PZCKM3500(k)=0.05*v(k,10)*Fn(k,10);
PZCKM3750(k)=0.05*v(k,11)*Fn(k,11);
PZCKM4000(k)=0.05*v(k,12)*Fn(k,12);
PZCKM4250(k)=0.05*v(k,13)*Fn(k,13);
PZCKM4500(k)=0.05*v(k,14)*Fn(k,14);
PZCKM4750(k)=0.05*v(k,15)*Fn(k,15);
PZCKM5000(k)=0.05*v(k,16)*Fn(k,16);
PZCKM5250(k)=0.05*v(k,17)*Fn(k,17);
PZCKM5500(k)=0.05*v(k,18)*Fn(k,18);
end
Pz_centricky1500rpm=mean(abs(PZCKM1500));
Pz_centricky1750rpm=mean(abs(PZCKM1750));
Pz_centricky2000rpm=mean(abs(PZCKM2000));
Pz_centricky2250rpm=mean(abs(PZCKM2250));
Pz_centricky2500rpm=mean(abs(PZCKM2500));
Pz_centricky2750rpm=mean(abs(PZCKM2750));
Pz_centricky3000rpm=mean(abs(PZCKM3000));
Pz_centricky3250rpm=mean(abs(PZCKM3250));
Pz_centricky3500rpm=mean(abs(PZCKM3500));
Pz_centricky3750rpm=mean(abs(PZCKM3750));
Pz_centricky4000rpm=mean(abs(PZCKM4000));
Pz_centricky4250rpm=mean(abs(PZCKM4250));
Pz_centricky4500rpm=mean(abs(PZCKM4500));
Pz_centricky4750rpm=mean(abs(PZCKM4750));
Pz_centricky5000rpm=mean(abs(PZCKM5000));
Pz_centricky5250rpm=mean(abs(PZCKM5250));
Pz_centricky5500rpm=mean(abs(PZCKM5500));


for i=2:1:18  
Fn_prumerna_centrickyKM(i)=mean(abs(Fn(:,i)));
end

%% Ztrátovy výkon 1500 rpm až po sílu v ojnici
for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n1500=1500;
        omega1500=2*pi*(n1500/60);
        
    for k=1:length(alfa)   
        v_e_1500(k,o)=r_excentricky*omega1500*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_1500(k,o)=r_excentricky*(omega1500^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end
nacteni_tlaky_excentricky_komplet_1500rpm=xlsread('tlaky_excentrickyKM_usporadane_1500rpm.xlsx');
tlaky_excentr_1500rpm_2az30=nacteni_tlaky_excentricky_komplet_1500rpm;
alfa_deg_exc=tlaky_excentr_1500rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;
ulozeni_tlaku_excentr_2az30mm_1500=tlaky_excentr_1500rpm_2az30(:,2:end);

for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_1500(k,o)=((ulozeni_tlaku_excentr_2az30mm_1500(k,o)*1e5-p0)*Spistu)-a_e_1500(k,o)*mpc;
            Fn_e_rozsah_1500(k,o)=Fnapistcelkova_exc_1500(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_1500(k,o)=Fn_e_rozsah_1500(k,o).*v_e_1500(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_1500(k,o)=Fnapistcelkova_exc_1500(k,o)/cos(beta_exc_rozsah(k,o));
           
    end                
end

for o=1:1:15
prumer_Fn_1500(o)=mean(abs(Fn_e_rozsah_1500(:,o)));
prumer_Pztratovy_1500(o)=mean(abs(Pztratovy_1500(:,o)));
end
%% Konec výpočtu ztrátový výkon 1500 rpm
%%

%% Ztrátovy výkon 1750 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n1750=1750;
        omega1750=2*pi*(n1750/60);
        
    for k=1:length(alfa)   
        v_e_1750(k,o)=r_excentricky*omega1750*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_1750(k,o)=r_excentricky*(omega1750^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_1750rpm=xlsread('tlaky_excentrickyKM_usporadane_1750rpm.xlsx');

tlaky_excentr_1750rpm_2az30=nacteni_tlaky_excentricky_komplet_1750rpm;
alfa_deg_exc=tlaky_excentr_1750rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_1750=tlaky_excentr_1750rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_1750(k,o)=((ulozeni_tlaku_excentr_2az30mm_1750(k,o)*1e5-p0)*Spistu)-a_e_1750(k,o)*mpc;
            Fn_e_rozsah_1750(k,o)=Fnapistcelkova_exc_1750(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_1750(k,o)=Fn_e_rozsah_1750(k,o).*v_e_1750(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_1750(k,o)=Fnapistcelkova_exc_1750(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_1750(o)=mean(abs(Fn_e_rozsah_1750(:,o)));
prumer_Pztratovy_1750(o)=mean(abs(Pztratovy_1750(:,o)));
end

%% Konec výpočtu ztrátový výkon 1750 rpm
%%

%% Ztrátový výkon 2000 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n2000=2000;
        omega2000=2*pi*(n2000/60);
        
    for k=1:length(alfa)   
        v_e_2000(k,o)=r_excentricky*omega2000*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_2000(k,o)=r_excentricky*(omega2000^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_2000rpm=xlsread('tlaky_excentrickyKM_usporadane_2000rpm.xlsx');

tlaky_excentr_2000rpm_2az30=nacteni_tlaky_excentricky_komplet_2000rpm;
alfa_deg_exc=tlaky_excentr_2000rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_2000=tlaky_excentr_2000rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_2000(k,o)=((ulozeni_tlaku_excentr_2az30mm_2000(k,o)*1e5-p0)*Spistu)-a_e_2000(k,o)*mpc;
            Fn_e_rozsah_2000(k,o)=Fnapistcelkova_exc_2000(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_2000(k,o)=Fn_e_rozsah_2000(k,o).*v_e_2000(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_2000(k,o)=Fnapistcelkova_exc_2000(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_2000(o)=mean(abs(Fn_e_rozsah_2000(:,o)));
prumer_Pztratovy_2000(o)=mean(abs(Pztratovy_2000(:,o)));
end

%% Konec výpočtu ztrátový výkon 2000 rpm
%%

%% Ztrátový výkon 2250 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n2250=2250;
        omega2250=2*pi*(n2250/60);
        
    for k=1:length(alfa)   
        v_e_2250(k,o)=r_excentricky*omega2250*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_2250(k,o)=r_excentricky*(omega2250^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_2250rpm=xlsread('tlaky_excentrickyKM_usporadane_2250rpm.xlsx');

tlaky_excentr_2250rpm_2az30=nacteni_tlaky_excentricky_komplet_2250rpm;
alfa_deg_exc=tlaky_excentr_2250rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_2250=tlaky_excentr_2250rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_2250(k,o)=((ulozeni_tlaku_excentr_2az30mm_2250(k,o)*1e5-p0)*Spistu)-a_e_2250(k,o)*mpc;
            Fn_e_rozsah_2250(k,o)=Fnapistcelkova_exc_2250(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_2250(k,o)=Fn_e_rozsah_2250(k,o).*v_e_2250(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_2250(k,o)=Fnapistcelkova_exc_2250(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_2250(o)=mean(abs(Fn_e_rozsah_2250(:,o)));
prumer_Pztratovy_2250(o)=mean(abs(Pztratovy_2250(:,o)));
end

%% Konec výpočtu ztrátový výkon 2250 rpm
%%

%% Ztrátový výkon 2500 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n2500=2500;
        omega2500=2*pi*(n2500/60);
        
    for k=1:length(alfa)   
        v_e_2500(k,o)=r_excentricky*omega2500*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_2500(k,o)=r_excentricky*(omega2500^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_2500rpm=xlsread('tlaky_excentrickyKM_usporadane_2500rpm.xlsx');

tlaky_excentr_2500rpm_2az30=nacteni_tlaky_excentricky_komplet_2500rpm;
alfa_deg_exc=tlaky_excentr_2500rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_2500=tlaky_excentr_2500rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_2500(k,o)=((ulozeni_tlaku_excentr_2az30mm_2500(k,o)*1e5-p0)*Spistu)-a_e_2500(k,o)*mpc;
            Fn_e_rozsah_2500(k,o)=Fnapistcelkova_exc_2500(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_2500(k,o)=Fn_e_rozsah_2500(k,o).*v_e_2500(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_2500(k,o)=Fnapistcelkova_exc_2500(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_2500(o)=mean(abs(Fn_e_rozsah_2500(:,o)));
prumer_Pztratovy_2500(o)=mean(abs(Pztratovy_2500(:,o)));
end

%% Konec výpočtu ztrátový výkon 2500 rpm
%%

%% Ztrátový výkon 2750 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n2750=2750;
        omega2750=2*pi*(n2750/60);
        
    for k=1:length(alfa)   
        v_e_2750(k,o)=r_excentricky*omega2750*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_2750(k,o)=r_excentricky*(omega2750^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_2750rpm=xlsread('tlaky_excentrickyKM_usporadane_2750rpm.xlsx');

tlaky_excentr_2750rpm_2az30=nacteni_tlaky_excentricky_komplet_2750rpm;
alfa_deg_exc=tlaky_excentr_2750rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_2750=tlaky_excentr_2750rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_2750(k,o)=((ulozeni_tlaku_excentr_2az30mm_2750(k,o)*1e5-p0)*Spistu)-a_e_2750(k,o)*mpc;
            Fn_e_rozsah_2750(k,o)=Fnapistcelkova_exc_2750(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_2750(k,o)=Fn_e_rozsah_2750(k,o).*v_e_2750(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_2750(k,o)=Fnapistcelkova_exc_2750(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_2750(o)=mean(abs(Fn_e_rozsah_2750(:,o)));
prumer_Pztratovy_2750(o)=mean(abs(Pztratovy_2750(:,o)));
end

%% Konec výpočtu ztrátový výkon 2750 rpm
%%

%% Ztrátový výkon 3000 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n3000=3000;
        omega3000=2*pi*(n3000/60);
        
    for k=1:length(alfa)   
        v_e_3000(k,o)=r_excentricky*omega3000*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_3000(k,o)=r_excentricky*(omega3000^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_3000rpm=xlsread('tlaky_excentrickyKM_usporadane_3000rpm.xlsx');

tlaky_excentr_3000rpm_2az30=nacteni_tlaky_excentricky_komplet_3000rpm;
alfa_deg_exc=tlaky_excentr_3000rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_3000=tlaky_excentr_3000rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_3000(k,o)=((ulozeni_tlaku_excentr_2az30mm_3000(k,o)*1e5-p0)*Spistu)-a_e_3000(k,o)*mpc;
            Fn_e_rozsah_3000(k,o)=Fnapistcelkova_exc_3000(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_3000(k,o)=Fn_e_rozsah_3000(k,o).*v_e_3000(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_3000(k,o)=Fnapistcelkova_exc_3000(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_3000(o)=mean(abs(Fn_e_rozsah_3000(:,o)));
prumer_Pztratovy_3000(o)=mean(abs(Pztratovy_3000(:,o)));
end

%% Konec výpočtu ztrátový výkon 3000 rpm
%%


%% Ztrátový výkon 3250 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n3250=3250;
        omega3250=2*pi*(n3250/60);
        
    for k=1:length(alfa)   
        v_e_3250(k,o)=r_excentricky*omega3250*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_3250(k,o)=r_excentricky*(omega3250^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_3250rpm=xlsread('tlaky_excentrickyKM_usporadane_3250rpm.xlsx');

tlaky_excentr_3250rpm_2az30=nacteni_tlaky_excentricky_komplet_3250rpm;
alfa_deg_exc=tlaky_excentr_3250rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_3250=tlaky_excentr_3250rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_3250(k,o)=((ulozeni_tlaku_excentr_2az30mm_3250(k,o)*1e5-p0)*Spistu)-a_e_3250(k,o)*mpc;
            Fn_e_rozsah_3250(k,o)=Fnapistcelkova_exc_3250(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_3250(k,o)=Fn_e_rozsah_3250(k,o).*v_e_3250(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_3250(k,o)=Fnapistcelkova_exc_3250(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_3250(o)=mean(abs(Fn_e_rozsah_3250(:,o)));
prumer_Pztratovy_3250(o)=mean(abs(Pztratovy_3250(:,o)));
end

%% Konec výpočtu ztrátový výkon 3250 rpm
%%

%% Ztrátový výkon 3500 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n3500=3500;
        omega3500=2*pi*(n3500/60);
        
    for k=1:length(alfa)   
        v_e_3500(k,o)=r_excentricky*omega3500*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_3500(k,o)=r_excentricky*(omega3500^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_3500rpm=xlsread('tlaky_excentrickyKM_usporadane_3500rpm.xlsx');

tlaky_excentr_3500rpm_2az30=nacteni_tlaky_excentricky_komplet_3500rpm;
alfa_deg_exc=tlaky_excentr_3500rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_3500=tlaky_excentr_3500rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_3500(k,o)=((ulozeni_tlaku_excentr_2az30mm_3500(k,o)*1e5-p0)*Spistu)-a_e_3500(k,o)*mpc;
            Fn_e_rozsah_3500(k,o)=Fnapistcelkova_exc_3500(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_3500(k,o)=Fn_e_rozsah_3500(k,o).*v_e_3500(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_3500(k,o)=Fnapistcelkova_exc_3500(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_3500(o)=mean(abs(Fn_e_rozsah_3500(:,o)));
prumer_Pztratovy_3500(o)=mean(abs(Pztratovy_3500(:,o)));
end

%% Konec výpočtu ztrátový výkon 3500 rpm
%%


%% Ztrátový výkon 3750 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n3750=3750;
        omega3750=2*pi*(n3750/60);
        
    for k=1:length(alfa)   
        v_e_3750(k,o)=r_excentricky*omega3750*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_3750(k,o)=r_excentricky*(omega3750^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_3750rpm=xlsread('tlaky_excentrickyKM_usporadane_3750rpm.xlsx');

tlaky_excentr_3750rpm_2az30=nacteni_tlaky_excentricky_komplet_3750rpm;
alfa_deg_exc=tlaky_excentr_3750rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_3750=tlaky_excentr_3750rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_3750(k,o)=((ulozeni_tlaku_excentr_2az30mm_3750(k,o)*1e5-p0)*Spistu)-a_e_3750(k,o)*mpc;
            Fn_e_rozsah_3750(k,o)=Fnapistcelkova_exc_3750(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_3750(k,o)=Fn_e_rozsah_3750(k,o).*v_e_3750(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_3750(k,o)=Fnapistcelkova_exc_3750(k,o)/cos(beta_exc_rozsah(k,o));
            
            
    end                
end

for o=1:1:15
prumer_Fn_3750(o)=mean(abs(Fn_e_rozsah_3750(:,o)));
prumer_Pztratovy_3750(o)=mean(abs(Pztratovy_3750(:,o)));
end

%% Konec výpočtu ztrátový výkon 3750 rpm
%%

%% Ztrátový výkon 4000 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n4000=4000;
        omega4000=2*pi*(n4000/60);
        
    for k=1:length(alfa)   
        v_e_4000(k,o)=r_excentricky*omega4000*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_4000(k,o)=r_excentricky*(omega4000^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_4000rpm=xlsread('tlaky_excentrickyKM_usporadane_4000rpm.xlsx');

tlaky_excentr_4000rpm_2az30=nacteni_tlaky_excentricky_komplet_4000rpm;
alfa_deg_exc=tlaky_excentr_4000rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_4000=tlaky_excentr_4000rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_4000(k,o)=((ulozeni_tlaku_excentr_2az30mm_4000(k,o)*1e5-p0)*Spistu)-a_e_4000(k,o)*mpc;
            Fn_e_rozsah_4000(k,o)=Fnapistcelkova_exc_4000(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_4000(k,o)=Fn_e_rozsah_4000(k,o).*v_e_4000(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_4000(k,o)=Fnapistcelkova_exc_4000(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_4000(o)=mean(abs(Fn_e_rozsah_4000(:,o)));
prumer_Pztratovy_4000(o)=mean(abs(Pztratovy_4000(:,o)));
end

%% Konec výpočtu ztrátový výkon 4000 rpm
%%

%% Ztrátový výkon 4250 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n4250=4250;
        omega4250=2*pi*(n4250/60);
        
    for k=1:length(alfa)   
        v_e_4250(k,o)=r_excentricky*omega4250*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_4250(k,o)=r_excentricky*(omega4250^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_4250rpm=xlsread('tlaky_excentrickyKM_usporadane_4250rpm.xlsx');

tlaky_excentr_4250rpm_2az30=nacteni_tlaky_excentricky_komplet_4250rpm;
alfa_deg_exc=tlaky_excentr_4250rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_4250=tlaky_excentr_4250rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_4250(k,o)=((ulozeni_tlaku_excentr_2az30mm_4250(k,o)*1e5-p0)*Spistu)-a_e_4250(k,o)*mpc;
            Fn_e_rozsah_4250(k,o)=Fnapistcelkova_exc_4250(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_4250(k,o)=Fn_e_rozsah_4250(k,o).*v_e_4250(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_4250(k,o)=Fnapistcelkova_exc_4250(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_4250(o)=mean(abs(Fn_e_rozsah_4250(:,o)));
prumer_Pztratovy_4250(o)=mean(abs(Pztratovy_4250(:,o)));
end

%% Konec výpočtu ztrátový výkon 4250 rpm
%%

%% Ztrátový výkon 4500 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n4500=4500;
        omega4500=2*pi*(n4500/60);
        
    for k=1:length(alfa)   
        v_e_4500(k,o)=r_excentricky*omega4500*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_4500(k,o)=r_excentricky*(omega4500^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_4500rpm=xlsread('tlaky_excentrickyKM_usporadane_4500rpm.xlsx');

tlaky_excentr_4500rpm_2az30=nacteni_tlaky_excentricky_komplet_4500rpm;
alfa_deg_exc=tlaky_excentr_4500rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_4500=tlaky_excentr_4500rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_4500(k,o)=((ulozeni_tlaku_excentr_2az30mm_4500(k,o)*1e5-p0)*Spistu)-a_e_4500(k,o)*mpc;
            Fn_e_rozsah_4500(k,o)=Fnapistcelkova_exc_4500(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_4500(k,o)=Fn_e_rozsah_4500(k,o).*v_e_4500(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_4500(k,o)=Fnapistcelkova_exc_4500(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_4500(o)=mean(abs(Fn_e_rozsah_4500(:,o)));
prumer_Pztratovy_4500(o)=mean(abs(Pztratovy_4500(:,o)));
end

%% Konec výpočtu ztrátový výkon 4500 rpm
%%

%% Ztrátový výkon 4750 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n4750=4750;
        omega4750=2*pi*(n4750/60);
        
    for k=1:length(alfa)   
        v_e_4750(k,o)=r_excentricky*omega4750*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_4750(k,o)=r_excentricky*(omega4750^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_4750rpm=xlsread('tlaky_excentrickyKM_usporadane_4750rpm.xlsx');

tlaky_excentr_4750rpm_2az30=nacteni_tlaky_excentricky_komplet_4750rpm;
alfa_deg_exc=tlaky_excentr_4750rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_4750=tlaky_excentr_4750rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_4750(k,o)=((ulozeni_tlaku_excentr_2az30mm_4750(k,o)*1e5-p0)*Spistu)-a_e_4750(k,o)*mpc;
            Fn_e_rozsah_4750(k,o)=Fnapistcelkova_exc_4750(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_4750(k,o)=Fn_e_rozsah_4750(k,o).*v_e_4750(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_4750(k,o)=Fnapistcelkova_exc_4750(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_4750(o)=mean(abs(Fn_e_rozsah_4750(:,o)));
prumer_Pztratovy_4750(o)=mean(abs(Pztratovy_4750(:,o)));
end

%% Konec výpočtu ztrátový výkon 4750 rpm
%%

%% Ztrátový výkon 5000 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n5000=5000;
        omega5000=2*pi*(n5000/60);
        
    for k=1:length(alfa)   
        v_e_5000(k,o)=r_excentricky*omega5000*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_5000(k,o)=r_excentricky*(omega5000^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_5000rpm=xlsread('tlaky_excentrickyKM_usporadane_5000rpm.xlsx');

tlaky_excentr_5000rpm_2az30=nacteni_tlaky_excentricky_komplet_5000rpm;
alfa_deg_exc=tlaky_excentr_5000rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_5000=tlaky_excentr_5000rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_5000(k,o)=((ulozeni_tlaku_excentr_2az30mm_5000(k,o)*1e5-p0)*Spistu)-a_e_5000(k,o)*mpc;
            Fn_e_rozsah_5000(k,o)=Fnapistcelkova_exc_5000(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_5000(k,o)=Fn_e_rozsah_5000(k,o).*v_e_5000(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_5000(k,o)=Fnapistcelkova_exc_5000(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_5000(o)=mean(abs(Fn_e_rozsah_5000(:,o)));
prumer_Pztratovy_5000(o)=mean(abs(Pztratovy_5000(:,o)));
end

%% Konec výpočtu ztrátový výkon 5000 rpm
%%

%% Ztrátový výkon 5250 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n5250=5250;
        omega5250=2*pi*(n5250/60);
        
    for k=1:length(alfa)   
        v_e_5250(k,o)=r_excentricky*omega5250*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_5250(k,o)=r_excentricky*(omega5250^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_5250rpm=xlsread('tlaky_excentrickyKM_usporadane_5250rpm.xlsx');

tlaky_excentr_5250rpm_2az30=nacteni_tlaky_excentricky_komplet_5250rpm;
alfa_deg_exc=tlaky_excentr_5250rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_5250=tlaky_excentr_5250rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_5250(k,o)=((ulozeni_tlaku_excentr_2az30mm_5250(k,o)*1e5-p0)*Spistu)-a_e_5250(k,o)*mpc;
            Fn_e_rozsah_5250(k,o)=Fnapistcelkova_exc_5250(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_5250(k,o)=Fn_e_rozsah_5250(k,o).*v_e_5250(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_5250(k,o)=Fnapistcelkova_exc_5250(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_5250(o)=mean(abs(Fn_e_rozsah_5250(:,o)));
prumer_Pztratovy_5250(o)=mean(abs(Pztratovy_5250(:,o)));
end

%% Konec výpočtu ztrátový výkon 5250 rpm
%%

%% Ztrátový výkon 5500 rpm

for o=1:1:15
        e_rozsah(o)=(o*2)/1000;
        lambda_e_rozsah(o)=e_rozsah(o)/l_exentricky;
        n5500=5500;
        omega5500=2*pi*(n5500/60);
        
    for k=1:length(alfa)   
        v_e_5500(k,o)=r_excentricky*omega5500*(sin(alfa(k))-lambda_e_rozsah(o)*cos(alfa(k))+((lambda_pro_excentricky/2)*sin(2*alfa(k))));
        a_e_5500(k,o)=r_excentricky*(omega5500^2)*(cos(alfa(k))+lambda_e_rozsah(o)*sin(alfa(k))+((lambda_pro_excentricky)*cos(2*alfa(k))));
    end
end

nacteni_tlaky_excentricky_komplet_5500rpm=xlsread('tlaky_excentrickyKM_usporadane_5500rpm.xlsx');

tlaky_excentr_5500rpm_2az30=nacteni_tlaky_excentricky_komplet_5500rpm;
alfa_deg_exc=tlaky_excentr_5500rpm_2az30(:,1);
alfa_exc=(alfa_deg_exc/180)*pi;

ulozeni_tlaku_excentr_2az30mm_5500=tlaky_excentr_5500rpm_2az30(:,2:end);


for k=1:length(alfa_exc)
    for o=1:1:15
            lambda_e_rozsah(o)=(0.002*o)/l_exentricky;    
            beta_exc_rozsah(k,o)=asin(lambda_pro_excentricky*sin(alfa_exc(k))-lambda_e_rozsah(o));

            Fnapistcelkova_exc_5500(k,o)=((ulozeni_tlaku_excentr_2az30mm_5500(k,o)*1e5-p0)*Spistu)-a_e_5500(k,o)*mpc;
            Fn_e_rozsah_5500(k,o)=Fnapistcelkova_exc_5500(k,o).*tan(beta_exc_rozsah(k,o));

            Pztratovy_5500(k,o)=Fn_e_rozsah_5500(k,o).*v_e_5500(k,o)*koeficient_treni_prumerny;
            
            Fo_exc_5500(k,o)=Fnapistcelkova_exc_5500(k,o)/cos(beta_exc_rozsah(k,o));
    end                
end

for o=1:1:15
prumer_Fn_5500(o)=mean(abs(Fn_e_rozsah_5500(:,o)));
prumer_Pztratovy_5500(o)=mean(abs(Pztratovy_5500(:,o)));
end

%% Konec výpočtu ztrátový výkon 5500 rpm
%% definování rozsahu excentricit pro grafy
excentricity=2:2:30;
excentricity_0az30=0:2:30;

%% graf závislosti síly působící na píst na excentricitě pro 5500rpm

figure(31)  
hold on
plot(alfa_deg,Fvp(:,18),'k--','LineWidth',2)
plot(alfa_deg,Fnapistcelkova_exc_5500(:,2),'g-','LineWidth',2)
plot(alfa_deg,Fnapistcelkova_exc_5500(:,4),'b-','LineWidth',2)
plot(alfa_deg,Fnapistcelkova_exc_5500(:,6),'y-','LineWidth',2)
plot(alfa_deg,Fnapistcelkova_exc_5500(:,8),'m-','LineWidth',2)
plot(alfa_deg,Fnapistcelkova_exc_5500(:,10),'c-','LineWidth',2)
plot(alfa_deg,Fnapistcelkova_exc_5500(:,15),'r-','LineWidth',2)
hold off
title('Průběh výsledné síly na píst při 5500min-1')
xlabel('Natočení klikového hřídele [°]')
ylabel('Výsledná síla působící na píst [N]')
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
legend('e=0mm','e=4mm','e=8mm','e=12mm','e=16mm','e=20mm','e=30mm')
grid on

figure(32)
hold on
plot(alfa_deg,Fnapistcelkova_exc_5500(:,10),'r','LineWidth',3)
plot(alfa_deg,((ulozeni_tlaku_excentr_2az30mm_5500(:,10)*1e5-p0)*Spistu),'b-','LineWidth',2)
plot(alfa_deg,-a_e_5500(:,10).*mpc,'g-','LineWidth',2)
plot(alfa_deg,Fnapistcelkova_exc_5500(:,5),'k','LineWidth',3)
plot(alfa_deg,((ulozeni_tlaku_excentr_2az30mm_5500(:,5)*1e5-p0)*Spistu),'m-','LineWidth',2)
plot(alfa_deg,-a_e_5500(:,5).*mpc,'c-','LineWidth',2)
hold off
title('Průběh složek výsledné síly na píst při 5500min-1')
xlabel('Natočení klikového hřídele [°]')
ylabel('Síla působící na píst [N]')
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
legend('Výsledná síla na píst e=20mm','Síla od tlaku plynů e=20mm','Setrvačná složka síly e=20mm','Výsledná síla na píst e=10mm','Síla od tlaku plynů e=10mm','Setrvačná složka síly e=10mm')

grid on

%% graf závislosti boční síly působící na píst na excentricitě pro 5500rpm

figure(34)  
hold on
plot(alfa_deg,Fn(:,18),'k--','LineWidth',2)
plot(alfa_deg,Fn_e_rozsah_5500(:,2),'g-','LineWidth',2)
plot(alfa_deg,Fn_e_rozsah_5500(:,4),'b-','LineWidth',2)
plot(alfa_deg,Fn_e_rozsah_5500(:,6),'y-','LineWidth',2)
plot(alfa_deg,Fn_e_rozsah_5500(:,8),'m-','LineWidth',2)
plot(alfa_deg,Fn_e_rozsah_5500(:,10),'c-','LineWidth',2)
plot(alfa_deg,Fn_e_rozsah_5500(:,15),'r-','LineWidth',2)

% plot(alfa_deg,mean(abs(Fn(:,18))),'k--','LineWidth',2)
% plot(alfa_deg,prumer_Fn_5500(2),'g--','LineWidth',1)
% plot(alfa_deg,prumer_Fn_5500(4),'b.','MarkerSize',4)
% plot(alfa_deg,prumer_Fn_5500(6),'y--','LineWidth',1)
% plot(alfa_deg,prumer_Fn_5500(8),'m--','LineWidth',5)
% plot(0:719,prumer_Fn_5500(10),'c--','LineWidth',1)
% plot(alfa_deg,prumer_Fn_5500(15),'r--','LineWidth',1)

hold off
title('Průběh výsledné boční síly na píst při 5500min-1')
xlabel('Natočení klikového hřídele [°]')
ylabel('Výsledná boční síla působící na píst [N]')
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
legend('e=0mm','e=4mm','e=8mm','e=12mm','e=16mm','e=20mm','e=30mm')
grid on

%% 3D graf rpm-e-Pztratovy
rpm_rozsah=[1500,1750,2000,2250,2500,2750,3000,3250,3500,3750,4000,4250,4500,4750,5000,5250,5500];
matice_ztratovych_vykonu_centr=[Pz_centricky1500rpm;Pz_centricky1750rpm;Pz_centricky2000rpm;Pz_centricky2250rpm;Pz_centricky2500rpm;Pz_centricky2750rpm;Pz_centricky3000rpm;Pz_centricky3250rpm;Pz_centricky3500rpm;Pz_centricky3750rpm;Pz_centricky4000rpm;Pz_centricky4250rpm;Pz_centricky4500rpm;Pz_centricky4750rpm;Pz_centricky5000rpm;Pz_centricky5250rpm;Pz_centricky5500rpm];
matice_ztratovych_vykonu_excentr=[prumer_Pztratovy_1500;prumer_Pztratovy_1750;prumer_Pztratovy_2000;prumer_Pztratovy_2250;prumer_Pztratovy_2500;prumer_Pztratovy_2750;prumer_Pztratovy_3000;prumer_Pztratovy_3250;prumer_Pztratovy_3500;prumer_Pztratovy_3750;prumer_Pztratovy_4000;prumer_Pztratovy_4250;prumer_Pztratovy_4500;prumer_Pztratovy_4750;prumer_Pztratovy_5000;prumer_Pztratovy_5250;prumer_Pztratovy_5500];

matice_ztratova_celkova=[matice_ztratovych_vykonu_centr,matice_ztratovych_vykonu_excentr];

figure(35)
surf(excentricity_0az30,rpm_rozsah,matice_ztratova_celkova,'MarkerSize',1)
colorbar
title('Ztrátový výkon v závislosti na otáčkách a excentricitě')
xlabel('Excentricita [mm]')
ylabel('Otáčky motoru [min-1]')
zlabel('Ztrátový výkon v kontaktu vložky válce s pístem [W]')
set(gca,'XTick',0:2:30)
set(gca,'ZTick',0:100:800)
set(gca,'YTick',1500:500:5500)
grid on

%% grafy průběhů třecích ztrát na excentricitě pro jednotlivé otáčky
figure(36)  
hold on
plot(0,Pz_centricky1500rpm,'g+','MarkerSize',13,'LineWidth',4)
plot(excentricity,prumer_Pztratovy_1500,'g-o','LineWidth',2)
plot(0,Pz_centricky2500rpm,'m+','MarkerSize',13,'LineWidth',4)
plot(excentricity,prumer_Pztratovy_2500,'m-o','LineWidth',2)
plot(0,Pz_centricky3750rpm,'b+','MarkerSize',13,'LineWidth',4)
plot(excentricity,prumer_Pztratovy_3750,'b-o','LineWidth',2)
plot(0,Pz_centricky4500rpm,'c+','MarkerSize',13,'LineWidth',4)
plot(excentricity,prumer_Pztratovy_4500,'c-o','LineWidth',2)
plot(0,Pz_centricky5500rpm,'r+','MarkerSize',13,'LineWidth',4)
plot(excentricity,prumer_Pztratovy_5500,'r-o','LineWidth',2)
hold off
title('Třecí výkon ve styku pístu s vložkou válce')
xlabel('Excentricita [mm]')
ylabel('Ztrátový výkon [W]')
set(gca,'XTick',0:2:30)
grid on 
legend('Centrický KM 1500min-1','Excentrický KM 1500min-1','Centrický KM 2500min-1','Excentrický KM 2500min-1','Centrický KM 3750min-1','Excentrický KM 3750min-1','Centrický KM 4500min-1','Excentrický KM 4500min-1','Centrický KM 5500min-1','Excentrický KM 5500min-1')                
      
%% graf síly působící v ojnici pro 1500min-1
figure(37)  
hold on
plot(alfa_deg,Fo(:,2),'k--','LineWidth',2)
plot(alfa_deg,Fo_exc_1500(:,2),'g-','LineWidth',2)
plot(alfa_deg,Fo_exc_1500(:,4),'b-','LineWidth',2)
plot(alfa_deg,Fo_exc_1500(:,6),'y-','LineWidth',2)
plot(alfa_deg,Fo_exc_1500(:,8),'m-','LineWidth',2)
plot(alfa_deg,Fo_exc_1500(:,10),'c-','LineWidth',2)
plot(alfa_deg,Fo_exc_1500(:,15),'r-','LineWidth',2)

hold off
title('Průběh výsledné síly psůobící v ojnici pro 1500min-1')
xlabel('Natočení klikového hřídele [°]')
ylabel('Výsledná síla působící v ojnici [N]')
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
legend('e=0mm','e=4mm','e=8mm','e=12mm','e=16mm','e=20mm','e=30mm')
grid on                
               
%% graf síly působící v ojnici celý rozsah otáček a volenou excentricita
figure(38)  
hold on
plot(alfa_deg,Fo(:,18),'k--','LineWidth',2)
plot(alfa_deg,Fo_exc_1500(:,7),'g-','LineWidth',2)
plot(alfa_deg,Fo_exc_1750(:,7),'b-','LineWidth',2)
plot(alfa_deg,Fo_exc_2000(:,7),'r-','LineWidth',2)
plot(alfa_deg,Fo_exc_2250(:,7),'c-','LineWidth',2)
plot(alfa_deg,Fo_exc_2500(:,7),'b-','LineWidth',2)
plot(alfa_deg,Fo_exc_2750(:,7),'y-','LineWidth',2)
plot(alfa_deg,Fo_exc_3000(:,7),'m-','LineWidth',2)
plot(alfa_deg,Fo_exc_3250(:,7),'g-','LineWidth',2)
plot(alfa_deg,Fo_exc_3500(:,7),'b-','LineWidth',2)
plot(alfa_deg,Fo_exc_3750(:,7),'r-','LineWidth',2)
plot(alfa_deg,Fo_exc_4000(:,7),'c-','LineWidth',2)
plot(alfa_deg,Fo_exc_4250(:,7),'k-','LineWidth',2)
plot(alfa_deg,Fo_exc_4500(:,7),'y-','LineWidth',2)
plot(alfa_deg,Fo_exc_4750(:,7),'m-','LineWidth',2)
plot(alfa_deg,Fo_exc_5000(:,7),'g-','LineWidth',2)
plot(alfa_deg,Fo_exc_5250(:,7),'b-','LineWidth',2)
plot(alfa_deg,Fo_exc_5500(:,7),'r-','LineWidth',2)

hold off
title('Průběh výsledné síly působící v ojnici pro celý rozsah provoních otáček')
xlabel('Natočení klikového hřídele [°]')
ylabel('Výsledná síla působící v ojnici [N]')
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
legend('Centrický 5500min-1','1500min-1','1750min-1','2000min-1','2250min-1','2500min-1','2750min-1','3000min-1','3250min-1','3500min-1','3750min-1','4000min-1','4250min-1','4500min-1','4750min-1','5000min-1','5250min-1','5500min-1')
grid on   
                
%% graf síly působící v ojnici - vybrané otáčky a volená excentricita

figure(39)  
hold on
% plot(alfa_deg,Fo(:,2),'k--','LineWidth',2)
% plot(alfa_deg,Fo_exc_1500(:,7),'k-','LineWidth',2)
% plot(alfa_deg,Fo(:,6),'k--','LineWidth',2)
% plot(alfa_deg,Fo_exc_2500(:,7),'k-','LineWidth',2)
plot(alfa_deg,Fo(:,11),'b--','LineWidth',2)
plot(alfa_deg,Fo_exc_3750(:,7),'b-','LineWidth',2)
plot(alfa_deg,Fo(:,13),'k--','LineWidth',2)
plot(alfa_deg,Fo_exc_4250(:,7),'k-','LineWidth',2)
plot(alfa_deg,Fo(:,18),'r--','LineWidth',2)
plot(alfa_deg,Fo_exc_5500(:,7),'r-','LineWidth',2)

hold off
title('Průběh výsledné síly působící v ojnici pro vybrané otáčky a zvolenou excentricitu 14mm')
xlabel('Natočení klikového hřídele [°]')
ylabel('Výsledná síla působící v ojnici [N]')
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
% legend('CKM 1500min-1','KM 1500min-1','CKM 2500min-1','EKM 2500min-1','CKM 3750min-1','EKM 3750min-1','CKM 4500min-1','EKM 4500min-1','CKM 5500min-1','EKM 5500min-1')
legend('CKM 3750min-1','EKM 3750min-1','CKM 4500min-1','EKM 4500min-1','CKM 5500min-1','EKM 5500min-1')

grid on   

%% Výpočty sil v ojnici a sil na ojniční čep; točivý moment

for k=1:length(alfa_exc)
    for o=1:1:15
            Ft_exc_1500(k,o)=Fo_exc_1500(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_1500(k,o)=Fo_exc_1500(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_1500(k,o)=-m2*r*omega1500^2;        
            Frc_exc_1500(k,o)=Fr_exc_1500(k,o)+Fro_exc_1500(k,o);
            Fvoc_exc_1500(k,o)=sqrt(Frc_exc_1500(k,o)^2+Ft_exc_1500(k,o)^2);
            Mtoc_exc_1500(k,o)=Ft_exc_1500(k,o)*r;
     
            Ft_exc_1750(k,o)=Fo_exc_1750(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_1750(k,o)=Fo_exc_1750(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_1750(k,o)=-m2*r*omega1750^2;        
            Frc_exc_1750(k,o)=Fr_exc_1750(k,o)+Fro_exc_1750(k,o);
            Fvoc_exc_1750(k,o)=sqrt(Frc_exc_1750(k,o)^2+Ft_exc_1750(k,o)^2);
            Mtoc_exc_1750(k,o)=Ft_exc_1750(k,o)*r;
            
            Ft_exc_2000(k,o)=Fo_exc_2000(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_2000(k,o)=Fo_exc_2000(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_2000(k,o)=-m2*r*omega2000^2;        
            Frc_exc_2000(k,o)=Fr_exc_2000(k,o)+Fro_exc_2000(k,o);
            Fvoc_exc_2000(k,o)=sqrt(Frc_exc_2000(k,o)^2+Ft_exc_2000(k,o)^2);
            Mtoc_exc_2000(k,o)=Ft_exc_2000(k,o)*r;
            
            Ft_exc_2250(k,o)=Fo_exc_2250(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_2250(k,o)=Fo_exc_2250(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_2250(k,o)=-m2*r*omega2250^2;        
            Frc_exc_2250(k,o)=Fr_exc_2250(k,o)+Fro_exc_2250(k,o);
            Fvoc_exc_2250(k,o)=sqrt(Frc_exc_2250(k,o)^2+Ft_exc_2250(k,o)^2);
            Mtoc_exc_2250(k,o)=Ft_exc_2250(k,o)*r;
            
            Ft_exc_2500(k,o)=Fo_exc_2500(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_2500(k,o)=Fo_exc_2500(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_2500(k,o)=-m2*r*omega2500^2;        
            Frc_exc_2500(k,o)=Fr_exc_2500(k,o)+Fro_exc_2500(k,o);
            Fvoc_exc_2500(k,o)=sqrt(Frc_exc_2500(k,o)^2+Ft_exc_2500(k,o)^2);
            Mtoc_exc_2500(k,o)=Ft_exc_2500(k,o)*r;
                      
            Ft_exc_2750(k,o)=Fo_exc_2750(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_2750(k,o)=Fo_exc_2750(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_2750(k,o)=-m2*r*omega2750^2;        
            Frc_exc_2750(k,o)=Fr_exc_2750(k,o)+Fro_exc_2750(k,o);
            Fvoc_exc_2750(k,o)=sqrt(Frc_exc_2750(k,o)^2+Ft_exc_2750(k,o)^2);
            Mtoc_exc_2750(k,o)=Ft_exc_2750(k,o)*r;
                                    
            Ft_exc_3000(k,o)=Fo_exc_3000(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_3000(k,o)=Fo_exc_3000(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_3000(k,o)=-m2*r*omega3000^2;        
            Frc_exc_3000(k,o)=Fr_exc_3000(k,o)+Fro_exc_3000(k,o);
            Fvoc_exc_3000(k,o)=sqrt(Frc_exc_3000(k,o)^2+Ft_exc_3000(k,o)^2);
            Mtoc_exc_3000(k,o)=Ft_exc_3000(k,o)*r;
                                    
            Ft_exc_3250(k,o)=Fo_exc_3250(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_3250(k,o)=Fo_exc_3250(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_3250(k,o)=-m2*r*omega3250^2;        
            Frc_exc_3250(k,o)=Fr_exc_3250(k,o)+Fro_exc_3250(k,o);
            Fvoc_exc_3250(k,o)=sqrt(Frc_exc_3250(k,o)^2+Ft_exc_3250(k,o)^2);
            Mtoc_exc_3250(k,o)=Ft_exc_3250(k,o)*r;
                                                
            Ft_exc_3500(k,o)=Fo_exc_3500(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_3500(k,o)=Fo_exc_3500(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_3500(k,o)=-m2*r*omega3500^2;        
            Frc_exc_3500(k,o)=Fr_exc_3500(k,o)+Fro_exc_3500(k,o);
            Fvoc_exc_3500(k,o)=sqrt(Frc_exc_3500(k,o)^2+Ft_exc_3500(k,o)^2);
            Mtoc_exc_3500(k,o)=Ft_exc_3500(k,o)*r;
                         
            Ft_exc_3750(k,o)=Fo_exc_3750(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_3750(k,o)=Fo_exc_3750(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_3750(k,o)=-m2*r*omega3750^2;        
            Frc_exc_3750(k,o)=Fr_exc_3750(k,o)+Fro_exc_3750(k,o);
            Fvoc_exc_3750(k,o)=sqrt(Frc_exc_3750(k,o)^2+Ft_exc_3750(k,o)^2);
            Mtoc_exc_3750(k,o)=Ft_exc_3750(k,o)*r;
                                    
            Ft_exc_4000(k,o)=Fo_exc_4000(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_4000(k,o)=Fo_exc_4000(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_4000(k,o)=-m2*r*omega4000^2;        
            Frc_exc_4000(k,o)=Fr_exc_4000(k,o)+Fro_exc_4000(k,o);
            Fvoc_exc_4000(k,o)=sqrt(Frc_exc_4000(k,o)^2+Ft_exc_4000(k,o)^2);
            Mtoc_exc_4000(k,o)=Ft_exc_4000(k,o)*r;
                                    
            Ft_exc_4250(k,o)=Fo_exc_4250(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_4250(k,o)=Fo_exc_4250(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_4250(k,o)=-m2*r*omega4250^2;        
            Frc_exc_4250(k,o)=Fr_exc_4250(k,o)+Fro_exc_4250(k,o);
            Fvoc_exc_4250(k,o)=sqrt(Frc_exc_4250(k,o)^2+Ft_exc_4250(k,o)^2);
            Mtoc_exc_4250(k,o)=Ft_exc_4250(k,o)*r;
                                    
            Ft_exc_4500(k,o)=Fo_exc_4500(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_4500(k,o)=Fo_exc_4500(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_4500(k,o)=-m2*r*omega4500^2;        
            Frc_exc_4500(k,o)=Fr_exc_4500(k,o)+Fro_exc_4500(k,o);
            Fvoc_exc_4500(k,o)=sqrt(Frc_exc_4500(k,o)^2+Ft_exc_4500(k,o)^2);
            Mtoc_exc_4500(k,o)=Ft_exc_4500(k,o)*r;
            
            Ft_exc_4750(k,o)=Fo_exc_4750(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_4750(k,o)=Fo_exc_4750(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_4750(k,o)=-m2*r*omega4750^2;        
            Frc_exc_4750(k,o)=Fr_exc_4750(k,o)+Fro_exc_4750(k,o);
            Fvoc_exc_4750(k,o)=sqrt(Frc_exc_4750(k,o)^2+Ft_exc_4750(k,o)^2);
            Mtoc_exc_4750(k,o)=Ft_exc_4750(k,o)*r;
            
            Ft_exc_5000(k,o)=Fo_exc_5000(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_5000(k,o)=Fo_exc_5000(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_5000(k,o)=-m2*r*omega5000^2;        
            Frc_exc_5000(k,o)=Fr_exc_5000(k,o)+Fro_exc_5000(k,o);
            Fvoc_exc_5000(k,o)=sqrt(Frc_exc_5000(k,o)^2+Ft_exc_5000(k,o)^2);
            Mtoc_exc_5000(k,o)=Ft_exc_5000(k,o)*r;
                                    
            Ft_exc_5250(k,o)=Fo_exc_5250(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_5250(k,o)=Fo_exc_5250(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_5250(k,o)=-m2*r*omega5250^2;        
            Frc_exc_5250(k,o)=Fr_exc_5250(k,o)+Fro_exc_5250(k,o);
            Fvoc_exc_5250(k,o)=sqrt(Frc_exc_5250(k,o)^2+Ft_exc_5250(k,o)^2);
            Mtoc_exc_5250(k,o)=Ft_exc_5250(k,o)*r;
                        
            Ft_exc_5500(k,o)=Fo_exc_5500(k,o).*(sin(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fr_exc_5500(k,o)=Fo_exc_5500(k,o).*(cos(alfa_exc(k)+beta_exc_rozsah(k,o)));
            Fro_exc_5500(k,o)=-m2*r*omega5500^2;        
            Frc_exc_5500(k,o)=Fr_exc_5500(k,o)+Fro_exc_5500(k,o);
            Fvoc_exc_5500(k,o)=sqrt(Frc_exc_5500(k,o)^2+Ft_exc_5500(k,o)^2);
            Mtoc_exc_5500(k,o)=Ft_exc_5500(k,o)*r;
            
    end
end

%% Grafy sil v ojnicích

figure(41)
hold on
plot(alfa_deg,Ft(:,18),'g--','LineWidth',2)
plot(alfa_deg, Ft_exc_5500(:,7),'g-','LineWidth',2)

plot(alfa_deg,Fr(:,18),'b--','LineWidth',2)
plot(alfa_deg,Fr_exc_5500(:,7),'b-','LineWidth',2)

plot(alfa_deg,Frc(:,18),'r--','LineWidth',2)
plot(alfa_deg,Frc_exc_5500(:,7),'r-','LineWidth',2)

plot(alfa_deg,Fvoc(:,18),'k--','LineWidth',2)
plot(alfa_deg,Fvoc_exc_5500(:,7),'k-','LineWidth',2)

title('Síly působící na ojniční čep, e=14mm, 5500min-1')
xlabel('Natočení klikového hřídele [°]')                               % popis osy x
ylabel('Síla působící na ojnici [N]')                               % popis osy y
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
grid on
legend('Tečná síla přenášená ojnicí CKM','Tečná síla přenášená ojnicí EKM','Radiální síla přenášená ojnicí CKM','Radiální síla přenášená ojnicí EKM','Výsledná radiální síla na ojničním čepu CKM','Výsledná radiální síla na ojničním čepu EKM','Výsledná zatěžující ojniční čep CKM','Výsledná zatěžující ojniční čep EKM')
hold off 

%% Výpočet točivého momentu na 1 zalomení

figure(42)
hold on
plot(alfa_deg,Mtoc_exc_5500(:,7),'r-','LineWidth',2)
plot(alfa_deg,Mkv1(:,18),'k--','LineWidth',2)
plot(alfa_deg,0,'k-','LineWidth',2)

title('Točivý moment na jednom zalomení klikové hřídele při 5500min-1')
xlabel('Natočení klikového hřídele [°]')                               % popis osy x
ylabel('Točivý moment na zalomení [Nm]')                               % popis osy y
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
grid on
legend('Točivý moment excentrického klikového mechanismu e=14mm','Točivý moment centrického klikového mechanismu')
hold off 


%% Výpočet točivého momentu celkově (tzn. u setrvačníku)

Mkv1_3750=circshift(Mtoc_exc_3750(:,7),[PV(1),0]);                   % posunutí průběhu točivého momentu 1. válce
Mkv2_3750=circshift(Mtoc_exc_3750(:,7),[PV(2),0]);                   % posunutí průběhu točivého momentu 2. válce
Mkv3_3750=circshift(Mtoc_exc_3750(:,7),[PV(3),0]);                   % posunutí průběhu točivého momentu 3. válce
Mkv4_3750=circshift(Mtoc_exc_3750(:,7),[PV(4),0]);                   % posunutí průběhu točivého momentu 4. válce

Mkv1_5500=circshift(Mtoc_exc_5500(:,7),[PV(1),0]);                   % posunutí průběhu točivého momentu 1. válce
Mkv2_5500=circshift(Mtoc_exc_5500(:,7),[PV(2),0]);                   % posunutí průběhu točivého momentu 2. válce
Mkv3_5500=circshift(Mtoc_exc_5500(:,7),[PV(3),0]);                   % posunutí průběhu točivého momentu 3. válce
Mkv4_5500=circshift(Mtoc_exc_5500(:,7),[PV(4),0]);                   % posunutí průběhu točivého momentu 4. válce

Mkc_3750=Mkv1_3750+Mkv2_3750+Mkv3_3750+Mkv4_3750;                    % celkový průběh točivého momentu
Mkc_5500=Mkv1_5500+Mkv2_5500+Mkv3_5500+Mkv4_5500;                    % celkový průběh točivého momentu


figure(43)
hold on
plot(alfa_deg,Mkc_3750,'r-','LineWidth',2)
plot(alfa_deg,Mkc(:,18),'k--','LineWidth',2)
plot(alfa_deg,0,'k-','LineWidth',2)

title('Výsledný točivý moment na setrvačníku při 5500min-1')
xlabel('Natočení klikového hřídele [°]')                               % popis osy x
ylabel('Točivý moment [Nm]')                               % popis osy y
xlim([0,720])                               % graf bude ukončen na x=720°
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
grid on
hold off 
legend('Točivý moment excentrického klikového mechanismu e=14mm','Točivý moment centrického klikového mechanismu')

 
strMk5550_14mm=mean(Mkc_5500)
strMk5550_0mm=mean(Mkc(:,18))

strabsMk5550_14mm=mean(abs(Mkc_5500))
strabsMk5550_0mm=mean(abs(Mkc(:,18)))


%% Točivý moment na setrvačníku při e=14mm, 5500min-1
figure(44)
hold on
plot(alfa_deg,Mkv1_5500,'k--','LineWidth',2)
plot(alfa_deg,Mkv2_5500,'b--','LineWidth',2)
plot(alfa_deg,Mkv3_5500,'m--','LineWidth',2)
plot(alfa_deg,Mkv4_5500,'g--','LineWidth',2)
plot(alfa_deg,Mkc_5500,'r-','LineWidth',2.5)
plot(alfa_deg,0,'k-','LineWidth',3)
plot(alfa_deg,mean(Mkc_5500),'r-','LineWidth',2.5)

title('Točivý moment na setrvačníku při e=14mm, 5500min-1')
xlabel('Natočení klikového hřídele [°]')                               % popis osy x
ylabel('Točivý moment [Nm]')                               % popis osy y
xlim([0,720])                               % graf bude ukončen na x=720°
ylim([-200,600])
set(gca,'XTick',0:180:720)                  % měřítko grafu jde po 180° skocích
grid on
hold off 
legend('Točivý moment od prvního válce','Točivý moment od druhého válce','Točivý moment od třetího válce','Točivý moment od čtvrtého válce','Výsledný točivý moment excentrického mechanismu')


%% Setrvačné síly posuvných částí I. řádu

for k=1:length(alfa)  
Fs_pIr_v1(k)=mpc*(omega5500^2)*r*(cos(alfa(k))+lambda_e_rozsah(7)*sin(alfa(k)));
Fs_pIr_v2(k)=mpc*(omega5500^2)*r*(cos(alfa(k)+3*pi)+lambda_e_rozsah(7)*sin(alfa(k)+3*pi));
Fs_pIr_v3(k)=mpc*(omega5500^2)*r*(cos(alfa(k)+pi)+lambda_e_rozsah(7)*sin(alfa(k)+pi));
Fs_pIr_v4(k)=mpc*(omega5500^2)*r*(cos(alfa(k)+2*pi)+lambda_e_rozsah(7)*sin(alfa(k)+2*pi));

Fs_pIr_vc(k)=Fs_pIr_v1(k)+Fs_pIr_v2(k)+Fs_pIr_v3(k)+Fs_pIr_v4(k);
end

figure(50)
hold on
plot(alfa_deg,Fs_pIr_v1,'k-','LineWidth',3)
plot(alfa_deg,Fs_pIr_v2,'b-','LineWidth',3)
plot(alfa_deg,Fs_pIr_v3,'m:','LineWidth',3)
plot(alfa_deg,Fs_pIr_v4,'g:','LineWidth',3)
plot(alfa_deg,Fs_pIr_vc,'r-','LineWidth',3)
set(gca,'XTick',0:90:360)
xlim([0,360]) 
grid on
title('Setrvačné síly posuvných částí I. řádu')
xlabel('Natočení klikového hřídele [°]')
ylabel('Setrvačná síla posuvných částí I. řádu [N]')
legend('Setrvačná síla posuvných částí I. řádu 1. válce','Setrvačná síla posuvných částí I. řádu 2. válce','Setrvačná síla posuvných částí I. řádu 3. válce','Setrvačná síla posuvných částí I. řádu 4. válce','Výslednice setrvačných sil posuvných částí I. řádu')
hold off

%% Setrvačné síly posuvných částí II. řádu
for k=1:length(alfa)  
Fs_pIIr_v1(k)=mpc*(omega5500^2)*r*(cos(2*(alfa(k))));
Fs_pIIr_v2(k)=mpc*(omega5500^2)*r*(cos(2*(alfa(k)+3*pi)));
Fs_pIIr_v3(k)=mpc*(omega5500^2)*r*(cos(2*(alfa(k)+pi)));
Fs_pIIr_v4(k)=mpc*(omega5500^2)*r*(cos(2*(alfa(k)+2*pi)));

Fs_pIIr_vc(k)=Fs_pIIr_v1(k)+Fs_pIIr_v2(k)+Fs_pIIr_v3(k)+Fs_pIIr_v4(k);
end

figure(51)
hold on
plot(alfa_deg,Fs_pIIr_v1,'k-','LineWidth',3)
plot(alfa_deg,Fs_pIIr_v2,'b-','LineWidth',3)
plot(alfa_deg,Fs_pIIr_v3,'m-.','LineWidth',3)
plot(alfa_deg,Fs_pIIr_v4,'g:','LineWidth',3)
plot(alfa_deg,Fs_pIIr_vc,'r-','LineWidth',3)
set(gca,'XTick',0:90:360)
xlim([0,360]) 
grid on
hold off
title('Setrvačné síly posuvných částí II. řádu')
xlabel('Natočení klikového hřídele [°]')
ylabel('Setrvačná síla posuvných částí II. řádu [N]')
legend('Setrvačná síla posuvných částí II. řádu 1. válce','Setrvačná síla posuvných částí II. řádu 2. válce','Setrvačná síla posuvných částí II. řádu 3. válce','Setrvačná síla posuvných částí II. řádu 4. válce','Výslednice setrvačných sil posuvných částí II. řádu')

figure(100)
hold on
plot(alfa_deg,Fs_pIIr_vc,'r--','LineWidth',3)
plot(alfa_deg,-0.5*Fs_pIIr_vc,'g--','LineWidth',3)
plot(alfa_deg,-0.5*Fs_pIIr_vc,'b:','LineWidth',3)
plot(alfa_deg,0*Fs_pIIr_vc,'k-','LineWidth',3)
set(gca,'XTick',0:90:360)
xlim([0,360]) 
grid on
hold off
title('Setrvačné síly posuvných částí II. řádu a jejich vyvážení')
xlabel('Natočení klikového hřídele [°]')
ylabel('Síla ve směru válce [N]')
legend('Výslednice setrvačných sil posuvných částí II. řádu','Vyvážovací síla 1. vyvažovací hřídele','Vyvážovací síla 2. vyvažovací hřídele','Výslednice')





%% ověření těch sin a cos :D
% 
% figure (300)
% hold on
% plot(alfa_deg,cosd(alfa_deg)+sind(alfa_deg),'k-','LineWidth',2)
% plot(alfa_deg,cosd(alfa_deg+540)+sind(alfa_deg+540),'b-','LineWidth',2)
% plot(alfa_deg,cosd(alfa_deg+180)+sind(alfa_deg+180),'m:','LineWidth',2)
% plot(alfa_deg,cosd(alfa_deg+360)+sind(alfa_deg+360),'g:','LineWidth',2)
% set(gca,'XTick',0:180:720)
% xlim([0,720]) 
% grid on

%% Setrvačné momenty posuvných částí I. řádu

% počítám k bodu A
for k=1:length(alfa)  
M_s_pIr_v1(k)=Fs_pIr_v1(k)*-(3/2)*a_vzdalenost_os_valcu;
M_s_pIr_v2(k)=Fs_pIr_v2(k)*+(1/2)*a_vzdalenost_os_valcu;
M_s_pIr_v3(k)=Fs_pIr_v3(k)*-(1/2)*a_vzdalenost_os_valcu;
M_s_pIr_v4(k)=Fs_pIr_v4(k)*+(3/2)*a_vzdalenost_os_valcu;

M_s_pIr_vc(k)=M_s_pIr_v1(k)+M_s_pIr_v2(k)+M_s_pIr_v3(k)+M_s_pIr_v4(k);
end

figure(52)
hold on
plot(alfa_deg,M_s_pIr_v1,'k--','LineWidth',2)
plot(alfa_deg,M_s_pIr_v2,'b--','LineWidth',2)
plot(alfa_deg,M_s_pIr_v3,'m--','LineWidth',2)
plot(alfa_deg,M_s_pIr_v4,'g--','LineWidth',2)
plot(alfa_deg,M_s_pIr_vc,'r-','LineWidth',3)
set(gca,'XTick',0:90:360)
title('Momenty setrvačných sil posuvných částí I. řádu')
xlabel('Natočení klikového hřídele [°]')
ylabel('Moment setrvačných sil posuvných částí I. řádu [Nm]')
legend('Moment setrvačných sil posuvných částí I. řádu 1. válce','Moment setrvačných sil posuvných částí I. řádu 2. válce','Moment setrvačných sil posuvných částí I. řádu 3. válce','Moment setrvačných sil posuvných částí I. řádu 4. válce','Výslednice momentu setrvačných sil posuvných částí I. řádu')


xlim([0,360]) 
grid on
hold off


%% Setrvačné momenty posuvných částí II. řádu

% počítám k bodu T
for k=1:length(alfa)  
M_s_pIIr_v1(k)=-Fs_pIIr_v1(k)*(3/2)*a_vzdalenost_os_valcu;
M_s_pIIr_v2(k)=-Fs_pIIr_v2(k)*(1/2)*a_vzdalenost_os_valcu;
M_s_pIIr_v3(k)=Fs_pIIr_v3(k)*(1/2)*a_vzdalenost_os_valcu;
M_s_pIIr_v4(k)=Fs_pIIr_v4(k)*(3/2)*a_vzdalenost_os_valcu;

M_s_pIIr_vc(k)=M_s_pIIr_v1(k)+M_s_pIIr_v2(k)+M_s_pIIr_v3(k)+M_s_pIIr_v4(k);
end

figure(53)
hold on
plot(alfa_deg,M_s_pIIr_v1,'k--','LineWidth',2)
plot(alfa_deg,M_s_pIIr_v2,'b--','LineWidth',2)
plot(alfa_deg,M_s_pIIr_v3,'m--','LineWidth',2)
plot(alfa_deg,M_s_pIIr_v4,'g--','LineWidth',2)
plot(alfa_deg,M_s_pIIr_vc*0,'r-','LineWidth',3)
set(gca,'XTick',0:90:360)
xlim([0,360]) 
grid on

title('Momenty setrvačných sil posuvných částí II. řádu')
xlabel('Natočení klikového hřídele [°]')
ylabel('Moment setrvačných sil posuvných částí II. řádu [Nm]')
legend('Moment setrvačných sil posuvných částí II. řádu 1. válce','Moment setrvačných sil posuvných částí II. řádu 2. válce','Moment setrvačných sil posuvných částí II. řádu 3. válce','Moment setrvačných sil posuvných částí II. řádu 4. válce','Výslednice momentu setrvačných sil posuvných částí II. řádu')

hold off

%% "Přídavný" klopný příčný moment od výslednice setrvačných sil posuvných částí II. řádu

for k=1:length(alfa) 

    for o=1:1:15
M_klopny_pridavny(k,o)=Fs_pIIr_vc(k)*(o*0.002);
    end


end

figure(54)
hold on
plot(alfa_deg,M_klopny_pridavny(:,5)*0,'k--','LineWidth',2)
plot(alfa_deg,M_klopny_pridavny(:,5),'g--','LineWidth',2)
plot(alfa_deg,M_klopny_pridavny(:,7),'b-','LineWidth',3)
plot(alfa_deg,M_klopny_pridavny(:,10),'m-','LineWidth',2)
plot(alfa_deg,M_klopny_pridavny(:,15),'r--','LineWidth',2)
hold off
title('Přídavný příčný klopný moment pro 5500min-1')
xlabel('Natočení klikového hřídele [°]')
ylabel('Přídavný příčný klopný moment [Nm]')
legend('e=0mm','e=10mm','e=14mm','e=20mm','e=30mm')
set(gca,'XTick',0:90:360)
xlim([0,360]) 
grid on

