Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék

11. FEJEZET. PÉLDAPROGRAMOK 44
function nyiltfelszin
global B n g Q ii
nyiltfelszin.m
Q = 0.298; % eloirt terfogataram
y0= 0.6; % szint a csatorna elején
ii = 0.001; % hidraulikai lejtes
n = 0.01; % csoerdesseg
B = 0.5; % csatorna szelesseg
g = 9.81; % gravitacio
L = 1000; % csohossz
Rh = y0*B/(2*y0+B);
C = Rh^(1/6)/n;
Qn = y0*B*C*sqrt(ii*Rh);
fprintf(’\n Nyiltfelszinu csatornaaramlas\n’);
fprintf(’\n %g m-es cs}oeleji felszinhez tartozó térfogatáram: %g m^3/s’,y0,Qn);
% Megoldo hivasa:
options=odeset(’events’,@rohanas);
[x,y]=ode45(@nyiltfelszin_de,[0 L],y0,options);
% A csatorna alja és az ahhoz képesti folydékszint (rajzoláshoz):
for i=1:length(x)
yy(i)=ii*(1-x(i)/L)*L;
yf(i)=yy(i)+y(i); v(i) =Q/y(i)/B;
end
figure(1)
subplot(2,1,1), plot(x,yy,x,yf);
title([’Nyiltfelszin}u csatornaáramlás: i=’,num2str(ii),’
ylabel(’folyadekfelszin, y [m]’,’FontSize’,12)
legend(’csatorna feneke’,’folyadék felszíne’)
axis([0 L 0 yf(1)*1.2])
Q=’,num2str(Q)],’FontSize’,12);
subplot(2,1,2), plot(x,v);
xlabel(’csohossz, x [m]’,’FontSize’,12);
ylabel(’aramlasi sebesseg, v [m/s]’,’FontSize’,12);
%---------------------------------------------
function dydx=nyiltfelszin_de(x,y)
global B n g Q ii
A=y(1)*B; Rh=y(1)*B/(2*y(1)+B); C=Rh^(1/6)/n;
dydx=[ (ii-Q^2/A^2/C^2/Rh)/(1-Q^2*B/A^3/g)];
%---------------------------------------------
function [val,ter,dir]=rohanas(x,y)
val=y(1); % Az esemény definiciója
ter=1; % Állítsuk meg az integrálást.
dir=0; % Mindkét irányban (1: + -> -, 2: - -> +)