%Vypocet zatizitelnosti tramu v POLI 

%MaxMy cell array, kde kazdy cell je pro jedno zatizeni (=excel).
%Kazdy cell ma 2D array
%       Radek je jeden tram. Je 1 az 5 tramu na peti radkach.
%       Kazdy sloupec je Moment v jednom prurezu.

% resime normalni zatizitelnost, pouze silnicni, sestava n1
%v01 - funkcni - provedena rucni kontrola rucni_kontrola.m pro 
%       Overeni vypoctu zatizitelnosti pro tram 5, prurez B1
%       pouze silnicni, tj. w=6 m, kombinace 6.10a, sestava zatizeni "n1"
%       moment únosnosti toho prurezu spocitan v excelu "stanoveni momentu
%       unosnosti"
%v02 - rozsireni na sestava zatizeni n2, r1, r2, e1
%v03 - rozsireni na zeleznici - B1 
%v04 - zmena vysek tramu v jednotlivych polich
%v05 - zmena v rovnicich, chodniky dany do citatele
%v06 - stale zatizeni z TDA, + M vodorovne z B1 u varianty C
%v06_redukceteploty_01 - v rovnicich pro stanoveni zatizitelnosti redukovany posuny a teplota na 75%
%v06_redukceteploty_02 - k linearni teplote na mostovku doplnen soucinitel ksur=0.5 - zohledneni tloustky vrstvy mostniho svrsku
%                       - upravena chyba u vyhradni zatizitelnosti - vybrano vozidlo pouze na leve strane mostu 
%v08_redukceteploty_01 - doplneno pouze zeleznicni
%v09 - + pro ruzne vatianty teploty a vetru
%v10 - + pro zeleznicni v zavislosti na silnicni
%    - pridana zaleznicni prechodnost
%v11 - pridana kombinace 6.10b

clc;clear all
%vyber kombinace
% TypKombinace=['6_10_a'];%
TypKombinace=['6_10_b'];%

%zakomentuj typ zatizeni
TypZatizeni=['RedukceTeploty'];%
% TypZatizeni=['BezTeploty'];%
% TypZatizeni=['BezVetru'];%
% TypZatizeni=['Bez_vetru_i_teploty'];%

ulozExcel=0; %maji se ulozit tabulky do excelu?
ulozExcel=logical(ulozExcel);

load('TramyPoleMaxM_Zat.mat')
load('TabulkyTDA.mat');%NACTENI DAT z TDA

Fyd=180;%MPa mez kluzu navrhova

gama_g=1.35;%soucinitel spolehlivosti asi - zatim jen pro silnici
gama_g_pos=1.2;%soucinitel spolehlivosti
gama_q=1.5;%soucinitel spolehlivosti
gama_v=1.35;%soucinitel spolehlivosti - pouze pro silnicni

psi_0i=0.6;%soucinitel kombinace pro promenne vedlejsi zatizeni (vitr, teplota)
psi_01=0.75;%soucinitel kombinace pro zatizitelnost
psi_11=0.75;%soucinitel kombinace pro sestavu n2

delta=1.05;%dynamicky soucinitel silnicni pro vyjimecnou zatizitelnost
delta_1=1.25;%dynamicky soucinitel silnicni pro vyhradni zatizitelnost
delta_2=1.20;%dynamicky soucinitel silnicni pro normalni zatizitelnost

%pro zeleznici
gama_zel=1.45;%dilci soucinitel pro zeleznicni dopravu B1
phiT1=[1.23*ones(1,4),1.32*ones(1,10),1.24*ones(1,4)];% DYNAmicky soucinitel pro standardne udrzovanou kolej - pro B1

switch TypKombinace
    case '6_10_a'
        ksi=1;
        ksi_zel=1;
        psi_01=0.75;%soucinitel kombinace pro zatizitelnost
        psi_01_zel=0.8;%dilci soucinitel pro zeleznicni dopravu
    case '6_10_b'
        ksi=0.85;
        ksi_zel=0.95;
        psi_01=1;%soucinitel kombinace pro zatizitelnost
        psi_01_zel=1;%dilci soucinitel pro zeleznicni dopravu
    otherwise
        error('Vybran spatny typ zatizeni')
end

