EQUAZIONI CONCATENATE DI ALLIEVI 2 2 i i1 i1 i1 i i z z 2 2 Al( z z) 0 2 / ) ( h t h z i i 0 / ) ( i i t 0 0 2 gh aU Al 0 0 0 a: celerità g: accelerazione di gravità U : velocità iniziale di efflusso h :altezza piezometrica iniziale all'efflusso ω :area iniziale dell'otturatore
16
Embed
EQUAZIONI CONCATENATE DI ALLIEVI. c Al=numero di Allievi (aU0/2gh0) c a=celirità dell'onda di pressione c U0=velocità iniziale di efflusso c g=accelerazione.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
EQUAZIONI CONCATENATEDI ALLIEVI
2 2i i 1 i 1 i 1 i iz z 2 2 Al ( z z )
02 /)( hthz ii
0/)( ii t
0
0
2gh
aUAl
0
0
0
a: celeritàg: accelerazione di gravitàU : velocità iniziale di efflussoh :altezza piezometrica iniziale all'efflussoω :area iniziale dell'otturatore
c Al=numero di Allievi (aU0/2gh0)
c a=celirità dell'onda di pressione
c U0=velocità iniziale di efflusso
c g=accelerazione di gravità
c h0=carico iniziale all'efflusso
c teta=tempo di fase (2*L/a)
c L=lungheza della condotta
c tc=tempo di chiusura dell'otturatore
c tempo=istante di tempo finale
c*************************************************** allievi.for ***************************************************c questo programma partendo da un istante iniziale t1 calcola mediante la formula iterativa di c Allievi l’incremento relativo della quota piezometrica all'otturatore, per diversi istanti di tempo c distanziati fra loro per intervalli pari alla durata di fase teta (2*L/a).c**********************************************************************************************************************Diagramma di flusso del programma
START
Definiscigrandezze
Dati di Input
Listato del programma in linguaggio Fortran
real e,tt,t0,ti,tempo,t1,tc,teta,Xi,U integer n,j,i,ii common tc
c richiesta dati write(*,*)'inserire il tempo di fase (2*L/a)?' read(*,*)teta write(*,*)'inserire il tempo di chiusura?' read(*,*)tc write(*,*)'inserire il numero di Allievi (aU0/2gh0)?' read(*,*)Al
c inizio ciclo per l'inserimento del tempo iniziale
do 3 ii=1,10
write(*,*) 'inserire il tempo iniziale?'
write(*,*)'(numero compreso fra 0 e Θ)'
read(*,*)t1
c calcolo il tempo finale
tempo=t1+((n-1)*teta)
c valori iniziali del carico X(0)
X0=1
t0=0
Diagramma di flusso del programma
Apri file di output
Inizio sequenza di calcolo
delle soluzioniconcatenate
per il genericotempo iniziale t1
Listato del programma in linguaggio Fortran
c calcolo di n (numero di intervalli) if (teta.ge.tc) then
n=8 else n=NINT((tc*8)/teta) endif c apro il file dei risultati open(5,file='out8') write(5,*)' Al=', Al write(5,*)' teta=',teta
write(5,*)' tc=',tc write(5,*)' n=',n write (5,*)
Calcolo intervalli n
seΘ>Tc
ovveroΘ<Tc
c scrittura su file del tempo iniziale e finalewrite(5,*)write(5,*)write(5,*)write(5,*)' t1=',t1 write(5,*)'tempo finale=',tempo
c srittura su file delle condizioni iniziali su carico, apertura e velocitàwrite(5,*)write (5,*)' t h/h0 e U/U0'
c write (5,*)' 0 1 1 1 '
c formula iterativado 2 i=1,nti=t1+(i-1)*teta
c controllo numerico val=(Al*e(ti))**2-(X0)**2+2+(2*Al*e(t0)*X0)
if (val.lt.0.0001) thenwrite(*,*)‘ERRORE NUMERICO!!!!'goto 100else goto 200endif