% (Tsentrifugaal-)pumba survekarakteristiku arvutused
% Q - jõudlus (tootlikkus) [m3/h]*0.28=[l/s]
% http://www.ene.ttu.ee/elektriajamid/teadus/artiklid/jjoller/sm_tasuvus.html
n0=1450; % nimipöörlemiskiirus [1/min]
n=1150; % pöörlemiskiirus [1/min]
A=pi*(0.025^2); % toru ristlõige [m2]
v0=29; %voolu kiirus [m/s]
g=9.81; %[m/s2]
roo=998.2; %vedeliku tihedus [kg/m3]
nyy=0.68; %pumba kasutegur(antud tootlikuse juures)*mootori kasutegur

Q0=A*v0; %vooluhulk [m3/s]
Q1=A*v0*3600 %vooluhulk [m3/h]
Q2=A*v0*1000 %vooluhulk [l/s]
Q=Q0*(n/n0) %vähendatud kiirusel [m3/s]

H0=(v0^2)/(2*g) %tõstekõrgus [m] 
H=H0*(n/n0)^2 %vähendatud kiirusel

P0=(Q0*H0*g*roo)/nyy % vajalik mootori nimivõimsus [W] (11.19)
P=P0*(n/n0)^3 %vähendatud kiirusel
nyy_sm=nyy*0.9*0.96; %pumba-*mootori-*sagedusmuunduri kasutegur vähendatud n-ga
dH=30 % tõstekõrguste vahe (pumba sise-ja väljundrõhu vahe ?)
P_sm=(Q*dH*g*roo)/(nyy_sm) %Tarbitav võimsus sagedusreguleerimisel

% H = Hst_staatiline_surve(kõrguste vahe)+ ht_survekadu_võrgus
k=2500; % torustiku voolutakistust iseloomustav kordaja
Hst=15; ht=k*Q0^2;  % [m]

Q_=0.01:0.01:0.09; %[m3/s] c 50...350m3/h
h_=Hst+k*Q_.^2; %torustiku QH-karakteristik
h2_=Hst+1.5*k*Q_.^2; %suurendatud takistusega torustiku QH-karakteristik
H_=50-2500*Q_.^2; %pumba QH-karakteristik
H_n=32-2500*Q_.^2; %pumba QH-karakteristik vähendatud n-ga
P_=0.001*(Q_.*H_*g*roo)/nyy; %[kw]PQ tunnusjoon
P_n=0.001*(Q_.*H_n*g*roo)/nyy; %[kw]PQ tunnusjoon vähendatud n-ga
plot(Q_*3600,h_,Q_*3600,P_,Q_*3600,H_,Q_*3600,H_n,Q_*3600,P_n);grid %[m3/h]
legend('torustiku QH','PQ','pumba QH','pumba QH vähendatud n','PQ vähend n')
xlabel('Q [m3/h]'); ylabel('H [m], P [kW]')
title('Pumba ja torustiku QH ning PQ karakteristikud') 
% figure, plot(Q_*3600,h_,Q_*3600,h2_,Q_*3600,H_,Q_*3600,H_n);grid %[m3/h]