MaxMy=data.MaxMy;
MaxMy=cellfun(@abs,MaxMy,'UniformOutput',false);%absolutni hodnota na vsechny bunky cellu
MmaxGm=DATA2{1}.D;          %Stale zatizeni prepocitane dle TDA
MmaxDeltaSup=DATA2{6}.D;    %Prirustech ostatniho staleho zatizeni k_sup
for i1 = 1:5 %postupne po jednotlivych tramech
    if or(i1==1,i1==5)
        Prurezy=[1:8,11:18];%ma jit az do 18 - napsat [1:8,11:18]
    else
        Prurezy=[1:18];%ma jit az do 18 - napsat [1:18]
    end
    for i2=Prurezy% do 18 postupne po jednotlivych prurezech
        if i1==5 & i2==1
            disp('Ted')
        end
        %% --------------stale a vedlejsi promenne zatizeni------------------------------------------------------------------- 
        M_RD_N(i1,i2)=MRD_tramyPole_02(i1,data.Pole{i1,i2},data.DelkaPole(i2),data.iPrutu{i1,i2},Fyd,data.H(i1,i2));%Stejna struktura jako As
%         Mg=MaxMy{1}(i1,i2); ze 3D
%         Most=MaxMy{3}(i1,i2); ze 3D
        Mg=MmaxGm(i1,i2)+MmaxDeltaSup(i1,i2); %z TDA      
        Mpos=0.75*(MaxMy{5}(i1,i2)+MaxMy{6}(i1,i2));%Posuny vodorovne a svisle  
        koefVitr=max([MaxMy{7}(i1,i2)+MaxMy{8}(i1,i2)+MaxMy{9}(i1,i2),MaxMy{10}(i1,i2)+MaxMy{11}(i1,i2)+MaxMy{12}(i1,i2)]);%koeficient vitr
        
        
        %Prevod na jednu teplotu 
        k_sur=0.5;% soucinitel zohlednujici tlousku mostniho svrsku
        DeltaT_m=k_sur*MaxMy{14}(i1,i2)+MaxMy{15}(i1,i2)+MaxMy{16}(i1,i2)+MaxMy{17}(i1,i2);
        om_m=0.75;
        DeltaT_n=MaxMy{13}(i1,i2);
        om_n=0.35;
        koefTeploty=0.75*(max([DeltaT_m+om_n*DeltaT_n,om_m*DeltaT_m+DeltaT_n]));
        switch TypZatizeni
            case 'RedukceTeploty'
                filename=['ZatTramyPole_redukceTeploty'];
                disp('RedukceTeploty')
            case 'BezTeploty'
                filename=['ZatTramyPole_bezTeploty'];
                koefTeploty=0;
                disp('BezTeploty')
            case 'BezVetru'
                filename=['ZatTramyPole_bezVetru'];
                koefVitr=0;
                disp('BezVetru')
            case 'Bez_vetru_i_teploty'
                filename=['ZatTramyPole_bez_vetru_i_teploty'];
                koefTeploty=0;
                koefVitr=0;
                disp('Bez_vetru_i_teploty')
            otherwise
                error('Spatne zadany typ zatizeni')
        end
        
        M_chod=MaxMy{26}(i1,i2)+MaxMy{27}(i1,i2);%zatizeni chodniku
        
        %% --------------moment od dopravy, chodniku a vodorovnych sil - normalni zatizitelnost-------------------------------
        delta_i=delta_2;%dynamicky soucinitel - zatim jen silnicni
        M_vodr_n=MaxMy{32}(i1,i2);
        M_vk_c1_n=max([MaxMy{18}(i1,i2),MaxMy{19}(i1,i2),MaxMy{20}(i1,i2)])+max([MaxMy{22}(i1,i2),MaxMy{23}(i1,i2),MaxMy{24}(i1,i2)]);
        M_vk_ci_n=M_vk_c1_n;%moment od dopravy
        
        %sestava zatizeni n1 - normalni zatizitelnost (pouze silnicni)
        v_c1a_n1(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty)-(gama_v*psi_01*M_chod))/(gama_v*psi_01*delta_i*M_vk_ci_n);
        %sestava zatizeni n2 - normalni zatizitelnost (pouze silnicni)
        v_c1a_n2(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty))/(gama_v*psi_01*(delta_i*psi_11*M_vk_ci_n+M_vodr_n));

        %% ---------------moment od dopravy, chodniku a vodorovnych sil - vyhradni zatizitelnost--------------------------------- 
        delta_i=delta_1;%dynamicku soucinitel - zatim jen silnicni
        M_vodr_r=MaxMy{34}(i1,i2);
        M_vk_c2_r=max([MaxMy{28}(i1,i2),MaxMy{29}(i1,i2)]);
        M_vk_ci_r=M_vk_c2_r;%moment od dopravy
         
       
        %sestava zatizeni r1 - vyhradni zatizitelnost (pouze silnicni)
        Vrw_c2_r1(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty)-(gama_v*psi_01*M_chod))/(gama_v*psi_01*delta_i*M_vk_ci_r);
        %sestava zatizeni r2 - vyhradni zatizitelnost (pouze silnicni)
        Vrw_c2_r2(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty))/(gama_v*psi_01*(delta_i*psi_11*M_vk_ci_r+M_vodr_r));
        
        %% ---------------moment od dopravy, chodniku a vodorovnych sil - vyjimecna zatizitelnost---------------------------------      
        delta_i=delta;%dynamicky soucinitel - zatim jen silnicni
        M_vk_c3_e=max([MaxMy{30}(i1,i2),MaxMy{31}(i1,i2)]);
        M_vk_ci_e=M_vk_c3_e;%moment od dopravy
        %sestava zatizeni e1 - vyjimecna zatizitelnost (pouze silnicni)
        Vew_c3_e1(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty))/(gama_v*psi_01*(delta_i*M_vk_ci_e));
        
