Hydraulika - Krzysztof Książyński
Hydraulika - Krzysztof Książyński
Hydraulika - Krzysztof Książyński
HYDRAULIKA
ZESTAWIENIE POJĘĆ I WZORÓW STOSOWANYCH W BUDOWNICTWIE
Kraków 2008
3
SPIS TREŚCI
1. WPROWADZENIE .......................................................................................................... 7
1.1. Słowo wstępne ..................................................................................................... 7
1.2. Wykaz oznaczeń .................................................................................................. 9
1.3. Oznaczenia pomocnicze .................................................................................... 11
2. HYDROSTATYKA ........................................................................................................ 12
2.1. Ciśnienie ............................................................................................................ 12
2.1.1. Pojęcie ciśnienia ................................................................................. 12
2.1.2. Prawo Eulera ...................................................................................... 12
2.1.3. Prawo Pascala..................................................................................... 12
2.1.4. Hydrostatyczny rozkład ciśnienia ....................................................... 13
2.1.5. Ciśnienie względne i bezwzględne ..................................................... 13
2.1.6. Wysokość ciśnienia ............................................................................ 14
2.1.7. Obliczanie ciśnień w naczyniach połączonych ................................... 14
2.1.8. Przenoszenie ciśnienia zewnętrznego do wnętrza cieczy ................... 15
2.1.9. PrzełoŜenie prasy hydraulicznej ......................................................... 16
2.2. Parcie hydrostatyczne na powierzchnie płaskie ................................................. 18
2.2.1. Parcie hydrostatyczne ......................................................................... 18
2.2.2. Wykres przestrzennego rozkładu parcia ............................................. 19
2.2.3. Bryła parcia ........................................................................................ 19
2.2.4. Środek parcia...................................................................................... 20
2.2.5. Równowaga sił z udziałem parcia ...................................................... 22
2.2.6. Rozkład parcia na składowe ortogonalne ........................................... 23
2.3. Parcie na ściany krzywoliniowe ........................................................................ 24
2.3.1. Składowe ortogonalne parcia na ściany krzywoliniowe ..................... 24
2.3.2. Bryła parcia poziomego...................................................................... 24
2.3.3. Bryła parcia pionowego...................................................................... 25
2.3.4. Wyznaczanie parcia elementarne powierzchnie krzywoliniowe ......... 25
2.3.5. Kierunek działania wypadkowej parcia na ściany w kształcie
powierzchni obrotowych .................................................................... 26
2.3.6. Wyznaczanie parcia na złoŜone powierzchnie krzywoliniowe ........... 27
2.4. Pływanie ciał ..................................................................................................... 29
2.4.1. Prawo Archimedesa ............................................................................ 29
2.4.2. Równowaga ciał pływających ............................................................ 29
2.4.3. Podstawowe pojęcia teorii pływania .................................................. 30
2.4.4. Stateczność ciał całkowicie zanurzonych w cieczy ............................ 31
2.4.5. Stateczność ciała pływającego po powierzchni cieczy ....................... 32
2.4.6. Wyznaczanie wysokości metacentrycznej .......................................... 34
3. HYDRODYNAMIKA ..................................................................................................... 36
3.1. Opis ruchu płynu ............................................................................................... 36
3.1.1. Ciągłość płynu .................................................................................... 36
3.1.2. Prędkość lokalna ................................................................................ 36
4
1. WPROWADZENIE
v – wektor prędkości,
v – prędkość średnia,
vo – prędkość początkowa,
v – wartość prędkości,
vh – prędkość pozioma,
vv – prędkość pionowa,
vx – składowa prędkości w kierunku x,
vy – składowa prędkości w kierunku y,
vz – składowa prędkości w kierunku z.
12
2. HYDROSTATYKA
2.1. CIŚNIENIE
pa pa
p p p’
ciśnienie ciśnienie
bezwzględne zewnętrzne ciśnienie
względne
γh pa γh
h
h
Rys. 1
Znaczna wartość ciśnienia atmosferycznego w porównaniu ze skalą ciśnień
spotykanych na co dzień powoduje, Ŝe wygodniej jest się posługiwać ciśnieniami
względnymi, tzn. liczonymi względem ciśnienia atmosferycznego:
14
p’ = γ h ,
np. ciśnienie słupa wody o wysokości 1 [m] wynosi:
– bezwzględne: p = pa + γ h = 1 013 [hPa] + 98.1 [hPa] = 1 112 [hPa]
( ≡ 1.0968 [atm]),
– względne: p’ = γ h = 98.1 [hPa] ( ≡ 0.0968 [atm]).
h=p/γ.
Na ogół podaje się wysokość słupa wody (ciecz najczęściej spotykana) lub rtęci
(najcięŜsza). Ten sposób podawania wartości ciśnienia ułatwia obliczenia
hydrauliczne, gdyŜ nie wymaga przeliczania masy na cięŜar.
γ1
h1
h2
γ2
A A
B ∆h
Rys. 2
pB = γ1 h1 + γ2 ∆h = γ2 (h2 + ∆h),
czyli:
γ2 h2 = γ1 h1 ,
które jest podstawą rozwiązywania zadań z naczyń połączonych.
2.1.8. Przenoszenie ciśnienia zewnętrznego do wnętrza cieczy
Ciśnienie pB na powierzchni tłoka spręŜającego ciecz w cylindrze (rys. 3)
odpowiada ciśnieniu na głębokości hP = pB / γ pod powierzchnią swobodnego
zwierciadła cieczy. Takie rozumowanie pozwala wyznaczyć połoŜenie tzw.
16
z w i e r c i a d ł a p i e z o m e t r y c z n e g o , czyli p o z o r n e g o z w i e r c i a d ł a
cieczy (porównaj rozdział 4.3.2).
p
F hP
pB
B
A γ h
C pC
pB γh
h
Rys. 3
Ciśnienie poniŜej tłoka oblicza się z zaleŜności (rys. 3):
pC = pB + γ h, pB = F / A,
stąd:
pC = F / A + γ h .
4
17
F2
F1
d2
d1
Rys. 4
Zatem p r z e ł o Ŝ e n i e p r a s y h y d r a u l i c z n e j wynosi:
2
F2 d 2
= .
F1 d1
18
p/γ
h 0
θ ξ
SP h
H P z
H xs
•
• x
Ξ
γ
S
ys
Ξ η
y
ηo
h
Rys. 5
Ciśnienie hydrostatyczne cieczy na ścianę stanowi obciąŜenie ciągłe, które
moŜna, w przypadku obliczania reakcji działających na podpory ściany, zastąpić
działaniem skupionej s i ł y p a r c i a o wartości:
P = ps A ,
gdzie ps – ciśnienie w środku geometrycznym S ściany o polu powierzchni A (rys.
5). Siła parcia P na ścianę płaską jest wektorem o kierunku prostopadłym do
ściany i zwrocie od cieczy ku ścianie, zaczepionym w punkcie Ξ, noszącym nazwę
środka parcia.
19
p 0
γ1 h1 P1
p1 = γ1 h1
γ2 h2 P2
p2 = p1 + γ2 h2
γ3 h3 P3
p3 = p2 + γ3 h3
y
γ3 > γ2 > γ1
Rys. 6
b a
SP h
ys
b a
Rys. 7
– wyznaczanie analityczne:
h 2a+ b
ys =
3 a+b
21
ys =
∑ y i Vi ,
V
gdzie: Vi, V – objętości składowych brył elementarnych i całej bryły,
yi, , ys – współrzędne środków geometrycznych składowych brył
elementarnych
i całej bryły.
Analityczne wyznaczanie współrzędnych środka parcia
(ξ, η ; rys. 5):
Jx
η= ,
Mx
gdzie: J x = ∫Ay
2
dA – moment bezwładności ściany względem krawędzi
zwierciadła cieczy x,
y = z / sin θ – współrzędna danego punktu ściany (odległość od
krawędzi zwierciadła x),
Mx = A ys – moment statyczny powierzchni ściany względem
krawędzi zwierciadła wody x.
Zastosowanie twierdzenia Steinera do wyznaczania momentu
bezwładności podczas obliczania połoŜenia środka parcia prowadzi do
wniosku, iŜ środek ten musi zawsze znajdować się poniŜej środka
geometrycznego powierzchni:
Ay s + J o
2
Jx J
η= = = ys + o > ys .
Mx Ay s Mx
gdzie: Jo
– moment bezwładności powierzchni ściany
względem osi równoległej do x , przechodzącej
przez jej środek geometryczny S.
Jedyny wyjątek stanowią ściany poziome, dla których środek parcia pokrywa
się ze środkiem geometrycznym.
Jo
η= .
Mx
γb
+
γw
SB h
H η
P G
ξ O
F
b
a
Rys. 8
Parcie cieczy stanowi jeden z rodzajów obciąŜeń budowli i konstrukcji
wodnych, decydujący o ich stateczności. Budowla moŜe stracić stateczność na
skutek przewrócenia się lub poślizgu po podłoŜu. ObciąŜenia budowli moŜna
ogólnie podzielić na naruszające i podtrzymujące jej równowagę. Do
podstawowych obciąŜeń naruszających równowagę naleŜy zazwyczaj parcie wody
P, do podtrzymujących – cięŜar budowli G (rys. 8).
Stateczność ze względu na przewrócenie:
Działanie sił naruszających równowagę moŜe spowodować obrót wokół
pewnej osi, której śladem na rysunku jest punkt O. NaleŜy zatem sprawdzić
równowagę momentów sił względem O, przyjmując za dodatni np. moment
przewracający:
MO = P ξ – G η .
23
P Pz
Py θ1
θ2 Px
Rys. 9
W przypadku ścian ukośnych dogodnie jest, ze względów obliczeniowych,
wyznaczać wypadkową parcia P na podstawie jej składowych ortogonalnych Px,
Py, Pz (rys. 9), tzn. pokrywających się z osiami układu współrzędnych
kartezjańskich o pionowej osi z. Trzem składowym parcia odpowiadają trzy bryły
parcia. W przypadku ścian płaskich układ współrzędnych dobiera się tak, by siła
parcia miała tylko dwie składowe – poziomą Ph = Px i pionową Pv = Pz
i odpowiednio dwie bryły parcia. Zasada wyznaczania tych brył dla składowych
jest taka sama jak w przypadku parcia na ściany krzywoliniowe (rozdziały 2.3.2
i 0).
24
z
Pz
P
Py
Ξ Px
x y
Rys. 10
B r y ł a p a r c i a p o z i o m e g o jest to (trójwymiarowy)
wykres rozkładu wysokości ciśnienia działającego na rzut danej
powierzchni na płaszczyznę pionową.
25
Sh Ph
Rys. 11
B r y ł a p a r c i a p i o n o w e g o jest to (trójwymiarowy)
wykres rozkładu wysokości ciśnienia działającego na rzut danej
powierzchni na płaszczyznę poziomą.
Najwygodniej jest wykorzystać do konstrukcji bryły płaszczyznę
zwierciadła cieczy – wysokość ciśnienia stanowi wtedy po prostu odległość danej
powierzchni krzywoliniowej od zwierciadła (rys. 11). Dla płaskiego zwierciadła
bryła parcia pionowego jest zawarta pomiędzy powierzchnią ściany,
płaszczyzną zwierciadła i pionowymi tworzącymi wystawionymi na
krawędziach ściany. W przypadku gdy zwierciadło nie jest płaszczyzną (np.
podczas ruchu cieczy) bryła parcia jest ograniczona od dołu powierzchnią ściany, a
od góry powierzchnią zwierciadła. Wyznaczonemu na podstawie bryły parciu
pionowemu nadaje się zwrot od cieczy ku ścianie.
2.3.4. Wyznaczanie parcia elementarne powierzchnie krzywoliniowe
Oznacza to, Ŝe dowolna linia pionowa lub pozioma moŜe przeciąć powierzchnię
elementarną tylko w jednym punkcie (rys. 12). Bryły parcia poziomego
i pionowego dla takiej powierzchni wyznacza się na zasadach podanych
w rozdziałach 2.3.2 i 0. Na podstawie brył parcia wylicza się wszystkie elementy
(wartość, zwrot i kierunek działania) wektorów parcia poziomego
i pionowego, po czym dodaje się je wektorowo (za pomocą wieloboku sił)
uzyskując wypadkową parcia (rys. 9):
Pz Px
tg θ1 = , tg θ 2 = .
Px2 +Py2 Py
= +
= +
Rys. 12
P
Pv O
Ph
Ξ
Rys. 13
B-B
A-A
B A
Rys. 14
W prezentowanym przykładzie bryła parcia poziomego jest połową ściętego
walca o osi poziomej przechodzącej przez środek kuli. Natomiast bryła parcia
pionowego jest złoŜona z dwóch części – półkuli (reprezentującej parcie
skierowane w dół) i połowy walca o osi pionowej przechodzącej przez środek kuli
z wyciętą jej ćwiartką (reprezentującej parcie skierowane w górę). Dla parcia
29
F = W – G = γ V – γc V = (γ – γc) V = (ρ – ρc) g V ≠ 0.
NaleŜy przy tym odróŜniać cięŜar objętościowy ciała, który jest wyliczany dla
całej objętości V, od jego cięŜaru właściwego, który moŜe być róŜny dla róŜnych
części ciała niejednorodnego. Jeśli gęstość ciała będzie większa od gęstości cieczy
(ρ c > ρ), ciało będzie tonąć, osiągając równowagę dopiero na dnie zbiornika,
gdzie siła wypadkowa zostanie zrównowaŜona przez reakcję podłoŜa. Ciało
o takiej samej gęstości jak ciecz (ρ c = ρ) będzie się w niej swobodnie unosiło nie
zmieniając połoŜenia.
Ciało o gęstości mniejszej od gęstości cieczy (ρ c < ρ) wypłynie na
powierzchnię, gdzie na skutek wynurzenia zmniejszy się objętość wypieranej przez
nie cieczy. W miarę wynurzania wartość wyporu będzie maleć aŜ do uzyskania
równowagi pomiędzy wyporem i cięŜarem. Ciało będzie zatem pływało swobodnie
po powierzchni cieczy a objętość jego zanurzonej części będzie funkcją głębokości
zanurzenia – V(h). Z warunku równowagi sił moŜna tę głębokość wyliczyć
podstawiając taką jej wartość (h), przy której część zanurzona osiągnie wymaganą
objętość (G/γ):
V (h ) = .
G
G = W = γ V(h) →
γ
Dla obiektu o płaskim dnie i pionowych ścianach bocznych rozwiązanie uzyskuje
się wprost:
V (h ) = h A =
G G
→ h= ,
γ γ A
gdzie A oznacza powierzchnię dna. Objętość V konieczna dla uzyskania
równowagi jest róŜna w cieczach o róŜnych gęstościach, róŜne są więc głębokości
zanurzenia. Na statkach pełnomorskich dopuszczalne zanurzenie w róŜnych
warunkach wyznacza tak zwany z n a k P l i m s o l l a ( z n a k w o l n e j
b u r t y ) a jego wartość średnią – z n a k ś r e d n i e j l i n i i w o d n e j .
2.4.3. Podstawowe pojęcia teorii pływania
W W W
S
Ξ Ξ ΞΞ
Ξ Ξ Ξ
S G
Ξ ΞS Ξ
Ξ Ξ
Ξ
G G
równowaga obojętna Ξ trwała
równowaga równowaga chwiejna
Rys. 15
z
Ξ
M
m
Ξ
Ξ W’
S θ W1
a
Ξ Ξ ’
W2 Ξ
’ Ξ Ξ’
Ξ r
Ξ y
Ξ Ξ
G
W Ξ
Ξ
Rys. 16
G M
S M Ξ W S
Ξ Ξ S W Ξ
M Ξ
Ξ Ξ Ξ
Ξ W Ξ Ξ G
G
Ξm > 0Ξ Ξ Ξ
m=Ξ
0 m<0
równowaga obojętna równowaga trwała równowaga chwiejna
Rys. 17
Moment sił W’ i G jest momentem prostującym dla m > 0 (M powyŜej S,
równowaga trwała), przewracającym dla m < 0 (M poniŜej S, równowaga
chwiejna), zaś nie występuje dla m =0 (M pokrywa się z S, równowaga obojętna).
34
b
M
m Ξ
Ξ
Ξ W’
S θ Q
Ξ Ξ
Ξ’
Ξ
G
Ξ
Rys. 18
Z przekształcenia tej równości otrzymuje się wzór na wysokość metacentryczną:
Qb
m= ,
Wθ
gdzie W = G + Q – całkowita wyporność wagowa obiektu. Wzór jest poprawny
dla niewielkich kątów przechyłu θ, rzędu kilku stopni (zobacz rozdz. 2.4.5), a
wtedy ramię siły Q – b moŜe być mierzone tak jak pokazano na rys. 18.
36
3. HYDRODYNAMIKA
1
v = ∫A v dA .
A
Pojęcie prędkości średniej v w przekroju strumienia A jest dla hydrauliki pojęciem
podstawowym, poniewaŜ pozwala na określenie przepływu Q przy jedno-
wymiarowym opisie ruchu:
Q=vA
(porównaj rozdział 3.1.3).
U w a g a : W następnych rozdziałach będzie stosowana jedynie prędkość średnia,
dla uproszczenia oznaczana jako v .
3.2.2. Lepkość
L e p k o ś ć jest to zdolność cieczy do stawiania oporu podczas
ruchu jej elementów względem siebie.
Miarą lepkości jest współczynnik lepkości zaleŜny od rodzaju, temperatury
i ciśnienia cieczy.
K i n e m a t y c z n y w s p ó ł c z y n n i k l e p k o ś c i ν jest to
stała proporcjonalności pomiędzy siłami tarcia wewnętrznego
w ruchu laminarnym a poprzecznym do kierunku ruchu
przyrostem prędkości cieczy.
Taka proporcjonalność, znana jako h i p o t e z a N e w t o n a , zapisywana jest
wzorem:
ν ∂v
τ = ,
ρ ∂y
gdzie napręŜenie ścinające τ reprezentuje siłę tarcia wewnętrznego przypadającą
na jednostkę powierzchni tarcia [N/m2 = Pa], zaś pochodna ∂v/∂y określa
zmienność prędkości cieczy w kierunku poprzecznym do ruchu. Współczynnik
lepkości cieczy maleje ze wzrostem temperatury (T ). Dla wody obowiązuje
zaleŜność:
1.78 [mm 2 /s]
ν = .
1 + 0.0337 [1/ o C] T + 0.000221 [1/ o C 2 ] T 2
A1
v1
v2
A2
Rys. 19
Przy obliczaniu przepływu w strumieniu o przekroju kołowym, np. w rurociągu
otrzymuje się zaleŜność:
π d2
Q =v A=v = const ,
4
a stąd:
2
d
v 2 = v 1 1 .
d2
Równanie ciągłości w tej postaci umoŜliwia m.in. obliczanie prędkości
przepływu w róŜnych przekrojach rurociągu o zmiennych średnicach (ruch
niejednostajny ustalony). Analogicznie oblicza się prędkości w zmieniających się
przekrojach koryt otwartych.
W ruchu nieustalonym poprzeczne rozmiary strumienia mogą zmieniać się
w czasie i równanie ciągłości musi być poszerzone o człon retencyjny
(gromadzenie cieczy w czasie przepływu między przekrojami).
p1 α1 v1 p2 α2 v2
z1 + + = z2 + + + hstrat = const .
γ 2g γ 2g
z
[linia energii całkowitej]
α1 v 1 2
linia energii
[mechanicznej] hstrat
2g
p1 linia ciśnień α 2 v 22
piezometrycznych
γ 2g
p2
v1 γ
v2
z1
z2
[poziom porównawczy]
x
Rys. 20
Zgodnie z prawem zachowania
energii suma wysokości energii
p α v 2
grawitacyjnej (z), spręŜystej , kinetycznej i cieplnej (h ) jest
strat
γ
2 g
wielkością stałą (linia energii całkowitej), czyli jest taka sama w kaŜdym przekroju
poprzecznym strumienia (rys. 20). Wielkość hstrat jest wysokością energii
rozproszonej w otoczeniu w postaci ciepła (i dźwięku) w trakcie pokonywania
oporów tarcia pomiędzy przekrojami 1 i 2.
Ek(v) = α Ek (v) .
Stąd wynika wzór na współczynnik:
∫A v
3
E (v ) dA
α= k = 3
.
E k (v ) v A
Współczynnik Saint-Venanta dla ruchu turbulentnego w rurociągach wynosi
α ≅ 1, dla laminarnego w przekroju kołowym α = 2, w ciekach skanalizowanych
α ≅ 1.1, w ciekach naturalnych moŜe sięgać do α = 2.5.
46
4. HYDRAULIKA RUROCIĄGÓW
4.1. DEFINICJE
W efekcie na części jego długości panuje podciśnienie (p < pa). Ruch w lewarze
wymaga zapoczątkowania.
S y f o n jest to rurociąg, którego oś przebiega na pewnym
odcinku poniŜej wylotu.
Woda pozostaje na tym odcinku nawet po przerwaniu zasilania. Na całej długości
syfonu panuje nadciśnienie (p > pa).
W przypadku, gdy oś rurociągu przebiega na całej długości pomiędzy
zwierciadłem wody w zbiorniku zasilającym a poziomem wylotu, ciśnienie na jego
długości jest niewiele wyŜsze od ciśnienia atmosferycznego.
c) e)
a)
b) d)
f)
Rys. 21
S t r a t y e n e r g i i k i n e t y c z n e j to straty powstające na
końcu rurociągu w wyniku unoszenia pędu przez wypływający
płyn.
Energia ta moŜe być wykorzystana mechanicznie, ale ostatecznie ulega zamianie
na ciepło na skutek rozproszenia w otoczeniu (w powietrzu lub w zbiorniku
dolnym).
Wzór Weissbacha
Wszystkie straty energii są proporcjonalne do prędkości przepływu.
PoniewaŜ miarą strat jest wysokość energii mierzona w [m], zatem przyjmuje się,
Ŝe wysokość strat jest proporcjonalna do wysokości energii kinetycznej:
α v2
h =ζ ,
2g
gdzie: ζ – współczynnik strat. Obliczenia strat naleŜy przeprowadzać zgodnie
z normą PN-76/M-34034. Rurociągi. Zasady obliczeń strat ciśnienia.
W s p ó ł c z y n n i k c h r o p o w a t o ś c i b e z w z g l ę d n e j k jest to średnia
wysokość chropowatości na ściankach wewnętrznych rury. Dla stalowych rur
wodociągowych i Ŝeliwnych kanalizacyjnych wynosi ona około 1.5 [mm], dla rur
betonowych około 2.5 [mm].
Współczynnik chropowatości względnej rur:
k
ε= .
d
W s p ó ł c z y n n i k t a r c i a r u r λ zaleŜy od chropowatości rur i panującej
w rurociągu liczby Reynoldsa. Jego wartość określa empiryczny wzór
Colebrooka-White’a:
1
λ= 2
.
ε
2.51
2 log10 +
3.71 Re λ
Wzór ma jednak charakter uwikłany i w praktyce λ wyznacza się za pomocą
wykresu Moody’ego (norma PN-76/M-34034 załącznik 3).
W y k r e s M o o d y ’ e g o jest to wykres zaleŜności λ(Re,ε) przedstawiony
w układzie logarytmicznym z poziomą osią Re i pionową λ .Tworzy go szereg
krzywych dla róŜnych wartości ε. Na wykresie moŜna wyróŜnić pięć stref
zmienności współczynnika λ:
– s t r e f a r u c h u l a m i n a r n e g o (Re < 2300), w której współczynnik λ
zaleŜy (zgodnie ze wzorem Hagena-Poisseile’a) tylko od liczby Re:
64
λ= ;
Re
– s t r e f a p r z e j ś c i o w a (2300 < Re < 4000) o nieokreślonym reŜimie ruchu,
w której współczynnik λ wyznacza się w zaleŜności od rzeczywistego reŜimu
panującego w rurociągu zgodnie z regułami dla stref sąsiednich. W obliczeniach
projektowych przyjmuje się istnienie reŜimu mniej korzystnego;
52
R – promień hydrauliczny,
d – średnica rury
Sf = hstrat / l – spadek linii energii,
l – długość rurociągu.
n – współczynnik szorstkości, przyjmujący dla rurociągów wartości od
0.011 (stalowe) do 0.014 [s/m1/3] (betonowe).
Z porównania ze wzorem Darcy’ego-Weissbacha wynika zaleŜność:
8 3 4 g n2 124.5 [m/s 2 ] n 2
λ = 3
=
3
.
d d
Ze względu na wymiar współczynnika szorstkości średnicę (promień hydrauliczny)
naleŜy podstawiać w [m]. Stosowanie wzoru Manninga do obliczania strat
w rurociągach dla całego obszaru ruchu turbulentnego nie jest zgodne z obecnie
obowiązującymi normami.
linia ciśnień v2
v1 II
l1
x
Rys. 22
2
v2
H= + hw + h1 + hR + h2 + hL + h3 + h< + h4 + hz + h5 ,
2g
gdzie: h1, h2, h3, h4, h5 – straty na długości odcinków l1 do l5 ,
hw – strata na wlocie,
hR – strata na łuku,
hL – strata na zmianie kierunku,
h< – strata na poszerzeniu,
hz – strata na zmianie kierunku;
5. Podział na rysunku wysokości H zgodnie z równaniem na poszczególne straty
w kolejności ich występowania;
6. Wyznaczenie na wykresie połoŜenia linii energii w kolejnych przekrojach
pionowych poczynając od poziomu zasilania. Przekroje te lokalizuje się na
końcach odcinków lub w tych punktach osi rurociągu, w których występują
straty lokalne;
7. Wyznaczenie połoŜenia linii ciśnień równolegle do linii energii z obniŜeniem
α v2
w pionie o wartość wysokości prędkości 2g
w danym przekroju.
– równanie ciągłości:
4Q
vi = ,
π di 2
– wzór na współczynnik tarcia:
v i di k
λi(Rei, εi ) , Rei = , εi = ,
ν di
– wzór na współczynnik strat lokalnych:
ζj = f(dj, typ armatury);
j oznacza kolejny numer straty lokalnej. Razem 1 + 2 i + j niewiadomych
(H, vi, λi, ζj) i tyleŜ równań. Prędkość końcowa vk jest jedną z liczby i .
3. Rozwiązanie:
– obliczenie prędkości średnich vi w odcinkach o stałych średnicach di ,
– obliczenie współczynników strat lokalnych ζj z odpowiednich tablic lub
wzorów,
– obliczenie współczynnika tarcia rur λi na podstawie wykresu Moody’ego lub
wzorów,
– obliczenie wysokości energii zasilania przez podstawienie wyliczonych
wielkości do wzoru Bernoulliego.
Rozwiązanie uzyskuje się zatem przez proste podstawienie zadanych wartości do
wzorów.
2gH
vk = 2
,
2
l dk d
∑ λi di + ∑ ζ j dk
i di j
π dk 2
– obliczenie wydatku: Q = v k .
4
– PrzybliŜenia następne: współczynnik tarcia λ(Re,ε) oblicza się dla przepływu
wyliczonego w poprzednim przybliŜeniu:
– obliczenie prędkości vi na podstawie wcześniej wyliczonego przepływu,
– obliczenie współczynników tarcia rur λi na podstawie wykresu
Moody’ego lub wzorów,
– obliczenie prędkości końcowej i nowej wartości wydatku jak
w przybliŜeniu początkowym.
– Kryterium zakończenia iteracji („kryterium stopu”):
jeśli wartości współczynników tarcia λi w kolejnych przybliŜenia róŜnią się
od siebie nie więcej niŜ o 3 %, obliczenia zostają zakończone, a wartość
przepływu przyjmuje się z ostatniego przybliŜenia.
Innym sposobem rozwiązania jest metoda prób i błędów. Polega ona na
rozwiązywaniu zadania typu H dla kolejnych przybliŜeń wartości wydatku
i porównywaniu wyniku z zadaną wartością H. Wartość początkową Q moŜna
przyjąć np. na podstawie prędkości ekonomicznej.
– PrzybliŜenia następne:
– przyjęcie wartości pozostałych średnic di na podstawie kryteriów
pozahydraulicznych,
– rozwiązanie zadania typu H przy załoŜonych średnicach,
– porównanie wyliczonej wartości H’ z wartością zadaną H,
– przyjęcie nowej wartości średnicy podstawowej d’’ według zasady:
Jeśli H’ > H , to d’’ > d’, jeśli H’ < H , to d’’ < d’.
Średnice dobiera się na podstawie wyliczenia spośród wartości
nominalnych, podanych w programie produkcji rur stosowanych
w rurociągu lub według PN-83/H-02651.
– Kryterium zakończenia iteracji:
jeśli dla sąsiadujących w wykazie średnic nominalnych róŜnica H’ - H
zmienia znak, jako ostateczną przyjmuje się średnicę większą. Takie
rozwiązanie zapewnia uzyskanie dla zadanej wysokości energii zasilania
przepływu w rurociągu nie mniejszego niŜ zadany.
61
h
p min
γ
Rys. 23
Obliczenia przepływu
W lewarze (rozdział 4.1.2 oraz rys. 23), jak w kaŜdym rurociągu, przepływ
odbywa się na skutek występowania spadku wysokości hydraulicznej pomiędzy
wlotem i wylotem. Spadek ten, odpowiadający wysokości energii zasilania H,
wywołuje w rurociągu o średnicy d przepływ o wielkości Q. Zatem konieczne jest
przeprowadzenie dla lewara takich obliczeń hydraulicznych jak dla innych
rurociągów, czyli rozwiązanie zadania typu H, Q lub d.
Obliczenie ciśnienia w głowicy
Ze względu na występowanie wewnątrz lewara podciśnienia, konieczne jest
sprawdzenie, czy w jakimś punkcie lewara ciśnienie nie spadło poniŜej punktu
wrzenia cieczy w danej jej temperaturze. Przy znanym wzniesieniu h ponad
zwierciadło w górnym zbiorniku wysokość ciśnienia p (bezwzględnego) w danym
przekroju lewara wyniesie (zgodnie z równaniem Bernoulliego):
p pa pw
= − h − hstrat > .
γ γ γ
Sprawdzenie naleŜy przeprowadzić dla głowicy lewara, czyli punktu
o najwyŜszej wysokości z = h (najniŜej w stosunku do osi rurociągu połoŜonym
zwierciadle piezometrycznym). Jak wynika ze wzoru, panuje tam najniŜsze
ciśnienie. Wysokość strat dotyczy odcinka pomiędzy wlotem a głowicą. Normalna
wysokość bezwzględnego ciśnienia atmosferycznego pa /γ wynosi 10.33 [m],
62
4Q
gdzie: v i = .
π di 2
Równanie Bernoulliego
– dla szeregowego układu gałęzi: spadek wysokości hydraulicznej (linii ciśnień)
jest równy sumie strat na długości w poszczególnych gałęziach:
H= ∑ hi ,
– dla równoległego układu gąłęzi: spadki linii ciśnień na wszystkich gałęziach
łączących dwa węzły (A i B) są równe i odpowiadają stracie na długości kaŜdej
z gałęzi:
H = HA – HB = hi ,
– dla układu gałęzi tworzących oczko: spadek linii ciśnień na zamkniętym
obwodzie oczka (tzn. po powrocie do punktu wyjścia) wynosi zero:
H = ∑ hi = 0 .
Wynika to z równości strat na dwu hipotetycznych gałęziach tworzących oczko;
w przypadku obliczania strat w kierunku przeciwnym niŜ przepływ wysokość
ciśnienia rośnie, a straty te są traktowane jako ujemne.
Analogiem równania Bernoulliego w elektrotechnice jest prawo Ohma dla sieci
elektrycznych.
Qo Qk
Qw
Rys. 24
64
H1
H2 H
Hp H3
H4
Q1 Q4 Q3
Q2
H5
Qw Q5
Rys. 25
H2
Q1 H4
H3
Q3
Q2
Rys. 26
Obliczenie przepływów w sieci równorzędnych rurociągów jest zadaniem
bardziej złoŜonym. Rozkład przepływów zaleŜy wtedy od rozdziału w węzłach
sieci. Najprostszy przypadek stanowi zadanie z pojedynczym węzłem, jest to tzw.
zadanie „trzech zbiorników” (rys. 26). Zasada postępowania jest następująca:
1. Postawienie zadania:
Dane: d1, d2, d3, H1, H2, H3.
Szukane: Q1, Q2, Q3, H4 .
2. Równania podstawowe:
– Równanie ciągłości dla w = 1 węzłów:
Q1 + Q2 = Q3 ,
– Równanie Bernoulliego dla n = 3 gałęzi:
2
l1 v 1
H1 − H 4 = λ1 itp.,
d1 2 g
67
4 Qi
gdzie: vi = , łącznie w+n=1+3=4 równania i tyleŜ
π di 2
niewiadomych.
3. Rozwiązanie:
Z zestawienia wynika, Ŝe układ równań jest rozwiązywalny i w zasadzie moŜna
tego dokonać metodą podstawień eliminujących kolejne niewiadome, jest to
jednak procedura trudna i pełna zasadzek numerycznych. W praktyce
rozwiązanie najłatwiej uzyskać metodą iteracyjną:
– PrzybliŜenie początkowe: załoŜenie wysokości hydraulicznej (H4) w węźle
wewnętrznym sieci:
H1 > H2 = H4 > H3 ,
czyli przyjęcie braku przepływu w gałęzi o pośredniej wysokości ciśnienia,
przez co unika się konieczności ustalenia kierunku ruchu;
– PrzybliŜenia następne:
– Obliczenie wydatku w dwóch (n – 1) gałęziach (przy znanej wartości
spadku jest to dla kaŜdej gałęzi zadanie typu Q – rozdz. 4.4.4):
π d1 2 2 g (H1 − H 4 )
Q1 =
4 l
λ1 1
d1
i analogicznie w gałęzi nr 2;
– Obliczenie wydatku w ostatniej gałęzi z równania ciągłości:
Q 3 = Q 1 + Q 2;
– Określenie wysokości hydraulicznej w węźle (przy znanej wartości
wydatku jest to zadanie typu H):
2
l3 v 3 4 Q3
H 4 ' = H 3 + λ3 , v3 =
d3 2 g π d32
moŜna teŜ ocenić H4’ na podstawie Q3 w przybliŜeniu;
– Powtórzenie iteracji dla H4 = H4’;
– Kryterium stopu: jeśli wyliczone w kolejnych iteracjach wydatki Q3 róŜnią
się mniej niŜ o 3 %, obliczenia zostają zakończone, a jako wynik przyjmuje
się ostatnio wyliczony układ wydatków i wysokości ciśnienia.
Ze względu na iteracyjny przebieg obliczeń współczynników λi mogą być
wyznaczane dla poprzednio wyliczonych wydatków (w pierwszym przybliŜeniu
dla rur chropowatych) bez konieczności wprowadzania dodatkowych przybliŜeń.
68
H2
Q2
Q Q2 Q
Q3
Rys. 27
Rurociągi równoległe (rys. 27, rozdz. 4.1.3) to najprostszy rodzaj sieci
pierścieniowej. Węzeł rozdzielczy moŜe być połączony z węzłem zbiorczym
dowolną liczbą (n) gałęzi, które tworzą od jednego do n – 1 oczek. Podstawowym
celem obliczeń jest wyznaczenie rozdziału przepływów pomiędzy gałęziami. Jeśli
jest on znany, zadanie moŜna podzielić na kilka prostych zadań (rozdz. 4.4)
odnoszących się do poszczególnych gałęzi. Rurociągi równoległe to jedyny
przypadek sieci pierścieniowej, który moŜe być rozwiązywany bez uŜycia
skomplikowanych metod iteracyjnych. Przebieg obliczeń jest następujący:
1. Postawienie zadania:
Dane: Q, H2, di, li, k, ν
Szukane: Qi, H1
2. Równania podstawowe:
– Równanie ciągłości (identyczne dla obu węzłów):
Q = ΣQi ,
69
Bi B
hi hi+1 1 1:m
h θ
m
b
Rys. 27
∑ (hi + hi +1 ) Bi ,
1
A=
2
72
– trapezowego:
A = b h + m h2 .
O b w ó d z w i l Ŝ o n y jest to ta część obwodu przekroju
poprzecznego, na którym występuje wymiana pędu
z otoczeniem.
Na powierzchni obwodu zwilŜonego występuje tarcie, w wyniku którego lokalna
prędkość wyraźnie spada. Ma to miejsce przede wszystkim na dnie i ścianach
koryta, gdzie prędkość spada do zera. Natomiast nie wlicza się do obwodu
szerokości swobodnego zwierciadła B.
Wartość obwodu zwilŜonego oblicza się wzorami:
– w korycie naturalnym:
N
U= ∑ (hi +1 − hi ) + Bi2 ,
2
i =1
– w korycie trapezowym:
U = b + 2 h 1 + m2 .
K o r y t o p r y z m a t y c z n e to takie, w którym zmienia się
jedynie napełnienie h, natomiast geometria dna (szerokość
i nachylenie skarp) pozostaje stała.
K o r y t o c y l i n d r y c z n e to takie, którego przekrój daje się
opisać funkcją ciągłą.
Przykładem takiego koryta jest koryto o przekroju kołowym lub parabolicznym.
v12
H1 2g
z1 h1 Sf hstrat
S
2
v2
H2
So 2g
h2
l
z2
θ
Rys. 28
Wartości spadków podaje się w promilach [‰ = m/km] lub w procentach
[%]. Dla potoków górskich mogą one wynosić nawet kilkanaście procent, dla rzek
nizinnych zwykle dziesiąte części promila. Z tego względu odległość przekrojów
( l ) moŜna mierzyć w planie, co pozwala wyznaczać ją na podstawie mapy.
3
R2 Sf 1 23 1
v= = R Sf 2 .
n n
gdzie: R – promień hydrauliczny [m];
spadek Sf podaje się w wartościach bezwymiarowych.
P r o m i e ń h y d r a u l i c z n y R określany jest na ogólnych zasadach (rozdz.
3.2.3) i w płytkich, a szerokich ciekach odpowiada on głębokości średniej:
A bh
R= ≅ → h (dla b → ∞) ,
U b+2h
nigdy tej wielkości nie przekraczając. W rzeczywistości, przy praktycznie
występujących proporcjach wymiarów koryt cieków, promień hydrauliczny spełnia
warunki:
h
≤ R < h.
2
U3
U1
n3
U2
Rys. 29
n=
∑Ui ni
=
U1 n1 + U 2 n 2 + U 3 n3
.
U U1 + U 2 + U 3
Wzór Darcy’ego-Weissbacha
W dokładniejszych obliczeniach uwzględniających przepływy na terasach
porośniętych wysoką roślinnością stosowane jest tzw. u n i w e r s a l n e p r a w o
p r z e p ł y w u oparte na wzorze Darcy’ego-Weissbacha (rozdział 4.2.3):
1
v = 8 g R Sf ,
λ
gdzie współczynnik tarcia λ dla koryt nie porośniętych określa się ze wzoru
Prandtla-Nikuradse’go (rozdz. 4.2.3):
1
λ= 2
.
k
2 log10
14.84 R
n2 n1 n3
h
II h2 h3 III
h1 I
b2 b3
R
b1
Rys. 30
W korytach otwartych o złoŜonym przekroju poprzecznym, szczególnie
w korytach terasowych, które słuŜą do przeprowadzenia wielkich wód (rys. 30),
obliczenia prowadzone dla całego przekroju mogą doprowadzić do uzyskania
skokowego spadku przepływu przy wzroście głębokości. Jest to efekt nieciągłości
promienia hydraulicznego przy napełnieniu odpowiadającym poziomowi terasy.
Następuje tu nagły przyrost obwodu zwilŜonego przy niewielkim przyroście
powierzchni. Koryta o takiej charakterystyce muszą być traktowane jako koryta
wielodzielne. Podział na części obliczeniowe przeprowadza się w ten sposób, by
oddzielić od siebie fragmenty przekroju, na których występują wyraźnie róŜne
prędkości przepływu.
K o r y t o w i e l o d z i e l n e jest to koryto, którego przekrój
został podzielony liniami pionowymi na części obliczeniowe
w celu uniknięcie nieciągłości w obliczeniach przepływu
(rys. 30).
Podział koryta ma uzasadnia się występowaniem wyraźnie róŜnych prędkości
przepływu w jego poszczególnych częściach. W praktyce linie podziału wystawia
się na krawędziach teras.
77
Q2
So = 2
;
3
Ri
2
∑ Ai ni
h o”
ho”’
h o’
Q
Q’ Q Q”
Rys. 31
– Zadanie typu h – obliczenie napełnienia kanału:
dane: Q, So, szukane: ho.
Ze względu na uwikłanie zmiennej ho we wzorach opisujących geometrię
przekroju (A, U, R) obliczenie polega na iteracyjnym zastosowaniu wzoru na
wydatek dla kolejnych przybliŜeń wartości ho . Wyniki przedstawia się zwykle
w tabeli:
Lp. ho U A R n v Q
1 1m
2
Obliczony przepływ porównuje się z zadanym i na tej podstawie przyjmuje się
nowe przybliŜenie napełnienia ho zakładając proporcjonalność Q ~ ho . MoŜna
posłuŜyć się przy tym pomocniczym wykresem (rys. 31). Iteracje kończy się, gdy
kolejne przybliŜenia h róŜnią się nie więcej niŜ o 1 %, a jako wynik przyjmuje
się wartość dającą wydatek bliŜszy zadanemu.
Ro 2 Ro
3 Ro
Ro Ro / 2 Ro
Rys. 32
Moduł przepływu i moduł prędkości
Przepływ w całkowicie wypełnionym przekroju zamkniętym jest funkcją
spadku dna So oraz rodzaju zastosowanego przekroju. Rodzaj przekroju określa
jego wymiary geometryczne i materiał ścian, a więc współczynnik szorstkości n.
Tak więc wzór na prędkość i przepływ moŜna zapisać w postaci:
vo = w So , Qo = K So ,
gdzie w i K są charakterystyczne dla danego rodzaju przekroju:
– moduł prędkości:
3
R2
w = ,
n
– moduł przepływu:
3
R2
K =A =w A .
n
Dla kolektora betonowego (n = 0.014 [s/m1/3]):
– ∅ 500 [mm]: K = 3.50 [m3/s],
– ∅1000 [mm]: K = 22.2 [m3/s].
80
– Zadanie typu So :
dane: Q, h, rozmiar kolektora H; szukane: So
– obliczenie napełnienia względnego ξ = h/H,
– wyznaczenie z tablic współczynnika sprawności FQ dla danego napełnie-
nia ξ ,
– obliczenie wydatku przy całkowitym napełnieniu: Qo = Q / FQ,
– wyznaczenie modułu K dla danego kolektora z tablic,
2
Q
– obliczenie spadku kolektora: So = o ;
K
– Zadanie typu h :
dane: Q, So, rozmiar kolektora H; szukane: h
– wyznaczenie z tablic modułów K i w dla danego kolektora,
– obliczenie Qo : Qo = K So ,
– obliczenie współczynnika sprawności FQ: FQ = Q / Qo,
– wyznaczenie z tablic napełnienia względnego ξ = h/H na podstawie
wartości FQ,
– obliczenie napełnienia: h = ξ H.
g hk
vk = .
α
H1 = H 2 +
S1 + S 2
2
l−
α
2g
v1 − v 2 .
2 2
( )
Ostatni człon wzoru, wyraŜający zamianę energii kinetycznej na potencjalną,
uwzględnia się tylko w ciekach o skoncentrowanym przekroju (rzeki, kanały).
W zbiornikach, w których wykształcają się wewnętrzne strumienie, podczas gdy
pozostała masa wody pozostaje niemal nieruchoma, człon kinetyczny pomija się
w obliczeniach. Zakłada się, Ŝe energia prędkości zostaje rozproszona na skutek
tarcia wewnętrznego.
K o l e j n o ś ć o b l i c z e ń jest następująca:
1. Wstępne załoŜenie nieznanej wartości stanu H1 w przekroju górnym. MoŜna ją
określić z warunku: z1 + ho < H1 < z1 + h2 .
85
∆H
H ∆H H H
h
h
Rys. 33
W y p ł y w p r z e z o t w ó r w c i e n k i e j ś c i a n c e to taki,
w którym strumień odrywa się od ściany na krawędzi otworu
(rys. 34).
Zjawisko to zachodzi przy grubości ścianki l < (3.5 ÷ 4) H .
l
l
H
Rys. 34
W y p ł y w p r z e z p r z y s t a w k ę to wypływ, w którym
strumień przylega do ściany otworu na całej grubości ścianki
(rys. 34).
Wypływ taki zachodzi przy grubości ścianki (3.5 ÷ 4) H < l < (6 ÷ 7) H .
Powstaje wtedy efekt ssania zwiększający wydatek. Przy większej grubości ścianki
straty na długości powodują konieczność obliczania takiego wypływu jak
rurociągu pod ciśnieniem.
W y p ł y w p r z e z o t w ó r m a ł y to taki, w którym zmiany
ciśnienia na wysokości otworu są niewielkie i moŜna je
pominąć przy obliczaniu wypływu.
MoŜna przyjąć, Ŝe ma to miejsce dla H > 4.5 e (rys. 35) i dla otworów w ścianie
poziomej (np. w dnie).
90
W y p ł y w p r z e z p r z e l e w o o s t r e j k r a w ę d z i to
taki, w którym strumień odrywa się od ściany na krawędzi
przelewu (rys. 37).
Zjawisko to zachodzi przy grubości ścianki l < (0.1 ÷ 0.5) H .
Wypływ przez przelew o kształcie
p r a k t y c z n y m to taki, w którym strumień przylega do
korony przelewu, ulegając przy tym silnemu zakrzywieniu.
Zjawisko to zachodzi przy grubości ścianki (0.1 ÷ 0.5) H < l < 2 H .
W y p ł y w p r z e z p r z e l e w o s z e r o k i e j k o r o n i e to
taki, w którym na koronie przelewu występuje płaski odcinek
strumienia (rys. 39).
Zjawisko to zachodzi przy grubości ścianki 2 H < l < (8 ÷ 15) H. Większe
wartości przyjmuje się dla gładszych progów. Przy większej długości (l ) korony
przelewu wypływ naleŜy obliczać jak przepływ w kanale otwartym (ze stratami na
długości) z poziomem wody na wlocie jak dla progu o szerokiej koronie.
Typ wypływu przez przelew nie zaleŜy od geometrii przelewu, a tylko od
wysokości słupa wody ponad jego koroną. Ten sam przelew moŜe być w róŜnych
warunkach kanałem lub przelewem o ostrej krawędzi.
H L>3e
zupełna
L>3b L<3b
e
zupełna
brak
niezupełna
kontrakcja strumienia
b
Rys. 35
W s p ó ł c z y n n i k k o n t r a k c j i ( d ł a w i e n i a ) ε jest to
stopień zmniejszenia powierzchni przekroju strumienia za
upustem wskutek kontrakcji:
As
ε= ,
Ao
gdzie As i Ao to odpowiednio powierzchnia przekroju strumienia i powierzchnia
czynna upustu. Z definicji wynika ograniczenie ε < 1 .
Warunki wystąpienia kontrakcji
92
y
2
vo
2g
θ
Ho H1 Ho
H
H2
b(y)
vo
Rys. 36
Prędkość wypływu cieczy przez otwór niezatopiony moŜna określić
z równania Bernoulliego (rys. 36):
2 g Ho
v= =ϕ 2 g Ho ,
1+ζ
93
gdzie Ho oznacza wysokość linii energii nad punktem, dla którego liczona jest
prędkość:
2
vo
Ho = H + ,
2g
zaś ζ stanowi współczynnik strat na wlocie otworu.
Współczynnik p r ę d k o ś c i ϕ jest to stopień
zmniejszenia prędkości wypływu na skutek strat lokalnych na
otworze:
v 1
ϕ= = .
vt 1+ζ
Q = As v = µ Ao 2 g Ho ,
Q = As v = µ z Ao 2 g ∆H o ,
α vo 2
gdzie: ∆H o = ∆H + 2g
.
Współczynnik wydatku ulega zmniejszeniu w stosunku do wypływu swobodnego:
µz = 0.986 µ .
Przy wypływie spod zasuwy bez progu, czyli dla otworu przydennego,
energia strumienia moŜe być na tyle duŜa, by odrzucić wodę dolną od otworu.
Głębokość wody w odpływie wynika wtedy z kontrakcji strumienia w otworze
i wynosi:
h = ε e,
gdzie e oznacza wysokość otworu (otwarcie zasuwy), natomiast współczynnik
kontrakcji jest zbliŜony do wartości ε = 0.62 . Współczynnik wydatku dla zasuwy
z wlotem o ostrych krawędziach wynosi µ = 0.58 .
(
Q = m b 2 g H 23 2 − H13 2 , )
gdzie: m = 2 µ / 3 – współczynnik otworu.
Współczynnik wydatku przybiera wartości od µ = 0.65 dla otworów średniej
wielkości z kontrakcją zupełną do µ = 0.85 dla otworów duŜych praktycznie bez
kontrakcji.
Przy w y p ł y w i e p o d t o p i o n y m dokonuje się podziału strumienia na
część swobodną i zatopioną (o wysokości h, rys. 33) opisując wypływ przez kaŜdą
z nich odpowiednimi wzorami:
3 3
2 v o2 2 v o2 2 v o2
Q = b 2 g µ H 2 − h + − H1 +
+ µz h H2 − h +
3 2 g 2 g 2g
.
Współczynniki wypływu dla obu części są prawie równe i moŜna przyjąć:
96
µ = µz = 0.65 .
Dla w y p ł y w u z a t o p i o n e g o rozkład ciśnienia wypadkowego (∆p)
i prędkości wypływu jest stały, a więc poprawny jest wzór dla małego otworu
(rozdział 6.2.2).
W praktycznych obliczeniach wydatku duŜego otworu często wykorzystuje
się wzór na wypływ swobodny z małego otworu uzupełniony o współczynnik
zatopienia σ :
Q =σ m Ao 2g H o ,
6.3. PRZELEWY
gdzie: m = 2
3
εϕ = 2
3
µ – współczynnik przelewu.
Często stosowany jest skrócony zapis współczynników wzoru w postaci:
97
C =m 2g .
cg l
cd
Rys. 37
Ze względu na występowanie w sąsiedztwie przelewu ruchu przyspie-
szonego (rozdział 3.2.1) pojawia się depresja (obniŜenie) zwierciadła, która ma
zasięg rzędu (3÷5) H . Pomiar głębokości H wody na przelewie naleŜy zatem
wykonywać w odległości L ≅ 3 H (rys. 37).
0.4 > b > 2.0 [m], 0.1 < H < 0.6 [m], 0.2 < cg < 2.0 [m], 0.1 < b/Bo < 1.0 .
W przypadku tzw. strugi nieswobodnej wydajność przelewu, a więc i m wzrastają.
Q = σ ε m b 2 g Ho
32
.
θ
bo
a
Rys. 38
θ
0.996
Q = 1.331 tg H 2.47 .
2
Wzór jest bardzo zbliŜony do teoretycznej całki (0.996 ≅ 1, 2.47 ≅ 5/2) i daje
dobre wyniki dla kąta 20o < θ < 118o.
P r z e l e w T o m p s o n a – jest to przelew trójkątny o kącie rozwarcia θ = 90o.
Do obliczeń stosowany jest wzór empiryczny:
Q = 1.343 H 2.47 .
Wzór daje dobre wyniki dla 0.05 < H < 0.25 [m].
P r z e l e w t r a p e z o w y – dla kąta nachylenia skrzydeł θ i długości korony
poziomej bo stosuje się wzór:
Q = 0.42 ε ( bo + 0.8 H ctg θ ) 2 g H 1.5 .
Wzór daje dobre wyniki dla szerokości koryta doprowadzającego Bo < 6.68 d
i cg > 0.2 d .
P r z e l e w p r o p o r c j o n a l n y ( R e t t g e r a ) – jego kształt zapewnia liniową
zaleŜność wydatku od głębokości:
Q = C H,
przy szerokości opisywanej funkcją:
100
b(y ) =
2C
, µ ≅ 0.61 .
π µ 2g y
ε = ∑b ,
Bo
ξ p + (N − 1) ξ f H o
ε = 1 − 0. 2 ,
N b
gdzie ξp i ξf to współczynniki kształtu odpowiednio przyczółków i filarów, zaś N
to ilość przęseł przelewu. Dla prostokątnego kształtu części wlotowej przyczółków
i filarów ξp = 1, ξf = 0.8, przy wlocie zaokrąglonym ξp = 0.7, ξf = 0.45, dla filaru
opływowego (z krawędzią tnącą) ξf = 0.25. Wzór moŜna stosować przy
Bo > N (b + d) i stopniu zatopienia h/Ho ≤ (0.85÷0.90); d oznacza tu szerokość
filara. Jeśli warunki te nie zostaną spełnione, stosuje się wzór uproszczony:
Ho
ε = 1 − 0. 2 ξ f .
b
H Ho
h
cg hz
cd
Rys. 39
Potencjalną głębokość na progu określa się za pomocą wzoru: h = k Ho .
Parametr k wyznacza się z tablic lub wykresów opracowanych na podstawie badań
empirycznych.
102
Wypływ przez próg niezatopiony oblicza się podobnie jak dla progu
o kształtach praktycznych:
Q=ε mb
32
2 g Ho .
Q = ε ϕ hz b 2 g (H o − hz ) .
α v 12
2g
H v1
h1 Ho
ho
h v
Rys. 40
Przepływ pod mostem ma charakter przepływu przez próg o szerokiej
koronie przy zerowej wysokości progu (cg = 0). Dla większych mostów opory
ruchu są jednak na tyle małe, Ŝe moŜna je pominąć i prowadzić obliczenia przy
załoŜeniu stałej wysokości linii energii H w otoczeniu mostu (rys. 40).
Obliczenie ś w i a t ł a m o s t u polega na wyznaczeniu sumarycznej
szerokości b przęseł koniecznej dla przepuszczenia spodziewanego przepływu Q
w cieku przy prędkości v pod mostem, nie przekraczającej prędkości
dopuszczalnej vd z uwagi na rozmycie dna. Ze względów ekonomicznych zakłada
się maksymalną dopuszczalną prędkość v = vd . Dla piasku o średnicy ziarn
103
2 g Ho
vk =
3 α
i w oparciu o relację pomiędzy vk a vd określić reŜim przepływu.
1.05 Q
b= .
Ho −
α v d2
vd
2g
α Q2 α vd 2
H ≡ h1 + =3 ,
2 g B 2 h12 2g
e
Ho
H
vo
cd
h v
Rys. 41
106
P r ę d k o ś ć v w d o l n y m s t a n o w i s k u moŜna wyznaczyć
z równania Bernoulliego analogicznie jak dla upustów (rozdz. 6.1.6) ze wzoru:
v =ϕ 2 g (H o − h ) .
hs
h ho
Rys. 42
Powrót do nadkrytycznej głębokości normalnej w dolnym stanowisku
zachodzi w tzw. odskoku hydraulicznym. Przybiera on róŜne formy w zaleŜności
107
H ho
H hs
hs c ho c
hs ho
szykany
Rys. 43
W praktyce inŜynierskiej odskok hydrauliczny projektuje się jako zatopiony,
który jako krótszy wymaga mniejszej niecki wypadowej. Cel ten moŜna osiągnąć
przez zastosowanie ścianki do rozpraszania energii, studni wypadowej lub szykan
(drastyczne zwiększenie szorstkości dna, rys. 43). W przypadku ścianki (liczonej
jak przelew o kształcie praktycznym) wysokość wody nad jej koroną winna
zapewnić zatopienie odskoku. Podobnie wylot ze studni traktuje się jako przelew
o szerokiej koronie zapewniający zatopienie odskoku. Natomiast hydrauliczne
obliczanie szykan nie jest moŜliwe. Ruch wody ma tu charakter trójwymiarowy
i w praktyce bada się odskok na modelach fizycznych (w skali). W przypadku
ścianki lub studni zapewnienie zatopienia odskoku polega na spełnieniu warunku:
c + H ≥ hs .
Obliczenia inŜynierskie wymagają spełnienia warunku zatopienia odskoku
z odpowiednim marginesem bezpieczeństwa (nadmiarem).
109
7. FILTRACJA
∆H
l v
Rys. 44
Parametr k nosi nazwę w s p ó ł c z y n n i k a f i l t r a c j i .
strop
k
H
strefa saturacji
m
h
spąg
filtracja ze swobodnym
filtracja pod ciśnieniem
zwierciadłem
Rys. 45
Jeśli zwierciadło pojawia się w obszarze warstwy wodonośnej, nosi nazwę
z w i e r c i a d ł a s w o b o d n e g o . Rozgranicza ono s t r e f ę s a t u r a c j i (strefę
pełnego nasycenia wodą, ośrodek dwufazowy) od s t r e f y a e r a c j i
(nienasyconej, ośrodek trójfazowy). Gdy w warstwie wodonośnej ciśnienie wody
jest wyŜsze od atmosferycznego, woda w piezometrze wznosi się ponad strop.
Wystąpienie wody po przewierceniu się przez strop warstwy nosi nazwę
z w i e r c i a d ł a n a w i e r c o n e g o , natomiast ustalony powyŜej stropu poziom
zwierciadła to tzw. z w i e r c i a d ł o n a p i ę t e .
Klasyczny opis filtracji obejmuje wyłącznie strefę saturacji. Przy
zwierciadle napiętym granice pola filtracji są stałe i pokrywają się z granicami
geometrycznymi warstwy. Natomiast zwierciadło swobodne stanowi górną granicę
pola filtracji, a jego połoŜenie zaleŜy od ciśnienia wody. Dla kaŜdego z tych
przypadków uzyskuje się inne rozwiązania i dlatego wyróŜnia się f i l t r a c j ę
pod ciśnieniem i filtrację ze swobodnym zwierciadłem.
Głębokość strumienia filtracji o swobodnym zwierciadle to tzw. m i ą Ŝ s z o ś ć
h y d r o d y n a m i c z n a h . W przypadku filtracji pod ciśnieniem miąŜszość
strumienia jset stała: h = m .
115
H ∆H2
∆H
∆H3
k1 k2 k3 q3 q
q q1 q2 m3
m q4
k4 m4
Rys. 46
Przy filtracji liniowej w warstwie wodonośnej zakłada się identyczność
obrazu ruchu w kolejnych przekrojach pionowych równoległych do kierunku
przepływu. MoŜna zatem przyjąć j e d n o s t k o w ą s z e r o k o ś ć s t r u m i e n i a
f i l t r a c j i . Wynika stąd nowa wielkość opisująca objętość przepływu –
– wydatek jednostkowy:
Q
q= [m2/s].
b
Z definicji prędkości filtracji wynika równanie ruchu:
v A
q= =v h = k hS.
b
W tym przypadku ciągłość przepływu odnosi się do wydatku jednostkowego:
116
q = const .
W warstwie niejednorodnej, której wodoprzepuszczalność zmienia się
skokowo, mogą wystąpić dwa podstawowe schematy przepływu (rys. 46):
szeregowy i równoległy.
W u k ł a d z i e s z e r e g o w y m równanie ciągłości przybiera formę:
∆ H1 ∆ H2
q = q1 = q 2 → q = k1 m = k2 m ,
l1 l2
natomiast suma strat:
∆H = ∆H1 + ∆H2 ;
w układzie równoległym:
q = (k 3 m3 + k 4 m 4 )
∆ H3
q = q3 + q 4 → ,
l3
poniewaŜ:
∆H = ∆H3 = ∆H4 lub S3 = S4 = S = ∆H3 / l3 ;
w układzie mieszanym:
q = q2 = q3 + q4 ,
∆H = ∆H2 + ∆H3 .
W analogiczny sposób sumują się wydatki i straty podczas przepływu
w kolumnie filtracyjnej (rys. 44), przy czym nie musi być ona pionowa. Widać tu
pełną analogię do obliczeń rurociągów szeregowych i równoległych.
l rs
R
Hg Ho Hs
H H
k m Hd k m
x r
ruch liniowy ruch osiowosymetryczny
Rys. 47
117
H = H g − (H g − H d )
x
l
lub H = Hg − Hg − Hd
2
( 2 2
) xl .
l
q( x ) = − w − x ,
2
uwzględniając narastanie wydatku od granicy wododziału. Jednostronny dopływ
do rowu wyniesie :
l
q =w .
2
H Hmax
Hr k Hr
Rys. 48
Opisując średni przepływ w warstwie (q /2) równaniem ruchu dla filtracji
swobodnej uzyskuje się wzór na rozstaw rowów:
l=2
k
w
(2
H max − H r2 )
oraz, analogicznie jak w poprzednim przypadku, równanie swobodnego
zwierciadła:
H = H r2 +
w
k
(
l x − x2 . )
Pominięto tutaj wyprowadzenie dla filtracji pod ciśnieniem, gdyŜ w takich
warunkach nie jest moŜliwe zasilanie przez nieprzepuszczalny nadkład.
8. LITERATURA
8.2. NORMY