Giáo trình bài tập clctt6 b

Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

Chương 5: DÒNG CHẢY ĐỀU TRONG KÊNH HỞÛ
5.1 Khái niệm chung:
Dòng chảy trong kênh hở là dòng chảy có mặt thoáng và không áp (áp suất tại điểm cao nhất

D
AT

EC

H

EN

G

VN

của mặt cắt ướt là áp suất khí trời).

H.5.1a
+ Tính chất dòng chảy đều:
(i) Chiều sâu, diện tích ướt và biểu đồ phân bố vận tốc tại các mặt cắt dọc theo dòng
chảy không đổi.
(ii) Đường năng, mặt thoáng và đáy kênh song song với nhau.

122

Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

H.5.1b
+ Điều kiện xuất hiện dòng đều:
(i) Kênh lăng trụ (mặt cắt kênh ω(h) không đổi dọc theo dòng chảy )

VN

(ii) Độ nhám lòng kênh là hằng số dọc kênh (n=const)

G

(iii) Chiều dài kênh khá dài (L >>1), đủ để dòng chảy đạt đến sự cân bằng giữa tổn thất
năng lượng dọc dòng chảy và sự giảm thế năng do độ dốc kênh tạo ra ⇒

EN

(iv) Độ dốc kênh i > 0.
(v) Không có lưu lượng vào ra dọc dòng
chảy (Q = const dọc theo chiều dòng

H

chảy)

EC

+ Dòng chảy đều rất ít khi xảy ra trong thực

tế. Tuy nhiên nó là khái niệm cơ bản được
dùng trong thiết kế kênh.

D
AT

- Chiều sâu:

h

B

h

1

m

ψ

m

b

H.5.2 Mặt
cắtMat
ứơtcatkênh
hình thang
H.1.2
uot kenh

- Bề rộng đáy kênh:

b

- Bề rộng mặt thoáng:

B

- Diện tích mặt cắt ướt:

A

- Chu vi ướt:

P

- Bán kính thủy lực:

R=

- Hệ số mái dốc:

m = cotg(ψ)

A
P

5.2 Công thức Chezy và Manning:
Công thức bán thực nghiệm hoặc thực nghiệm để tính vận tốc dòng đều:

V = C.Rx.iy

(5.1)
123


Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

Với

V: vận tốc trung bình mặt cắt
R: bán kính thủy lực
i : độ dốc đáy kênh
C: hệ số ma sát
x, y: hằng số

+ Công thức Chezy (1769):
- Xét đoạn dòng chảy có chiều dài L, góc nghiêng đáy kênh θ (rất nhỏ)
- Lực gây chuyển động:

VN

Fm = W.sin(θ) = γ.A.L.sin(θ)
Với:
: diện tích ướt

L

: chiều dài đoạn dòng chảy

θ

: góc nghiêng đáy kênh rất nhỏ ⇒ sin(θ) ≅ tg(θ) ≅ i

EN

G

A

(5.2)

H

- Lực cản trở chuyển động:

EC

Dòng chảy trong kênh thường là dòng rối khu sức cản bình phương (thành hoàn toàn
nhám) ⇒ lực ma sát tỉ lệ với bình phương vận tốc.
⇒

D
AT

Gọi τo là ứng suất ma sát giữa dòng chảy và thành kênh (lực cản trên một đơn vò diện tích)
τo = k.V2

(5.3)

Với k = const
⇒ lực cản:

FR = L.P.k.V2

(5.4)

r
r
- Chuyển động đều ⇒ ∑ F = 0

⇒

Fm - FR = 0 ⇒
γ.A.L.i - L.P.k.V2 = 0

(5.5)

γ.A.L.i = L.P.k.V2

(5.6)

⎛γ ⎞
V =⎜ ⎟
⎝k⎠

1/ 2

R.i

(5.7)
124


Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

Đặt:

⎛γ ⎞
C =⎜ ⎟
⎝k⎠
⇒

k=

1/ 2

(5.8a)

γ

(5.8b)

C2

C được gọi là hệ số Chezy ⇒ công thức Chezy cho dòng đều:
(5.9)

V = C R.i
+ Công thức Manning (1889):

C=

1 1/ 6
R
n

(5.10)

G

R2/3
i
n

(5.11)

EN

V=

VN