%% --------------moment od dopravy + B1, chodniku a vodorovnych sil - normalni zatizitelnost-------------------------------
        delta_i=delta_1;%dynamicky soucinitel - zatim silnicni + B1 - je jen jeden pruh silnicni
        
        psi_11_zel=0.8;
        M_vodr_n=MaxMy{33}(i1,i2);
        M_vk_c1_n=MaxMy{21}(i1,i2)+MaxMy{25}(i1,i2);%- silnicni v jednom pruhu 
        M_vk_ci_n=M_vk_c1_n;%moment od dopravy silnicni
        M_B1=MaxMy{42}(i1,i2);%moment od dopravy B1
        M_B1_vodr=MaxMy{44}(i1,i2);%momenst od vodorovnych sil B1
        
        
        %sestava zatizeni n1 - normalni zatizitelnost ( silnicni + B1)
        v_c1a_n1_B1(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty)-(gama_zel*psi_01_zel*(phiT1(i2)*M_B1))-(gama_v*psi_01*M_chod))/(gama_v*psi_01*delta_i*M_vk_ci_n);
        %sestava zatizeni n2 - normalni zatizitelnost ( silnicni + B1)
        v_c1a_n2_B1(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty)-gama_zel*psi_01_zel*(psi_11_zel*phiT1(i2)*M_B1+M_B1_vodr))/(gama_v*psi_01*(delta_i*psi_11*M_vk_ci_n+M_vodr_n));

        %% ---------------moment od dopravy + B1, chodniku a vodorovnych sil - vyhradni zatizitelnost--------------------------------- 
        delta_i=delta_1;%dynamicku soucinitel - zatim  silnicni + B1
        M_vodr_r=MaxMy{34}(i1,i2);
        M_vk_c2_r=MaxMy{28}(i1,i2);%opraveno, vybrana pouze doprava v leve strane
        M_vk_ci_r=M_vk_c2_r;%moment od dopravy
        M_B1=MaxMy{42}(i1,i2);%moment od dopravy B1 
        M_B1_vodr=MaxMy{44}(i1,i2);%momenst od vodorovnych sil B1
        
        %sestava zatizeni r1 - vyhradni zatizitelnost (silnicni + B1)
        Vrw_c2_r1_B1(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty)-(gama_zel*psi_01_zel*(phiT1(i2)*M_B1))-(gama_v*psi_01*M_chod))/(gama_v*psi_01*delta_i*M_vk_ci_r);
        %sestava zatizeni r2 - vyhradni zatizitelnost (silnicni + B1)
        Vrw_c2_r2_B1(i1,i2)=(M_RD_N(i1,i2)-(gama_g*ksi*Mg+gama_g_pos*ksi*Mpos)-(gama_q*psi_0i*koefVitr+gama_q*psi_0i*koefTeploty)-gama_zel*psi_01_zel*(psi_11_zel*phiT1(i2)*M_B1+M_B1_vodr))/(gama_v*psi_01*(delta_i*psi_11*M_vk_ci_r+M_vodr_r));
         
