%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Parameters of nonlinear mathematical model of HPS.                    %
%                                                                         %
%                                                                         %
%                                                                         %
%                                        Datum posledni zmeny: 20.06.2010 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Konstanty soustavy
R=287;          %[J/K.kg]   Merna plynova konstanta vzduchu
T=298.15;       %[K]        Teplota
rok=1000;       %[kg/m3]    Hustota vody
g=9.81;         %[m/s2]     Tihove zrychleni 
cp=1005;        %[J/(K*kg)] Merna tepelna kapacita vzduchu 
m=1.35;         %[]         Poissonova konstanta pro vzduch

%--------------------------------------------------------------------------
%Geometricke rozmery soustavy

%Plocha prurezu nadrzi
rL=0.02;        %[m]        Polomer levych nadrzi 
rR=0.025;       %[m]        Polomer pravych nadrzi 
SL=pi*rL^2;     %[m2]       Plocha prurezu levych nadrzi 
SR=pi*rR^2;     %[m2]       Plocha prurezu pravych nadrzi 
clear rL rR

%Pneumaticky objem
H0=0.35;        %[m]        Celkova vyska hydraulickych nadrzi
gamaH=0.85;   %[]         Korekce horniho volneho pneumat. objemu 
gamaL=0.74;   %[]         Korekce dolniho volneho pneumat. objemu

r=0.025;        %[m]        Polomer dopl. pneumatickeho objemu
H=0.17;         %[m]        Vyska dopl. pneumatickeho objemu 
V0=pi*r^2*H;    %[m3]       Objem dopl. pneumatickeho objemu 
clear r H

%Vytokove clony
pk=2.9e-3;      %[m]        Prurez vytokovych otvoru pro kapalinu 
sLH=pi*pk^2/4;  %[m2]       Plocha vytokove clony leve horni nadrze
sLL=sLH;        %[m2]       Plocha vytokove clony leve dolni nadrze
sRH=sLH;        %[m2]       Plocha vytokove clony prave horni nadrze
sRL=sLH;        %[m2]       Plocha vytokove clony prave dolni nadrze
alfaLH=0.6940;  %[]         Korekcni clen pro levou horni vytok. clonu
alfaRH=0.7365;  %[]         Korekcni clen pro pravou horni vytok. clonu
alfaLL=0.6680;  %[]         Korekcni clen pro levou dolni vytok. clonu
alfaRL=0.7339;  %[]         Korekcni clen pro pravou dolni vytok. clon
clear pk

%Plocha odvzdusnovaci clony
s0H=1.96e-7;    %[m2]       Plocha clony horniho pneumat. objemu 
s0L=1.96e-7;    %[m2]       Plocha clony dolniho pneumat. objemu
betaH=0.2367;   %[]         Korekce plochy clony horniho pneumat. objemu 
betaL=0.6748;   %[]         Korekce plochy clony dolniho pneumat. objemu
KH=0.0161;
KL=0.1535;

%VENTILY
%Levy ventil
aLV=0.6300;
auL=1.3342;
sLV=5e-6;       %[m2]       Plocha sterbiny v levem ventilu
suL=5e-008;     %[m2]       Podtok leveho ventilu
uLmin=0.2919;     %[]         Rozsah ovladatelnosti leveho ventilu
uLmax=0.7828;     %[] 

%Pravy ventil
aRV=0.6061;
auR=2.3974;
sRV=5e-6;       %[m2]       Plocha sterbiny v praveho ventilu
suR=suL;        %[m2]       Podtok praveho ventilu
uRmin=0.2234;     %[]         Rozsah ovladatelnosti praveho ventilu
uRmax=0.7765;     %[]

%Vyskova cidla
%Leve horni cidlo
kLH=0.031860;   % [m/V] smernice
ofsLH=0.043;    % [m] ofset

%Prave horni cidlo
kRH=0.030421;       % [m/V] smernice
ofsRH=0.043;    % [m] ofset

%Leve dolni cidlo
kLL=0.032269;       % [m/V] smernice
ofsLL=0.0385;    % [m] ofset

%Prave dolni cidlo
kRL=0.030889;       % [m/V] smernice
ofsRL=0.0385;    % [m] ofset