Manning dùng phương pháp thực nghiệm để tìm ra công thức sau:

n: hệ số nhám Manning.

1,49 2 / 3

1,49 1 / 6
R
i=
R
R.i
n
n

(5.12)

C=

1,49 1 / 6
R
n

(5.13)

EC

V=

D
AT

⇒

H

Trong hệ thống đơn vò Anh (foot-pound-second) ⇒

τo = k.V2 =

γ

C

2

.V 2 = γRJ

(5.3a)

5.3 Xác đònh hệ số nhám:

+ Trong đường ống kín, hệ số nhám n phụ thuộc:

- Số Reynolds Re
- Độ nhám đường ống (ε)
- Hình dạng mặt cắt ướt (D).
Trong thực tế vì dòng chảy trong kênh có số Reynolds rất lớn ⇒ nên thuộc khu chảy rối
thành nhám hoàn toàn ⇒ hệ số nhám chỉ phụ thuộc vào hình đạng và bề mặt của kênh.
+ Các yếu tố ảnh hưởng đến hệ số nhám:
125


Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

- Độ nhám bề mặt là yếu tố quan trọng ảnh hưởng đến hệ số nhám. Vật liệu đáy kênh

mòn, nhỏ ⇒ n nhỏ & ít thay đổi khi mực nước thay đổi, và ngược lại.

- Lớp phủ thực vật (cỏ dại, cây nhỏ) có thể làm gia tăng hệ số nhám kênh thiên nhiên.
- Hình dạng mặt cắt kênh ảnh hưởng đến giá trò hệ số nhám. Chẳng hạn đối với kênh
thiên nhiên, ở những nơi kênh bò bồi thành bậc ⇒ n gia tăng.

- Vật cản như cây gỗ làm gia tăng hệ số nhám n.
- Tuyến kênh uốn cong với bán kính cong nhỏ làm gia tăng hệ số nhám n.
- Sự bồi xói làm hệ số nhám thay đổi.
- Mực nước và lưu lượng ảnh hưởng đến độ sâu dòng chảy. Do đó có thể ảnh hưởng đến

VN

hệ số nhám:

) Kênh không bò bao phủ: n giảm khi mực nước hoặc lưu lượng tăng.

G

) Kênh bò bao phủ bởi lớp thực vật ở mái dốc: n tăng khi mực nước hoặc lưu lượng

EN

tăng.

5.3.1 Trường hợp mặt cắt kênh đơn giản:

H

+ Phương pháp SCS (Soil Conversation Service method):

vật liệu ⇒ n1 .

EC

- Ước đònh hệ số nhám n cơ bản cho kênh thẳng, mặt cắt lăng trụ, đáy trơn và chỉ có một

D
AT

- Cộng hoặc nhân với hệ số hiệu chỉnh tùy điều kiện thực tế.
Ví dụ:

n1 = 0,02

(kênh đất)

n2 = 0,005

(phủ cỏ)

n3 = 0,01

(kênh có tiết diện thay đổi không đều, loại nhỏ)

n = n1 + n2 + n3 = 0,035
+ Phương pháp dùng bảng:
Dựa vào thực nghiệm và kinh nghiệm, n được lập thành bảng. Khi sử dụng, n được tra dựa
vào bảng lập sẵn.
+ Phương pháp hình ảnh:
Đo đạc xác đònh n ⇒ chụp ảnh và xếp loại. Khi sử dụng, chọn n thông qua hình ảnh.

+ Phương pháp dùng biểu đồ lưu tốc:

126


Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

Dựa vào phương trình phân bố lưu tốc:

n=

( x − 1)h 1 / 6
6,78.( x + 0.95)

(5.14)

Với
h: chiều sâu dòng chảy.
x=

U 0.2
U 0.8

Với U0.2: vận tốc ở vò trí 2/10 của chiều sâu

VN

+ Phương pháp dùng công thức thực nghiệm:

n = 0,047d 1 / 6

EN

Với d: đường kính cở hạt của lòng kênh (mm)

(5.15a)

G

a) Simon và Sentruk (1976):

b) Raudkivi (1976):

EC

c) Meyer and Peter (1948):

(5.15b)

H

1/ 6
n = 0,013d 65

p dụng đối với kênh có vật liệu đáy hỗn hợp, hạt có kích thước nhỏ:
(5.15c)