%% ------------- zeleznicni zatizitelnost - normalni 
        M_nv=max([MaxMy{35}(i1,i2),MaxMy{36}(i1,i2),MaxMy{37}(i1,i2)]);     %moment od kolovych sil
        M_UDL=max([MaxMy{38}(i1,i2),MaxMy{39}(i1,i2),MaxMy{40}(i1,i2)]);    %moment od rovnomerneho spojiteho zatizeni
        M_vodr_zel=MaxMy{41}(i1,i2);        %vodorovne sily od zeleznice
        M_1_chod=MaxMy{26}(i1,i2);          %zatizeni jednoho chodniku na strane kde nejezdi vlak
        
        [Z_LM71(i1,i2), koef_Zel_zat] = Zel_zat_01(i2,M_RD_N(i1,i2),ksi_zel*Mg,ksi_zel*Mpos,koefVitr,koefTeploty,M_1_chod,M_nv,M_UDL,M_vodr_zel,TypZatizeni,psi_01_zel);%zeleznicni zatizitelnost - normalni 
%%      zeleznizni prechodnost

        psi=phiT1(i2)/koef_Zel_zat.phi3(i2); %součinitel dynamické redukce
        %ucinnost provozniho zatizeni
        M_LM71=koef_Zel_zat.gama_zel*(0.488*M_nv+M_UDL+M_vodr_zel);%navrhovy moment od LM71
        M_B1_zel=koef_Zel_zat.gama_zel*(M_B1+M_B1_vodr); %navrhovy moment od provozniho zatizeni B1
        LambdaLM71=M_B1_zel/M_LM71;         %ucinnost provozniho zatizeni
        
        Pr_zel(i1,i2)=psi*LambdaLM71;
        
        if Z_LM71(i1,i2)>=psi*LambdaLM71
            prechodnost_zel{i1,i2}='Ok';
        else
            prechodnost_zel{i1,i2}='Ne';
        end
%% ------------- zeleznicni zatizitelnost + silnice - normalni 
        M_vodr_n=MaxMy{33}(i1,i2);
        M_vk_c1_n=MaxMy{21}(i1,i2)+MaxMy{25}(i1,i2);%- silnicni v jednom pruhu    
        [Z_LM71_ZS(i1,i2), koef_ZS_zat] = Zel_zat_v_zavislosi_na_silnicni_01(i2,M_RD_N(i1,i2),ksi_zel*Mg,ksi_zel*Mpos,koefVitr,koefTeploty,M_1_chod,M_nv,M_UDL,M_vodr_zel,M_vodr_n,M_vk_c1_n,TypZatizeni,psi_01_zel);
%%      zeleznizni prechodnost v zavislosti na silnicni

        psi_ZS=phiT1(i2)/koef_ZS_zat.phi3(i2); %součinitel dynamické redukce
        %ucinnost provozniho zatizeni
        M_LM71_ZS=koef_ZS_zat.gama_zel*(0.488*M_nv+M_UDL+M_vodr_zel);%navrhovy moment od LM71
        M_B1_ZS=koef_ZS_zat.gama_zel*(M_B1+M_B1_vodr); %navrhovy moment od provozniho zatizeni B1
        LambdaLM71_ZS=M_B1_ZS/M_LM71_ZS;         %ucinnost provozniho zatizeni
        
        Pr_ZS(i1,i2)=psi_ZS*LambdaLM71_ZS;
        
        if Z_LM71_ZS(i1,i2)>=psi_ZS*LambdaLM71_ZS
            prechodnost_ZS{i1,i2}='Ok';
        else
            prechodnost_ZS{i1,i2}='Ne';
        end
   
    end
end

%% Normalni zatizitelnost
fprintf('\n ------------ Normalni zatizitelnost ----------\n')
%stanovení sily na zadni napravu
V_c1a_n1=100*v_c1a_n1;%kN
V_c1a_n2=100*v_c1a_n2;%kN

