6.1.19 11:59 C:\Users\svenm\Docum...\JankyII_dokmit.m 1 of 4
% Jánky Gibsonova metoda na celém přivaděči L=1800 m
% 2 2tryskové horizontální Peltonovy turbíny
% 1 tlakový senzor, druhý se bere hladina na vtoku, kolísání hladiny
% neuvažováno
clc
close all clear all format compact
%%
load data_jankyII % načtení měření
c_pst=9328.981; % pentstock factor =suma Li/Ai [m-1]
delta_t=1; % časový krok [s]
L=1715; % délka úseku mezi tlakovými senzory [m]
Q0=0; % průsak po dovření RK [m3/s]
ro=998; % hustota vody [kg/m3]
g=9.81; % gravitační zrychlení [m/s-2]
t=data_jankyII(:,1); % čas s
P1=data_jankyII(:,2); % výkon turbiny č. 1 kW P2=data_jankyII(:,3); % výkon turbiny č. 2 kW
p=data_jankyII(:,4); % tlak před turbínou m v. sl.
tr11=data_jankyII(:,5); % otevření trysky 1 tubiny 1 tr12=data_jankyII(:,6); % otevření trysky 2 tubiny 1 tr21=data_jankyII(:,7); % otevření trysky 1 tubiny 2 tr22=data_jankyII(:,8); % otevření trysky 2 tubiny 2
figure;
plot(t,P1); %vizuální výběr doby integrace hold on
plot(t,P2);
plot(t,p);
plot(t,tr11);
plot(t,tr12);
plot(t,tr21);
plot(t,tr22);
legend('P1','P2','p','t11','t12','t21','t22')
% legend('P1','p','t11') grid on
legend show
% title({'p'});
% xlabel('t [s]');
% ylabel('p [m v. sl.]');
%%
Tstart = 1146; % doba manipulace Tstop = 1159;
udalost_start=t(1); % časový rozsah dat
6.1.19 11:59 C:\Users\svenm\Docum...\JankyII_dokmit.m 2 of 4
udalost_stop=t(length(t));
TPstart = Tstop; % doba pro stanovení p_offset, zároveň dokmit TPstop = udalost_stop;
p_offset_vyber=p(t>Tstop); % interval pro stanovení statického tlaku t_offset_vyber=t(t>Tstop);
P1_vyber=P1((t>Tstart) & (t<Tstop)); % výkon turbiny č. 1 kW P2_vyber=P2((t>Tstart) & (t<Tstop)); % výkon turbiny č. 2 kW
p_offset=mean(p_offset_vyber) % statický tlak jako průměr tlaku po uzavření
figure;
plot(t_offset_vyber,p_offset_vyber);
title({'p offset'});
xlabel('t [s]');
ylabel('p [m v. sl.]');
grid on
dp=(p-p_offset)*ro*g; % převod na Pa
figure;
plot(t,dp,'g');
title({'dp'});
xlabel('t [s]');
ylabel('p [Pa]');
dp_vyber = dp((t>Tstart) & (t<TPstop)); % t_vyber = t((t>Tstart) & (t<TPstop)); % grid on
figure;
plot(t_vyber,dp_vyber,'r') title({'dp vyber'});
xlabel('t [s]');
ylabel('p [Pa]');
grid on
p_poc_vyber=dp(t<Tstart);
t_poc_vyber=t(t<Tstart);
p_poc=mean(p_poc_vyber)
dp_vyber_orez = dp((t>Tstart) & (t<Tstop))./(ro*g); % checkpoint Q=Q0 t_vyber_orez = t((t>Tstart) & (t<Tstop)); % oříznuté pro výběr dp ořez
%%
Q_odhad=0.1; % odhad počátečního průtoku
Q_konc=Q_odhad; % pro spuštění cyklu Q_odhad musí být většínez Q0 k=1; % koeficient tlumení iterace, aby nedošlo k oscilaci dummy=0; % čítač iterací
%%
while ((abs(Q_konc-Q0))> 0.00001) % dokud je větší, cyklus probíhá
6.1.19 11:59 C:\Users\svenm\Docum...\JankyII_dokmit.m 3 of 4
R=-(1/(Q_odhad^2))*p_poc; % úprava R na aktuální Q_odhad
Q(1)=Q_odhad-1/(ro*c_pst)*R*abs(Q_odhad)*Q_odhad*delta_t-(1/(ro*c_pst))
*dp_vyber(1)*delta_t;
Hzt(1)=R*abs(Q_odhad)*Q_odhad; % ztráty, abs kvůli oscilaci
for n=2:length(dp_vyber) % je jedno, jestli výběr nebo celý
Q(n)=Q(n-1)-1/(ro*c_pst)*R*abs(Q(n-1))*Q(n-1)*delta_t-(1/(ro*c_pst))*dp_vyber (n)*delta_t;
Hzt(n)=R*abs(Q(n-1))*Q(n-1);
end
Q_konc=Q(length(dp_vyber_orez)); % koncový průtok musí být Q0 Q_odhad=Q_odhad-(Q_konc-Q0)*k;
dummy = dummy + 1;
Q_iterace(dummy,:)=Q; % vykreslení iterace Q
end
%%
Q_vyber_delka=length(t_vyber_orez);
Q_vyber=Q(1:Q_vyber_delka);
% figure;
% plot(t_vyber_orez,Q_vyber);
% hold on
% plot(t_vyber_orez,P1_vyber);
% plot(t_vyber_orez,P2_vyber);
% plot(t_vyber_orez,dp_vyber_orez);
% title({'souvislosti'});
% xlabel('t [s]');
% ylabel('Q [m3/s]');
figure;
plot(t_vyber,Q_iterace);
title({'Q_iterace'});
xlabel('t [s]');
ylabel('Q [m3/s]');
grid on
pocet_iteraci=dummy
konec=length(dp_vyber);
zacatek=length(dp_vyber_orez)+1;
Q_dokmit=Q(zacatek:konec); % výběr pro zprůměrování koncového průtoku
% figure;
% plot(Q_dokmit);
% title({'Q dokmit'});
% xlabel('t [s]');
% ylabel('Q [m3/s]');
6.1.19 11:59 C:\Users\svenm\Docum...\JankyII_dokmit.m 4 of 4
Q_kontrola=mean(Q_dokmit)
Hzt=Hzt/(ro*g);
figure;
plot(Q,Hzt);
title({'Hzt'});
xlabel('Q [m3/s]');
ylabel('Hzt [m v. sl.]');
grid on
figure;
plot(t_vyber,Q) title({'Q'});
xlabel('t [s]');
ylabel('Q [m3/s]');
grid on
dp_vyber=dp_vyber/1000000;
P1=P1/1000;
tr11=tr11/1000;
figure;
plot(t_vyber,dp_vyber,'r');
hold on
plot(t_vyber,Q);
plot(t,P1);
plot(t,tr11);
title({'dokmit'});
xlabel('t [s]');
ylabel('p [Pa]');
legend('dp','Q','P1','tr11') grid on
Q_prubeh=Q;
Q=Q_odhad % ustálený průtok před začátkem zavírání m3/s