D
AT

1/ 6
n = 0,038d 90

Với d65 là đường kính cở hạt mà trọng lượng của các hạt có đường kính nhỏ hơn hay bằng nó
chiếm 65% trọng lượng toàn thể mẫu.
5.3.2 Trường hợp mặt cắt kênh phức tạp:
(i) Công thức Horton, Einstein và Blank (vận tốc trung bình bằng nhau):
⎡ N
3/ 2 ⎤
⎢ ∑ Pi .n i ⎥
⎥
ne = ⎢ i = 1
⎥
⎢
P
⎥
⎢
⎦
⎣

2/3

(5.16)

Với: ne : hệ số nhám tương đương.

n1

n2

n3

H.1.3
H.5.3

Pi : chu vi ướt của từng diện tích đơn giản
P : chu vi ướt của toàn bộ mặt cắt

127


Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

N : số lượng các mặt cắt đơn giản
(ii) Lực ma sát trên toàn mặt cắt bằng tổng lực ma sát trên từng mặt cắt:

⎡ N
2 ⎤
⎢ ∑ Pi .n i ⎥
⎥
n e = ⎢ i =1
⎢ P ⎥
⎢
⎥
⎣
⎦

1/ 2

(5.17)

(iii) Lưu lượng trên mặt cắt tổng bằng tổng lưu lượng các mặt cắt đơn giản:

P.R 5 / 3
N
Pi .Ri5 / 3
∑
ni
i =1

(5.18)

VN

ne =

tổng.

EN

(iii) Từ thí nghiệm COX (1973) đề nghò:
N

∑ n .A
i

i =1

i

A

(5.19)

H

ne =

G

Với Ri bán kính thủy lực của từng mặt cắt đơn giản; R bán kính thủy lực của mặt cắt

EC

Ai: diện tích mặt cắt ướt của từng mặt cắt đơn giản.
A : diện tích mặt cắt ướt tổng.

D
AT

5.4 Tính toán dòng đều:

Bài toán về dòng đều có liên quan đến năm yếu tố: Lưu lượng (Q), mặt cắt kênh (b, m), độ
nhám (n), độ dốc (i), độ sâu dòng đều (ho). Công thức sử dụng:

V=
hay Q =

K=

1 2/3
R
i
n

(5.20)

1
A.R 2 / 3 i
n

(5.21)

1
A.R 2 / 3
n

(5.22a)

Q=K i

(5.22b)

K: modul lưu lượng (m3/s), phụ thuộc tính chất và kích thước kênh.
5.4.1 Xác đònh độ sâu dòng đều:
p dụng đối với kênh có mặt cắt phức tạp:

128

Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

+ Phương pháp thử dần:
Thay đổi giá trò của h sao cho phương trình sau đây thỏa mãn:

A.R 2 / 3 =

n.Q
i

(5.23)

= f ( h)

+ Phương pháp đồ thò:

- Đặt y = A.R 2 / 3 = f (h) , vẽ hàm y = f(h).
- Tính

n.Q
i

và đặt lên trục y ⇒ ho.

VN

+ Phương pháp biểu đồ:

- Tùy hình dạng mặt cắt kênh, người ta lập sẵn biểu đồ: K = f(ho).

G

- Biết K, dựa vào biểu đồ suy ra ho.

D
AT

EC

H

EN

- Ví dụ:

H.5.4
Đối với mặt cắt hình tròn có đường kính lớn là D. Người ta lập sẵn biểu đồ các quan hệ sau
đây:

129


Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

K
h

= f1 ( )
K ng
D

W
h
= f2 ( )
W ng
D

hay

Với

K=

Q
1
A.R 2 / 3 =
n
i

(5.24)

1
Ang .R ng2 / 3
n

(5.25)

K ng =

và modul vận tốc:

1 2/3 V
R
=
n
i

W=

1 2/3
R ng
n

(5.26)

VN

W=

G

Với Ang và Rng là diện tích ướt và bán kính thủy lực khi chảy ngập đầy ống h = D.

hoặc

W
W ng

D
AT

+ Phương pháp số:

⎛h⎞
⎛h⎞
⇒ ⎜ ⎟
⇒ h =⎜ ⎟
.D
⎝ D ⎠Trabang
⎝ D ⎠Trabang

EC

K
K ng

H