%tíha vozidla
Vnw_c1_n1=(4/3)*V_c1a_n1;%kN
Vnw_c1_n2=(4/3)*V_c1a_n2;%kN

%normalni zatizitelnost
Vn_c1_n1=(1/10)*Vnw_c1_n1%t
Vn_c1_n2=(1/10)*Vnw_c1_n2%t

%najdi min
Vn_c1_n1_pom=Vn_c1_n1;
Vn_c1_n1_pom(Vn_c1_n1==0)= NaN;%nahrazeni nul NaN
Vn_c1_n2_pom=Vn_c1_n2;
Vn_c1_n2_pom(Vn_c1_n2==0)= NaN;%nahrazeni nul NaN

Vn_c1_n1_min=min((min(Vn_c1_n1_pom,[],'omitnan')))%[t] normalni zatizitelnost
Vn_c1_n2_min=min((min(Vn_c1_n2_pom,[],'omitnan')))%[t] normalni zatizitelnost
%% Vyhradni zatizitelnost
fprintf('\n ------------ Vyhradni zatizitelnost ----------\n')
Vr_c2_r1=(1/10)*Vrw_c2_r1%t
Vr_c2_r2=(1/10)*Vrw_c2_r2%t
%najdi min
Vr_c2_r1_pom=Vr_c2_r1;
Vr_c2_r1_pom(Vr_c2_r1==0)= NaN;%nahrazeni nul NaN
Vr_c2_r2_pom=Vr_c2_r2;
Vr_c2_r2_pom(Vr_c2_r2==0)= NaN;%nahrazeni nul NaN

Vr_c2_r1_min=min((min(Vr_c2_r1_pom,[],'omitnan')))%[t] vyhradni zatizitelnost
Vr_c2_r2_min=min((min(Vr_c2_r2_pom,[],'omitnan')))%[t] vyhradni zatizitelnost

%% Vyjimecna zatizitelnost
fprintf('\n ------------ Vyjimecna zatizitelnost ----------\n')
Ve_c3_e1=(1/10)*Vew_c3_e1%t

%najdi min
Ve_c3_e1_pom=Ve_c3_e1;
Ve_c3_e1_pom(Ve_c3_e1==0)= NaN;%nahrazeni nul NaN


Ve_c3_e1_min=min((min(Ve_c3_e1_pom,[],'omitnan')))%[t] vyjimecna zatizitelnost
%% Normalni zatizitelnost po odecteni B1
fprintf('\n ------------ Normalni zatizitelnost po odecteni B1----------\n')
%stanovení sily na zadni napravu
V_c1a_n1_B1=100*v_c1a_n1_B1;%kN
V_c1a_n2_B1=100*v_c1a_n2_B1;%kN

%tíha vozidla
Vnw_c1_n1_B1=(4/3)*V_c1a_n1_B1;%kN
Vnw_c1_n2_B1=(4/3)*V_c1a_n2_B1;%kN

%normalni zatizitelnost
Vn_c1_n1_B1=(1/10)*Vnw_c1_n1_B1%t
Vn_c1_n2_B1=(1/10)*Vnw_c1_n2_B1%t

%najdi min
Vn_c1_n1_pom_B1=Vn_c1_n1_B1;
Vn_c1_n1_pom_B1(Vn_c1_n1_B1==0)= NaN;%nahrazeni nul NaN
Vn_c1_n2_pom_B1=Vn_c1_n2_B1;
Vn_c1_n2_pom_B1(Vn_c1_n2_B1==0)= NaN;%nahrazeni nul NaN

Vn_c1_n1_B1_min=min((min(Vn_c1_n1_pom_B1,[],'omitnan')))%[t] normalni zatizitelnost po odecteni B1
Vn_c1_n2_B1_min=min((min(Vn_c1_n2_pom_B1,[],'omitnan')))%[t] normalni zatizitelnost po odecteni B1
%% Vyhradni zatizitelnost po odecteni B1
fprintf('\n ------------ Vyhradni zatizitelnost po odecteni B1----------\n')
Vr_c2_r1_B1=(1/10)*Vrw_c2_r1_B1%t
Vr_c2_r2_B1=(1/10)*Vrw_c2_r2_B1%t
%najdi min
Vr_c2_r1_pom_B1=Vr_c2_r1_B1;
Vr_c2_r1_pom_B1(Vr_c2_r1_B1==0)= NaN;%nahrazeni nul NaN
Vr_c2_r2_pom_B1=Vr_c2_r2_B1;
Vr_c2_r2_pom_B1(Vr_c2_r2_B1==0)= NaN;%nahrazeni nul NaN

