% DANE OGOLNE
rho = 1000;   %kg/m3
g = 9.81;
Ec = 2030e6; % modul sprezystosci cieczy WODA
Er = 3200e6; % modul sprezystosci rurociagu PE [Pa]
lambda = 0.02; % czy nie powinno tobyæ wyliczone 64/Re?
pb = 101300; % cisnienie barometryczne
pv = 2338;   % cisnienie preznosci

Q = 500/3600;
pmin = pb-pv; %minimalne cisnienie w rurociagu

%wsp. strat miejscowych
Z1 = 0;
Z2 = 0;
Z3 = 0;
Z4 = 0;
Z5 = 0;

% Dlugosci
L1 = 1000;
L2 = 1000;
L3 = 1000;
L4 = 1000;
L5 = 1000;

% Rzedne wezlów
RzW_1 = 0;
RzW_2 = 0;
RzW_3 = 50;
RzW_4 = 100;
RzW_5 = 150;
RzW_6 = 200;
RzW_7 = 250;

% Srednica
D1 = 0.267;
D2 = 0.2906;
D3 = 0.3106;
D4 = 0.3274;
D5 = 0.3412;

% Scianka
e1 = 0.0665;
e2 = 0.0547;
e3 = 0.0447;
e4 = 0.0363;
e5 = 0.0294;

% Przekroj
A1 = 3.14*D1^2/4;
A2 = 3.14*D2^2/4;
A3 = 3.14*D3^2/4;
A4 = 3.14*D4^2/4;
A5 = 3.14*D5^2/4;

% Objetosci
Ob1 = L1*A1;
Ob2 = L2*A2;
Ob3 = L3*A3;
Ob4 = L4*A4;
Ob5 = L5*A5;

% Wsp. masy A^2/m
Am1 = A1/L1/rho;
Am2 = A2/L2/rho;
Am3 = A3/L3/rho;
Am4 = A4/L4/rho;
Am5 = A5/L5/rho;

% Roznica rzednych
H1 = RzW_3 - RzW_2;
H2 = RzW_4 - RzW_3;
H3 = RzW_5 - RzW_4;
H4 = RzW_6 - RzW_5;
H5 = RzW_7 - RzW_6;

% Cisnienie od ciezkosci
pg1 = H1*rho*g;
pg2 = H2*rho*g;
pg3 = H3*rho*g;
pg4 = H4*rho*g;
pg5 = H5*rho*g;

% Rezystancja (od predkosci) Zastepczy wsp strat
Rv1 = 8*rho/3.14^2/D1^4*(lambda*L1/D1 + Z1);
Rv2 = 8*rho/3.14^2/D2^4*(lambda*L2/D2 + Z2);
Rv3 = 8*rho/3.14^2/D3^4*(lambda*L3/D3 + Z3);
Rv4 = 8*rho/3.14^2/D4^4*(lambda*L4/D4 + Z4);
Rv5 = 8*rho/3.14^2/D5^4*(lambda*L5/D5 + Z5);

% Straty
ps1 = Rv1*abs(Q)*Q;
ps2 = Rv2*abs(Q)*Q;
ps3 = Rv3*abs(Q)*Q;
ps4 = Rv4*abs(Q)*Q;
ps5 = Rv5*abs(Q)*Q;

% Cisnienie poczatkowe
p5 = pg5 + ps5;
p4 = pg4 + ps4 + p5; 
p3 = pg3 + ps3 + p4; 
p2 = pg2 + ps2 + p3; 
p1 = pg1 + ps1 + p2; 


% Sprezystosc objetosciowa ciecz+rura
KV1 = 1/(Ob1*(D1/Er/e1+1/Ec));
KV2 = 1/(Ob2*(D2/Er/e2+1/Ec));
KV3 = 1/(Ob3*(D3/Er/e3+1/Ec));
KV4 = 1/(Ob4*(D4/Er/e4+1/Ec));
KV5 = 1/(Ob5*(D5/Er/e5+1/Ec));

% Zmiana objetosci segmenu
DV1 = p1/KV1;
DV2 = DV1 + p2/KV2;
DV3 = DV2 + p3/KV3;
DV4 = DV3 + p4/KV4;
DV5 = DV4 + p5/KV5;




