Dòng ổn định không đều là dòng chảy trong đó vận tốc không đổi theo thời gian
(?u/?t = 0), nhưng không phải là hằng số trong không gian (?u/?x ?0, ?u/?y ?0). Tính
không đều có thể do những thay đổi không gian của mặt cắt ngang hoặc bởi những vật
chắn (đập tràn) trong dòng chảy.
Hai loại dòng không đều được xem xét: dòng chảy biến đổi dần (chậm) và dòng chảy
biến đổi gấp (nhanh).
46 trang |
Chia sẻ: lylyngoc | Lượt xem: 1607 | Lượt tải: 1
Bạn đang xem trước 20 trang tài liệu Chương 7. Dòng chảy ổn định không đều, để xem tài liệu hoàn chỉnh bạn click vào nút DOWNLOAD ở trên
95
Chương 7. dòng chảy ổn định không đều
7.1. Mở đầu
Dòng ổn định không đều là dòng chảy trong đó vận tốc không đổi theo thời gian
(u /t = 0), nhưng không phải là hằng số trong không gian (u /x 0, u /y 0). Tính
không đều có thể do những thay đổi không gian của mặt cắt ngang hoặc bởi những vật
chắn (đập tràn) trong dòng chảy.
Hai loại dòng không đều được xem xét: dòng chảy biến đổi dần (chậm) và dòng chảy
biến đổi gấp (nhanh).
Trong dòng biến đổi dần, vận tốc thay đổi dần dần từ mặt cắt này đến mặt cắt
khác; những đường dòng thực chất là song song và có thể giả thiết áp suất thủy tĩnh.
Giả thiết khác là tổn thất năng lượng cũng được xem như đối với dòng chảy đều.
Trong dòng biến đổi gấp, những thay đổi độ sâu, bề rộng và do đó biến đổi vận tốc
xảy ra trong những đoạn lòng dẫn ngắn; những đường dòng uốn cong mạnh và áp suất
trong dòng chảy không phải là thủy tĩnh (ví dụ: nước nhảy thuỷ lực, chảy tràn tự do,
dòng chảy qua đập tràn). Thông thường, tổn thất năng lượng do ma sát đáy có thể bỏ
qua so với những tổn thất năng lượng khác (tổn thất do mở rộng).
Những dòng không đều cũng có thể nghiên cứu bằng việc giả thiết dòng thế, có
nghĩa là bỏ qua những ứng suất trượt do nhớt (độ nhớt bằng không).
Các lực dòng chảy tác động lên công trình do những biến đổi vận tốc và áp suất
cũng được xem xét trong chương này.
7.2. Dòng thế
7.2.1. Mở đầu
Dòng thế là dòng chảy không quay với những thành phần vận tốc có thể dẫn xuất
từ một hàm thế (). Vì dòng chảy không quay (có nghĩa là xoáy bằng không, xem mục
4.6), phương trình Bernoulli hợp lệ cho toàn bộ trường dòng chảy (xem mục 5.4.3).
Từ bức tranh đường dòng và phương trình Bernoulli, có thể nhận được những biến
đổi vận tốc và áp suất trong trường dòng chảy. Những đường dòng cùng với đường thế
liên quan với chúng hình thành một lưới dòng (lưới thuỷ động lực) có thể xây dựng dễ
96
dàng bằng phép thử sai.
Nhiều chất lỏng thực có thể biểu thị như những dòng thế khi hiệu ứng ma sát nội
là nhỏ để bỏ qua, thông thường là trường hợp cho dòng chảy tăng tốc (gấp). Ví dụ, hình
7.3 cho thấy một lưới dòng trong khu vực tăng tốc của ống dẫn. ứng dụng cách tiếp cận
lưới dòng không cho ta những kết quả chính xác đối với những dòng chảy giảm tốc
nhanh.
7.2.2. Dòng thế hai chiều
Dòng thế có thể mô tả dưới dạng thế vận tốc () và hàm dòng (). ứng dụng của
hàm thế () và hàm dòng () chỉ thích hợp cho một trường dòng chảy hai chiều: u =
f(x,y) hoặc u = f(x,z).
Đối với trường dòng chảy hai chiều thẳng đứng không quay, có thể xác định thế
vận tốc (), sao cho những thành phần vận tốc u và w là:
x
u
và
z
w
. (7.2.1)
Phương trình liên tục như sau:
0
z
w
x
u
(7.2.2)
hoặc
0
2
2
2
2
zx
. (7.2.3)
Phương trình (7.2.3) được gọi là phương trình Laplace.
Đối với một hàm liên tục nó dẫn đến:
xzzx
22
hoặc (7.2.4)
x
w
z
u
(7.2.5)
có nghĩa là dòng chảy không quay (xem mục 4.6) và cho thấy rằng những thành phần
vận tốc đương nhiên có thể thể hiện như những gradient của đại lượng vô hướng
(phương trình 7.2.1).
Trong mục 4.3 đã chỉ ra rằng những thành phần vận tốc (u, w) cũng liên quan đến
hàm dòng (), như sau:
z
u
và
x
w
. (7.2.6)
Một đường mà dọc theo đó không đổi biểu thị một đường dòng. Những giá trị
khác nhau thể hiện những đường dòng khác nhau.
Thay phương trình (7.2.6) vào phương trình (7.2.2) dẫn đến:
97
0
22
xzzx
(7.2.7)
nói rằng những thành phần vận tốc cũng có thể biểu thị như gradient của một đại
lượng vô hướng (phương trình 7.2.6).
Thế vận tốc và hàm dòng đều là những hàm số của cả x lẫn z. Như vậy,
wdzudxdz
z
dx
x
d
(7.2.8)
udzwdxdz
z
dx
x
d
. (7.2.9)
Thế không đổi dọc theo một đường đẳng thế, có nghĩa là d = 0.
Hàm dòng không đổi dọc theo một đường dòng, có nghĩa là d = 0.
Hai họ đường có thể xây dựng trong mặt phẳng x - z:
đường đẳng thế: udx + wdz = 0 (7.2.10)
đường dòng: wdx - udz = 0. (7.2.11)
Những phương trình (7.2.10) và (7.2.11) thể hiện một hệ trực giao của nhiều đường,
mà có nghĩa là taị bất kỳ giao điểm nào những đường dòng đều thẳng góc với những
đường đẳng thế (xem hình 7.1). Những đường đẳng thế cũng thẳng góc với những biên
của khu vực dòng chảy, bởi vì những biên cũng là những đường dòng.
Hình 7.1. Đường dòng thẳng góc với đường đẳng thế
7.2.3. Lưới dòng (lưới thuỷ động lực)
Lưới dòng là một họ các đường dòng và đường đẳng thế tạo nên các hình vuông
cong. Trong một lưới dòng khoảng cách giữa các đường dòng () và khoảng cách giữa
các thế vận tốc () bằng nhau trong toàn bộ lưới. Điều này có thể thấy bằng việc áp
dụng một hệ toạ độ tự nhiên (s, n), xem hình 7.1. Các vận tốc trong hệ là:
98
ns
vs
(7.2.12)
0
sn
vn
. (7.2.13)
Từ phương trình (7.2.13) thấy rằng:
n
s
. (7.2.14)
Bằng việc lấy hình vuông (s = n), dẫn đến = .
Từ phương trình (7.2.12) dẫn đến:
vsn = = q = const. (7.2.15)
Như vậy, lưu lượng giữa 2 đường dòng bằng nhau và không đổi.
Hình 7.2. Lưới dòng (lưới thuỷ động lực)
Ví dụ trong hình 7.2, lưu lượng q giữa 2 đường dòng là q = 3 m2/s.
Khi biết lưu lượng toàn bộ q, lưu lượng q có thể rút ra từ q = q/m, trong đó m là
số lượng các khoảng đường dòng giữa đáy và mặt nước. Bằng cách đo những giá trị n,
có thể xác định vận tốc tại bất kỳ điểm nào:
vs1n1 = vs2n2 = q. (7.2.16)
Trường áp suất có thể xác định từ phương trình Bernoulli:
2
2
2
2
1
1
2
1
22
z
g
p
g
v
z
g
p
g
v ss
. (7.2.17)
Với nhiều mục đích những lưới dòng có thể vẽ tay. Những đường dòng được vẽ cách
đều nhau một khoảng nào đó nơi dòng chảy song song. Số lượng đường dòng phụ thuộc
vào độ chính xác mong muốn. Khoảng cách càng nhỏ, độ chính xác càng cao và đòi hỏi
nỗ lực lớn hơn khi vẽ lưới dòng. Sau khi những đường dòng đã được vẽ bằng mắt, vẽ
những đường đẳng thế. Thực hiện những điều chỉnh liên tiếp cho cả những đường dòng
lẫn những đường đẳng thế, cho đến khi những hình vuông thích hợp xuất hiện. Để
kiểm tra, có thể vẽ những đường chéo qua đỉnh các hình vuông, chúng cũng phải lập
99
thành một hệ trực giao.
7.2.4. ứng dụng
Dòng thế chảy qua một ống dẫn
Hình 7.3. Dòng thế qua một ống (Thijsse, 1951)
Hình 7.3 cho thấy lưới dòng đối với dòng chất lỏng chảy từ một hồ chứa có mực
nước cao thông qua một ống dẫn đến một hồ chứa có mực nước thấp. Chênh lệch cột
nước là h = 2 m.
áp dụng phương trình Bernoulli từ điểm 1 đến điểm 2 dẫn đến:
2251575100
22
21
21
2
1
2
2 ,,zz
g
p
g
p
g
u
g
u
u2 = 6,26 m/s.
Lưu lượng (q2) giữa những đường dòng dưới điểm 2 là q2 = u2n2. Giá trị n2 có
thể xác định bằng cách đo từ lưới dòng, kết quả là n2 = 0,22 m và do đó q2 = 6,26 x
0,22 = 1,41 m3/s. Lưu lượng toàn bộ q = 5q2 = 7,05 m
2/s.
Bây giờ có thể xác định vận tốc và áp suất tại những điểm khác bằng việc sử dụng
phương trình liên tục và phương trình Bernoulli.
Những kết quả cho trong bảng sau:
100
Điểm z
(m)
n
(m)
u
(m/s)
u2/2g
(m)
p/g
(m)
z+p/g
(m)
1 5,00 0 0 0 5,00
2 1,25 0,22 6,26 2,00 1,75 3,00
3 0 0,32 4,29 0,92 4,08 4,08
4 1,15 0,11 12,52 8,00 -4,15 -3,00
5 0 0 0 5,00 5,00
6 3,00 - 0 0 0 3,00
Vận tốc lớn nhất xuất hiện tại điểm nơi hình vuông nhỏ nhất, gần điểm 4 (u4 = 12,5
m/s. áp suất tại điểm 4 là số âm, có nghĩa là thấp hơn áp suất không khí. Những ví dụ
khác của những lưới dòng cho trong hình 7.4.
Hình 7.4. Ví dụ của lưới dòng
7.3. Dòng chảy rối biến đổi dần
7.3.1. Mở đầu
Trong trường hợp của dòng biến đổi dần (ổn định) độ sâu nước thay đổi chậm theo
chiều dài lòng dẫn. Những đường dòng thực tế là song song nên áp suất chất lỏng là
thủy tĩnh. Giả thiết cơ bản cho những loại dòng chảy này là: có thể xác định ứng suất
trượt tại đáy cho mỗi đoạn bằng việc áp dụng công thức sức cản cho dòng đều. Như vậy,
2
2
C
u
gb (7.3.1)
cũng hợp lệ cục bộ đối với dòng biến đổi dần. Điều này có nghĩa là tổn thất cột nước
trong một đoạn cũng như đối với một dòng chảy đều.
Hai loại đường cong mặt nước có thể phân biệt trong dòng không đều:
1. đường nước dâng khi độ sâu dòng chảy tăng theo hướng dòng chảy (dh/dx > 0),
101
và
2. đường nước hạ khi độ sâu dòng chảy giảm theo hướng dòng chảy (dh/dx < 0).
Một đường nước dâng phát sinh khi dòng chảy được ngăn bởi một đập tràn (hình
7.7.16). Một đường nước hạ phát sinh trong trường hợp chảy tràn tự do (hình 7.7.16).
7.3.2. Phương trình Belanger
Những đường cong nước dâng và nước hạ có thể xác định bằng việc áp dụng
phương trình Belanger, rút ra từ phương trình (5.4.54), cho một mặt cắt ngang tuỳ ý
hoặc bằng phương trình (5.4.55) cho một mặt cắt ngang hình chữ nhật rộng (b >> h).
Phương trình Belanger cũng có thể dẫn xuất bằng việc áp dụng cân bằng động
lượng đối với một phần tử chất lỏng có chiều dài x và chiều cao h như trong hình 7.5
đối với lòng dẫn chữ nhật rộng.
Lực áp suất thực tế trên đơn vị bề rộng (Fp) theo hướng s là:
s
ds
dh
ghFp . (7.3.2)
Các lực khác theo hướng s là trọng lực FG,H và lực ma sát đáy Fw:
F G,H = gh s sin (7.3.3)
Fw = -bs. (7.3.4)
Hình 7.5. Các lực trong dòng không đều
Gia tốc của phần tử chất lỏng theo hướng s là:
ds
ud
u . (7.3.5)
Phương trình chuyển động dẫn đến:
Fs = mas
102
))((sin
ds
ud
ushssgh
ds
dh
sgh b (7.3.6)
hoặc
0sin
h
g
ds
dh
g
ds
ud
u b
. (7.3.7)
Trong trường hợp độ dốc nhỏ sin = ib và phương trình (7.3.7) trong hệ thống toạ độ
x – z là:
0
h
gi
ds
dh
g
ds
ud
u bb
. (7.3.8)
Phương trình liên tục dẫn đến:
0
)(
dx
ud
h
dx
dh
u
dx
hud
. (7.3.9)
Thay phương trình (7.3.9) vào phương trình (7.3.8) cho ta:
h
gi
ds
dh
h
u
g bb
)(
2
đối với ib 0. (7.3.10)
Cuối cùng thay phương trình (7.3.1) dẫn đến phương trình Belanger:
b
b i
gh
u
ihC
q
dx
dh
2
32
2
1
1
. đối với ib > 0 (7.3.11)
Phương trình (7.3.11) cũng có thể biểu thị như sau:
b
c
e i
hh
hh
dx
dh
33
33
đối với ib > 0 (7.3.12)
trong đó:
3/2)(
b
e
iC
q
h = độ sâu cân bằng tại cùng lưu lượng q ứng với công thức Chezy
3/1
2
)(
g
q
hc = độ sâu phân giới tại cùng q.
Đối với những giá trị đã cho của q, ib, và C, biết được những độ sâu he và hc, ta có
phương trình vi phân đối với độ sâu h. Nếu cần thiết, hệ số cũng được xét đến bằng
việc áp dụng hc = (q
2/g)1/3.
Khi đã biết giá trị ks, độ sâu cân bằng he có thể tính toán từ:
q = Che(heib)
1/2 = [18log(12he/ks)]he
3/2 ib
1/2 .
7.3.3. Phân loại những đường cong mặt nước
103
Hình 7.6. Phân loại những đường cong mặt nước (De Vries, 1985)
104
Phương trình (7.3.12) có thể giải được khi biết q, C, ib và những điều kiện biên ở x =
0. Trong trường hợp dòng dưói phân giới (Fr < 1) cần biết độ sâu nước hạ lưu. Trong
trường hợp dòng trên phân giới cần biết độ sâu nước thượng lưu, bởi vì sóng mặt không
thể lan truyền ngược hướng dòng chảy trong trường hợp này. Nhiều lời giải của phương
trình (7.3.12) đã cho trong hình 7.6. Phân biệt những loại sau:
A – Loại đường cong cho độ dốc đáy ngược
H – Loại đường cong cho đáy nằm ngang
M – Loại đường cong cho độ dốc đáy thuận vừa phải (he > hc)
C – Loại đường cong cho dốc đáy thuận tới hạn (he = hc)
S – Loại đường cong cho độ dốc đáy thuận rất dốc (he < hc).
Độ dốc phân giới xuất hiện khi he = hc hoặc ic = g/C
2, có thể dẫn xuất từ những biểu
thức đối với he và hc (phương trình 7.3.12). Dốc vừa phải là ib ic.
Nhận xét
1. Đường cong mặt nước được biết nhiều nhất là đường nước dâng thượng lưu một
đập tràn trong trường hợp độ dốc đáy vừa phải (ib = 10
-4). Nó phát sinh đường cong M1
vì h > he > hc. Độ sâu nước giảm theo hướng thượng lưu (xem hình 7.6).
2. Trong trường hợp đáy ngang độ sâu cân bằng he lớn vô hạn (he = ). Như vậy
không thể h > he. Như vậy, đường cong loại H1 không tồn tại.
3. Trong trường hợp he = hc thấy rằng dh/dx = ib dẫn đến mặt nước nằm ngang đối
với h > he = hc và h < he = hc.
Các ví dụ:
Vấn đề trong tính toán những đường cong mặt nước là tìm được độ sâu nước thích
hợp, hoặc hạ lưu trong trường hợp dòng dưói phân giới (Fr < 1) hoặc thượng lưu trong
trường hợp dòng trên phân giới (Fr > 1).
1. Lối ra tự do hoặc chảy tràn tự do
Trong trường hợp đáy dốc vừa phải, có nghĩa là ib < ic với ic = g/C
2 và độ sâu nước h0
tại đầu cuối lòng dẫn bằng độ sâu phân giới, như vậy là h0 = hc (xem hình 7.7.1). Đó là
quy luật tự nhiên mà mặt nước tìm được vị trí thấp nhất có thể của nó ứng với độ cao
năng lượng nhỏ nhất Hemin, như đã cho trong mục 6.6.2. Vì dòng chảy là phân giới, cần
sử dụng độ sâu hạ lưu h0 = hc = (q
2/g)1/3 làm điều kiện biên. Đường mặt nước là một
đường cong loại M2. Nếu mực nước hA trong thuỷ vực tại hạ lưu lòng dẫn lớn hơn hc nh-
ưng nhỏ hơn he (trạng thái 2), sẽ nhận được đường cong M2 khác. Nếu mực nước hA lớn
hơn he (trạng thái 3), sẽ được nhận được đường cong M1.
Trong trường hợp đáy rất dốc, có nghĩa là ib > ic và hc > he, độ sâu nước tại đầu cuối
lòng dẫn sẽ bằng độ sâu cân bằng he = (q/Cib
0,5)2/3 khi lòng dẫn đủ dài. Vì dòng chảy trên
phân giới (Fr > 1), độ sâu nước trong lòng dẫn phụ thuộc vào điều kiện biên thượng lưu
105
(xem hình 7.7.2). Nếu lòng dẫn ngắn, thì có thể không đạt đến độ sâu cân bằng, và độ
sâu h0 có thể lớn hơn hoặc nhỏ hơn he phụ thuộc vào độ sâu nước thượng lưu.
Hình 7.7.2. Lối ra tự do trong trường hợp rất dốc (Fr > 1)
Nếu độ sâu hA lớn hơn hc, nước nhảy thủy lực sẽ phát sinh ở giao điểm của đường
cong S1 và độ sâu h2 (phương trình 7.4.11).
2. Lối vào tự do
Trong trường hợp độ dốc vừa phải có thể tính toán độ sâu h0 ở lối vào của lòng dẫn
bằng việc áp dụng phương trình Bernoulli (xem hình 7.7.3):
g
u
hh A
2
2
0
0 .
Vì dòng chảy dưới phân giới, sẽ hình thành một bậc nước.
Trong trường hợp rất dốc độ sâu h0 ở lối vào bằng độ sâu phân giới hc. Vì dòng chảy
trên phân giới (Fr > 1), điều kiện biên áp dụng ở thượng lưu phát sinh đường cong S2
(hình 7.7. 4).
Hình 7.7.1. Lối ra tự do trong trường hợp độ dốc vừa phải (Fr < 1)
106
Hình 7.7.3. Lối vào tự do trong trường hợp độ dốc vừa phải (Fr < 1)
Hình 7.7.4. Lối vào tự do trong trường hợp rất dốc (Fr > 1)
3. Đập tràn chảy trên mặt
Hình 7.7.5. Đập tràn trong trường hợp độ dốc vừa phải (Fr < 1)
Thông thường, những đập tràn được sử dụng để làm tăng độ sâu nước thượng lưu
đập cho mục đích giao thông thuỷ. Hình 7.7.5 cho thấy một đập tràn tràn trong trường
hợp độ dốc vừa phải (ib < ic). Độ sâu nước hạ lưu đập tràn sẽ bằng độ sâu cân bằng. Độ
sâu nước ở trên đập tràn sẽ xấp xỉ bằng độ sâu phân giới hc. Độ sâu nước ngay tại
thượng lưu đập tràn sẽ bằng h0 = d + H, trong đó d là chiều cao đập tràn và H là mực
nước thượng lưu đập tràn. Mực nước H dẫn xuất từ công thức lưu lượng cho đập tràn
(xem mục 7.4.7) là H = q2/3 với là hệ số liên quan đến loại đập tràn và q là lưu lượng
107
trên đơn vị bề rộng cần biết.
Vì dòng chảy dưới phân giới, điều kiện biên để tính toán mặt nước cho tại hạ lưu
h0 = d + H, đường mặt nước là đường cong M1.
4. Đập chảy dưới sâu
Có thể phân biệt hai trạng thái: dòng dưới phân giới (hình 7.7.6) hoặc dòng trên
phân giới (hình 7.7.7) ở hạ lưu đập tràn trong trường hợp độ dốc vừa phải.
Hình 7.7.6. Đập chảy dưới sâu (dốc vừa phải)
Hình 7.7.6 cho thấy dòng dưới phân giới ở hạ lưu đập tràn. Những độ sâu nước hạ
lưu và thượng lưu đập tràn liên quan với nhau qua công thức tính lưu lượng như sau:
q = a[2g(h0 - he)]
0,5
trong đó là hệ số liên quan đến loại đập tràn và a là chiều cao mở (xem mục 7.4.7). Độ
sâu nước h0 có thể tính toán khi biết q, , a và he. Vì dòng chảy dưới phân giới, có thể
tính toán đường cong mặt nước M1 bằng việc áp dụng h0 làm điều kiện biên hạ lưu.
Hình 7.7.7. Đập chảy dưới sâu (dốc vừa phải)
Hình 7.7.7 cho thấy dòng trên phân giới hạ lưu đập tràn. áp dụng phương trình
Bernoulli và phương trình cân bằng động lượng, lưu lượng có thể biểu thị như sau:
0
0
2
ha
g
ahq
108
trong đó a là chiều cao mở và là hệ số co hẹp (h1 = a, = 0,6). Độ sâu nước h0 có thể
tính toán khi q, và a được biết.
Mặt nước giữa đập tràn và nước nhảy thủy lực là đường cong M3. Vị trí của nước
nhảy thủy lực đã cho trong mục 7.4.3.
5. Sự quá độ từ đáy trơn đến đáy nhám
a) dốc vừa phải: ib < ic và hc < he và Fr < 1
Hình 7.7.8.
b) rất dốc: ib > ic và hc > he và Fr > 1
Hình 7.7.9.
6. Sự quá độ độ dốc đáy
a) dốc vừa phải: ib < ic và hc < he và Fr < 1
Hình 7.7.10.
b) rất dốc: ib > ic và hc > he và Fr > 1
109
Hình 7.7.11.
7. Sự quá độ từ dốc vừa phải đến rất dốc
Hình 7.7.12.
8. Sự quá độ từ rất dốc đến dốc vừa phải
Nước nhảy sẽ phát sinh trong lòng dẫn rất dốc khi độ sâu he2 lớn hơn độ sâu h2 của
nước nhảy (phương trình 7.1.1.1). Nếu không, nước nhảy sẽ phát sinh trong lòng dẫn
dốc vừa phải.
Hình 7.7.13.
110
9. Sự quá độ trong bề rộng dòng chảy
a) dốc vừa phải: ib < ic và hc < he và Fr < 1
b) rất dốc: ib > ic và hc > he và Fr > 1
Nếu đoạn 2 ngắn, thì đường cong S2 không thể tự nó phát triển trong đoạn đó; thay
vào đó phát sinh một bậc nước.
10. Kết hợp những đường cong mặt nước
Hình 7.7.16 cho thấy các kết hợp của những đường cong mặt nước cho độ dốc đáy
vừa phải và rất dốc.
Hình 7.7.14.
Hình 7.7.15.
111
Hình 7.7.16. (Henderson, 1970)
7.3.4. Tính toán giải tích những đường cong mặt nước
Có thể thực hiện tính toán giải tích cho những trạng thái đặc biệt, như dòng chảy
trên một đáy nằm ngang (ib = 0). Phương trình (7.3.12) có thể biểu thị như sau:
33
2
2
3
c
b
hh
C
q
ih
dx
dh
. (7.3.13)
Thay ib = 0 dẫn đến:
33
2
2
chh
C
q
dx
dh
hoặc dhhhdx
C
q
c )(
33
2
2
(7.3.14)
)(
000
33
2
2
dhhdhh
q
C
xd
h
h
c
h
h
x
x
(7.3.15)
)](
4
1
)([ 40
4
101
3
2
2
01 hhhhh
q
C
xx c (7.3.16)
)](
4
1
)([ 40
4
101
3
2
2
01 hhhhh
q
C
xx c (7.3.17)
trong đó:
x0 = vị trí nơi điều kiện biên h = h0 được biết (hình 7.8)
112
x1 = vị trí nơi độ sâu nước h bằng h1 (hình 7.8).
Hệ số Chezy giả thiết không đổi, nhưng cũng có thể tính theo C = 18log(12 h/ks) với
h =1/2(h0 + h1) = độ sâu nước trung bình trên khoảng cách x1 – x0.
Hình 7.8. Tính toán giải tích đường cong mặt nước
Những tính toán giải tích cho nhiều trạng thái phức tạp hơn có thể thực hiện bằng
phương pháp của Bresse.
Phương pháp của Bresse
Độ sâu nước được biểu thị như sau: h = he.
Phương trình Belanger biểu thị như sau:
b
e
c
e
e i
h
h
h
h
h
h
dx
dh
3
3
3
3
3
3
1
(7.3.18)
b
b
e i
g
iCdx
d
h
2
3
3 1
d
g
iC
i
h
dx
b
b
e
13
2
3
d
g
iC
i
h
dx
b
b
e
1
11
3
2
3
113
d
g
iC
i
h
dx
b
b
e ]
1
1
1[
3
2
]
1
1
)1([
3
2
d
g
iC
d
i
h
dx b
b
e
. ( 7.3.19)
Tích phân cho ta ( = 1 – C2ib/g):
]
1
)[(
1
0
30101
d
i
h
xx
b
e (7.3.20)
))]()(()[( 010101
b
e
i
h
xx . (7.3.21)
Những giá trị do Bresse (xem Bảng 7.1) đưa ra cho một lòng dẫn rộng (b >> h).
Hướng x dương tính theo hướng dòng chảy. Hệ số Chezy giả thiết không đổi, nhưng
cũng có thể biểu thị bằng C = 18log(12h/ ks) với h =1/2(h0 + h1) = độ sâu nước trung bình
trên khoảng cách x1 – x0.
Hình 7.9. Tính toán đường co