Vr_c2_r1_B1_min=min((min(Vr_c2_r1_pom_B1,[],'omitnan')))%[t] vyhradni zatizitelnost po odecteni B1
Vr_c2_r2_B1_min=min((min(Vr_c2_r2_pom_B1,[],'omitnan')))%[t] vyhradni zatizitelnost po odecteni B1
%% shrnuti vysledku do jedne promenne
zatizitelnostMSU.TramyPole.a{1}{1}=['Normalni zatizitelnost - sestava n1'];
zatizitelnostMSU.TramyPole.a{1}{2}=Vn_c1_n1_min;
zatizitelnostMSU.TramyPole.a{1}{3}=Vn_c1_n1;

zatizitelnostMSU.TramyPole.a{2}{1}=['Normalni zatizitelnost - sestava n2'];
zatizitelnostMSU.TramyPole.a{2}{2}=Vn_c1_n2_min;
zatizitelnostMSU.TramyPole.a{2}{3}=Vn_c1_n2;

zatizitelnostMSU.TramyPole.a{3}{1}=['Vyhradni zatizitelnost - sestava r1'];
zatizitelnostMSU.TramyPole.a{3}{2}=Vr_c2_r1_min;
zatizitelnostMSU.TramyPole.a{3}{3}=Vr_c2_r1;

zatizitelnostMSU.TramyPole.a{4}{1}=['Vyhradni zatizitelnost - sestava r2'];
zatizitelnostMSU.TramyPole.a{4}{2}=Vr_c2_r2_min;
zatizitelnostMSU.TramyPole.a{4}{3}=Vr_c2_r2;

zatizitelnostMSU.TramyPole.a{5}{1}=['Vyjimecna zatizitelnost'];
zatizitelnostMSU.TramyPole.a{5}{2}=Ve_c3_e1_min;
zatizitelnostMSU.TramyPole.a{5}{3}=Ve_c3_e1;

zatizitelnostMSU.TramyPole.b{1}{1}=['Normalni zatizitelnost zeleznicni - sestava gr11'];
zatizitelnostMSU.TramyPole.b{1}{2}='Nevyplneno, zatim nepotrebne';
zatizitelnostMSU.TramyPole.b{1}{3}=Z_LM71;

zatizitelnostMSU.TramyPole.b{2}{1}=['Zeleznicni prechodnost'];
zatizitelnostMSU.TramyPole.b{2}{2}='Nevyplneno, zatim nepotrebne';
zatizitelnostMSU.TramyPole.b{2}{3}=Pr_zel;
zatizitelnostMSU.TramyPole.b{2}{4}=prechodnost_zel;

zatizitelnostMSU.TramyPole.c{1}{1}=['Normalni zatizitelnost po odecteni B1 - sestava n1'];
zatizitelnostMSU.TramyPole.c{1}{2}=Vn_c1_n1_B1_min;
zatizitelnostMSU.TramyPole.c{1}{3}=Vn_c1_n1_B1;

zatizitelnostMSU.TramyPole.c{2}{1}=['Normalni zatizitelnost po odecteni B1 - sestava n2'];
zatizitelnostMSU.TramyPole.c{2}{2}=Vn_c1_n2_B1_min;
zatizitelnostMSU.TramyPole.c{2}{3}=Vn_c1_n2_B1;

zatizitelnostMSU.TramyPole.c{3}{1}=['Vyhradni zatizitelnost po odecteni B1 - sestava r1'];
zatizitelnostMSU.TramyPole.c{3}{2}=Vr_c2_r1_B1_min;
zatizitelnostMSU.TramyPole.c{3}{3}=Vr_c2_r1_B1;

zatizitelnostMSU.TramyPole.c{4}{1}=['Vyhradni zatizitelnost po odecteni B1 - sestava r2'];
zatizitelnostMSU.TramyPole.c{4}{2}=Vr_c2_r2_B1_min;
zatizitelnostMSU.TramyPole.c{4}{3}=Vr_c2_r2_B1;