⇒ Nếu biết:

EN

D
D2
Ang = π
; Rng =
4

4

Sử dụng máy tính thông qua phương pháp tính đúng dần để xác đònh ho của mặt cắt bất kỳ.
5.4.2 Thiết kế kênh:

+ Những yếu tố ảnh hưởng đến lưu lượng chảy:

- Diện tích ướt (A)
- Bán kính thủy lực (R).
- Độ dốc (i).
- Hệ số nhám của kênh (n).
+ Mặt cắt có lợi nhất về mặt thủy lực:

- Mặt cắt có một diện tích ướt nhất đònh nhưng cho lưu lượng lớn nhất ⇒ hình tròn là m/c
có lợi nhất về thủy lực.
130


Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

- Mặt cắt có lợi nhất về mặt thủy lực chưa hẳn là mặt cắt có lợi nhất về mặt kinh tế.
+ Mặt cắt lợi nhất về mặt thủy lực đối với kênh hình thang:
Xét kênh hình thang có đáy b, chiều sâu h và mái dốc m. Với lưu lượng Q cho trước tìm
quan hệ giữa b, m và h sao cho đạt được mặt cắt có lợi nhất về mặt thủy lực ⇒ mặt cắt ướt
và chu vi ướt nhỏ nhất.
A = (b+m.h).h = (β+m).h2
Với:

β=

b
h

P = b + 2h. 1 + m 2 = h.( β + 2 1 + m 2 )

2h 2 .( β + m)

EC

thế vào (5.28) ⇒

(β + 2 1 + m 2 )

=0

D
AT

h2 -

(5.29)

H

dh
h
=dβ
(β + 2 1 + m 2 )

⇒

G

dP dh
=
(β + 2 1 + m 2 ) + h = 0
dβ d β

(5.28)

EN

và

VN

dh
dA
= h2 + 2.h.(β+m).
=0
dβ
dβ

⇒

⇒

β + 2 1 + m 2 = 2 (β+m)

⇒

βln = 2 ( 1 + m 2 -m)

⇒

Rln =

⇒

Rln =

(5.30)

( β ln + m).h
A
=
P ( β ln + 2 1 + m 2 )
(2 1 + m 2 − m).h
2(2 1 + m 2 − m)

=

h
=>
2

(5.31)

Rln =

h
2

(5.32)

Phương trình (5.30) và (5.32) là điều kiện để cho mặt cắt hình thang có lợi nhất về mặt
thủy lực.
+ Các dạng bài toán mặt cắt lợi nhất về mặt thủy lực đối với kênh hình thang:

131


Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

a) Giữ A = const:
1/ 2

⎛ A ⎞
⎟⎟
hln = ⎜⎜
⎝ β ln + m ⎠

và

bln = βln.hln.

(5.32a)

và

bln = βln.hln.

(5.32b)

b) Giữ Q = const:
⎛ 22 / 3.n.Q
hln = ⎜⎜
⎝ (β ln + m ).

⎞
⎟
i ⎟⎠

3/8

+ Mặt cắt lợi nhất về mặt thủy lực của hình chữ nhật:

b
=2
h

(5.33)

VN

Với m/c chữ nhật ⇒ m = 0 ⇒ β = 2 ⇒

⇒ chiều sâu bằng nửa chiều rộng ⇒ đào kênh quá sâu ⇒ không cho phép đối với kênh

G

lớn. Nếu thiết kế với R < Rln vài % ⇒ h giảm đáng kể và hình dạng kênh có thể chấp
nhận được về mặt thực tế.

EN

+ Kiểm soát vận tốc:

Wmax
0,065.i 1 / 4

EC

Vkl =

H

Vkl < V < Vkx

(5.34)

(5.35)

D
AT

Vkl : vận tốc không lắng.

Wmax: tốc độ lắng chìm của hạt có kích thước lớn nhất.
Vkx : vận tốc không xói, tra bảng lập sẵn.

Vkx = Kx.Q 0,1
Trong đó :

Kx : Hệ số phụ thuộc vào đất lòng kênh, xác định theo tài liệu;
Q : Lưu lượng của kênh, m3/s.

132


EN

G

VN

Trường Đại Học Bách Khoa – ĐHQG TP. HCM
PGS. TS. Lê Văn Dực

D
AT

EC

H

H.5.5

133