Corso di Infrastrutture Idrauliche ing Stefano Alvisi stefano

Corso di Infrastrutture Idrauliche ing. Stefano Alvisi stefano. alvisi@unife. it Esercitazione 1 Sistema acquedottistico di adduzione

Dato il sistema acquedottistico di adduzione schematizzato in figura si vuole dimensionare la condotta di adduzione in acciaio che possa convogliare la portata Qr in condizioni di regime dal serbatoio di monte al serbatoio di valle e che possa sopportare idraulicamente i fenomeni di moto vario indotti da una valvola a galleggiante posta in prossimità del serbatoio di valle. L’elevata capacità del serbatoio di monte consente di assumere che la quota del pelo libero dell’acqua in esso contenuta zm non vari significativamente nel tempo. Il serbatoio di valle ha forma cilindrica con area di base Av ed il pelo libero dell’acqua in esso contenuta zv risulta generalmente variabile nel tempo in relazione alla portata entrante ed uscente dal serbatoio stesso. La portata Q convogliata dalla condotta è regolata automaticamente in funzione della quota zv secondo la curva caratteristica della valvola a galleggiante.

Dati del problema Qr=0. 4 m 3 s-1 zm=62 m s. l. m. zvr=40 m s. l. m. L=1500 m Av=100 m 2 portata di regime quota pelo libero di monte quota pelo libero di valle a regime lunghezza condotta area serbatoio di valle Parametri delle perdite distribuite =2 mm Coefficienti di perdita di carico concentrate 1=0. 5 all’ imbocco della condotta 2=0. 3 al gomito in 2 3=0. 3 al gomito in 3 4 funzione del grado di apertura della valvola 5=1 allo sbocco della condotta

Dimensionamento della condotta: formulazione generale

Dimensionamento della condotta: formulazione generale

Dimensionamento della condotta: assunzione di moto assolutamente turbolento

Dimensionamento della condotta: assunzione di moto assolutamente turbolento

% ESERCITAZIONE 1 : sistema acquedottistico di adduzione clear all %Dati Qr=0. 4; zm=62; zvr=40; regime L=1500; Av=100; xi=[0. 5 0. 3 concentrata % [m^3 s^-1] % [m] slm quota del serbatoio di valle a % [m] lunghezza della condotta % [m^2] area del serbatoio di valle 0 1]; % coefficienti di perdita di carico %xi 1=>imbocco; xi 2, xi 3=>spigoli; xi 4=>valv ola; xi 5=>sbocco epsilon=0. 002; % [m] mu=1. 206 e-6; % [m g=9. 81; % [m Vmin=1; % [m Vmax=7; % [m scabrezza tubazione esercizio s^-2] viscosità cinematica s^-2] s^-1] velocità minima in condotta s^-1] velocità massima in condotta

%Parte prima: %determinazione del diametro della condotta Dmin=sqrt((4*Qr)/(pi*Vmax)); % [m] diametro min. della condotta Dmax=sqrt((4*Qr)/(pi*Vmin)); % [m] diametro max. della condotta D=[Dmin: (Dmax-Dmin)/1000: Dmax];

%Calcolo la funzione phi per ogni valore di D for cont=1: length(D); % ciclo sui diametri %Calcolo iterativo del lambda per ciascun D %inizializzo il ciclo while con moto ass. turbolento lambda 0=(2*log 10(3. 71*D(cont)/epsilon))^(-2); Re=… err=1; while err>=0. 001 lambda 1=. . err=abs(lambda 1 -lambda 0); lambda 0=lambda 1; end lambda. D(cont)=lambda 0; end phi=zm-zvr-((lambda. D. /D)*…

Dd=interp 1(phi, D, 0) %interpolo la funzione phi per calcolare il diametro di progetto figure(1), plot(D, phi), grid on title('Funzione Phi(Q)') xlabel('diametri [m]'), ylabel('Phi(Q)') Dc=…; % Diametro commerciale scelto

%Calcolo ora le portate che otterrei con il Dc Lambda 0=. . . Q 0 es=… Re=… err=1; while err>=0. 001 vv. prima ma attenzione: Q 0 es=… Re=… end lambdac=… Q 0 es=…;