zatizitelnostMSU.TramyPole.d{1}{1}=['Normalni zatizitelnost zeleznicni a silnicni - sestava gr11'];
zatizitelnostMSU.TramyPole.d{1}{2}='Nevyplneno, zatim nepotrebne';
zatizitelnostMSU.TramyPole.d{1}{3}=Z_LM71_ZS;

zatizitelnostMSU.TramyPole.d{2}{1}=['Zeleznicni prechodnost v zavislosti na silnicni'];
zatizitelnostMSU.TramyPole.d{2}{2}='Nevyplneno, zatim nepotrebne';
zatizitelnostMSU.TramyPole.d{2}{3}=Pr_ZS;
zatizitelnostMSU.TramyPole.d{2}{4}=prechodnost_ZS;

zatizitelnostMSU.TramyPole.Zatizeni=TypZatizeni;
zatizitelnostMSU.TramyPole.TypKombinace=TypKombinace;
%% export do ecxelu
if ulozExcel
    warning('Pred ukladanim je potreba smazat puvodni excel!  Jinak hozi ze zustanoustare hodnoty na mistech, kde se nezapisuje.')
    
    w=cd
    switch TypKombinace
        case '6_10_a'
            cd([w,'/vysledky/6_10_a']);
        case '6_10_b'
            cd([w,'/vysledky/6_10_b']);
    end
    %vyber zadanych radku a sloupecku do excelu
    sL1=[5:8];  %sloupecky porvni vyber
    sL2=[11:14];%sloupecky druhy vyber
    
    Vn_c1_n1=[Vn_c1_n1(:,sL1),Vn_c1_n1(:,sL2)];   
    Vn_c1_n2=[Vn_c1_n2(:,sL1),Vn_c1_n2(:,sL2)];
    Vr_c2_r1=[Vr_c2_r1(:,sL1),Vr_c2_r1(:,sL2)];
    Vr_c2_r2=[Vr_c2_r2(:,sL1),Vr_c2_r2(:,sL2)]; 
    Ve_c3_e1=[Ve_c3_e1(:,sL1),Ve_c3_e1(:,sL2)];
    Vn_c1_n1_B1=[Vn_c1_n1_B1(:,sL1),Vn_c1_n1_B1(:,sL2)];
    Vn_c1_n2_B1=[Vn_c1_n2_B1(:,sL1),Vn_c1_n2_B1(:,sL2)];
    Vr_c2_r1_B1=[Vr_c2_r1_B1(:,sL1),Vr_c2_r1_B1(:,sL2)];
    Vr_c2_r2_B1=[Vr_c2_r2_B1(:,sL1),Vr_c2_r2_B1(:,sL2)];  
    Z_LM71=[Z_LM71(:,sL1),Z_LM71(:,sL2)];
    Pr_zel=[Pr_zel(:,sL1),Pr_zel(:,sL2)];
    Z_LM71_ZS=[Z_LM71_ZS(:,sL1),Z_LM71_ZS(:,sL2)];
    Pr_ZS=[Pr_ZS(:,sL1),Pr_ZS(:,sL2)];

    %převede matice vysledku do listů v excelu
    xlswrite(filename, Vn_c1_n1, 'Vn_c1_n1');
    xlswrite(filename, Vn_c1_n2, 'Vn_c1_n2');
    xlswrite(filename, Vr_c2_r1, 'Vr_c2_r1');
    xlswrite(filename, Vr_c2_r2, 'Vr_c2_r2');
    xlswrite(filename, Ve_c3_e1, 'Ve_c3_e1');
    
    xlswrite(filename, Vn_c1_n1_B1, 'Vn_c1_n1_B1');
    xlswrite(filename, Vn_c1_n2_B1, 'Vn_c1_n2_B1');
    xlswrite(filename, Vr_c2_r1_B1, 'Vr_c2_r1_B1');
    xlswrite(filename, Vr_c2_r2_B1, 'Vr_c2_r2_B1');
      
    xlswrite(filename, Z_LM71, 'Z_LM71');
    xlswrite(filename, Pr_zel, 'Pr_zel');
    xlswrite(filename, Z_LM71_ZS, 'Z_LM71_ZS');
    xlswrite(filename, Pr_ZS, 'Pr_ZS');
    cd([w]);
end