Các bài toán cơ bản kênh hở hình thang,tính toán theo phương pháp đối chiếu mặt cắt có lợi nhất về thủy lực và dòng chảy trong ống
CÁC BÀI TOÁN CƠ BẢN KÊNH HỞ HÌNH THANG.
Ta xét thấy: Q=f(n, i, b, h, m)
Tính kênh đã biết.
Bài toán 1: khi có n, i, b, h, m ta cần tìm Q
Ta tính những trị số A, C, R rồi thay vào (1-10) tìm được Q.
Bài toán 2: khi có n, Q, b, h ta cần tìm i.
Ta tính những trị số A, C, R rồi thay vào (1-9) tìm được theo công thức:
(1-29)
Bài toán 3: Khi có Q, i, b, h ta cần tìm n .
Thiết kế kênh mới.
Khi thiết kế kênh, cần tính chiều rộng và độ sâu mực nước kênh (b, h), cần thu thập các số liệu sau:
- Xác định độ dốc đáy kênh i, từ tuyến kênh theo bản đồ địa hình.
- Xác định hệ số nhám n và hệ số mái dốc m, căn cứ vào vật liệu lòng dẫn.
- Xác định lưu lượng Q, căn cứ vào nhu cầu sử dụng nước hay tiêu thoát nước được xác định ở các bài toán thủy nông, thủy văn công trình, cân bằng nước, v.v...
Sau khi xác định Q, m, n, i và chọn một trong các thông số, tùy từng trường hợp, thường gặp các bài toán có cách giải khác nhau như sau :
Bài toán 1 : Chọn β.
Từ công thức (1-10), tính theo Manning ta được:
, (m3/s) (1-30)
Kết hợp các công thức(1-15), (1-17) và (1-8) thay vào ta tính được:
, (m) (1-31)
b=βh, (m) (1-31a)
Bài toán 2 : Chọn R hay v.
Từ (1-14) và (1-16), ta có:
(a)
(b)
Để giải bài toán, tìm nghiệm b và h từ hệ phương trình trên, cần xác định A và P
+ Nếu biết R, từ (1-28) ta tính :
, (m2) (1-32)
, (m) (1-33)
+ Nếu biết v, từ (1-9) theo Manning ta có:
, (m/s) (1-34)
Nên: , (m) (1-35)
, (m2) (1-36)
Từ hệ phương trình, dùng phương pháp suy ra được như (1-26), sau đó khử h, ta được phương trình bậc hai:
m0h2 - Ph + A = 0 (1-37)
ở đó:
Giải phương (1-35) ta tìm được h.
(1-38)
Từ h1 và h2 thay vào (1-26), ta chọn nghiệm dương, chiều rộng b và độ sâu mực nước hợp lý làm nghiệm.
Chú ý : Bài toán có nghiệm khi :
- Điều kiện của (1-38) là P2 > 4m0A
- Ngoài ra ta biết rằng khi mặt cắt có lợi nhất về thủy lực, thì bán kính thủy lực và vận tốc là lớn nhất và diện tích mặt cắt là nhỏ nhất. Như vậy bài toán chỉ có lời giải khi R và v cho trước nhỏ hơn R và v lợi nhất về thủy lực.
Bài toán 3 : Chọn b (hay h). Tính h (hay b)
Từ (1-12), ta tính (1-39)
Từ (1-11) ta cũng có thể truy tìm nghiệm bằng cách lập bảng hoặc bằng đồ thị. Dùng cách lập trình trong Visual basic, Pascal hay dùng phần mềm Mathcad để tính.
TÍNH TOÁN THEO PHƯƠNG PHÁP ĐỐI CHIẾU MẶT CẮT CÓ LỢI NHẤT VỀ THỦY LỰC. (Phương pháp của AGƠRÔTSKIN)
Bài toán có b tìm h hay có h tìm b, thường phải giải đúng dần, cho nên việc tính toán dùng máy tính tay gặp khó khăn về thời gian và mức độ chính xác phụ thuộc người tính. Vì vậy trong phần này giới thiệu phương pháp tính của Agơrôtskin bằng cách lập bảng tra đối với mặt cắt hình thang.
Agơrôtskin đặt hệ sốđặc trưng mặt cắt hình thang, không thứ nguyên, biểu thị quan hệ giữa b, h, m, nghĩa là biểu thị hình dạng mặt cắt.
Từ đó xác định các yếu tố thuỷ lực theo đặc trưng mặt cắt, điều quan trọng mặt cắt hình thang lợi nhất về thuỷ lực, có giá trị đặc trưng mặt cắt lợi nhất bằng một.
Từ đó xác định được bán kính lợi nhất thuỷ lực, đặc biệt quan hệ mặt cắt lợi nhất về thuỷ lực và mặt cắt bất kỳ là hàm số phụ thuộc vào đặc trưng mặt cắt.
Quan hệ hình dạng mặt cắt.
Từ (1-14), đặt bề rộng trung bình hình thang:
(1-40)
nên: (1-41)
Từ (1-40) rút b thay vào (1-16) xắp xếp lại ta được :
(1-42)
ở đó đặt : (1-43)
Tính bán kính thuỷ lực theo(1-41) và (1-42), ta được
(1-44)
ở đó đặt: (1-45)
Từ các công thức trên, nếu ta biết hệ số đặc trưng mặt cắt, thì quan hệ giữa các yếu tố của mặt cắt sẽ được xác định như sau:
Từ (1-44) rút h ta được :
h=(1+ σ)R (1-46)
Từ (1-45) rút chiều rộng trung bình và thay (1-46) vào, ta được:
(1-47)
Từ (1-40) rút chiều rộng và thay (1-47) vào, ta được :
(1-48)
Từ (1-41) thay (1-46) và (1-47) tính lại diện tích theo công thức :
(1-49)
Suy ra (1-50)
Từ (1-46) và (1-48) ta tình được hệ số:
hay (1-51)
Đặc trưng của mặt cắt có lợi nhất về thủy lực.
Cũng như ở 1.3, xét mặt cắt lợi nhất, theo (1-50) ta biết rằng diện tích mặt cắt và mái dốc cho trước, nên mặt cắt lợi về thủy lực khi có R lớn nhất. Để R đạt gía trị lớn nhất ta xét đạo hàm sau :
Tính đạo hàm và giải ra ta được σ=1. Vậy điều kiện để có mặt cắt lợi nhất về thủy lực của hình thang là khi :
σLn=1 (1-52)
Từ (1-51) cho bằng 1, và chú ý công thức (1-43), ta sẽ tìm được công thức (1-27). Điều này cho thấy mặt cắt lợi nhất thuỷ lực hình thang có thể biểu thịquan hệ khác nhau nhưng bản chất là như nhau.
Quan hệ giữa mặt cắt có lợi nhất về thủy lực và mặt cắt bất kỳ.
Xét phương trình cơ bản, ta có:
Ta tính hệ số C theo công thức (1-5) của Pavơlôpski; còn A tính theo (1-49) thay vào công thức trên, chuý thay σLN=1 ứng với mặt cắt lợi nhất. Sau đó, tính tỉ số bán kính bất kỳ trên mặt cắt lợi nhất về thuỷ lực và rút gọn ta được:
(1-53)
Nếu xem y là hằng số, ứng với σ cho trước, ta tính được công thức (1-52). Nếu chia hai vế công thức (1-46) và (1-48) cho RLn ta được:
(1-54)
(1-55)
Theo Phoocơrâyme lấy y = 0.2, ta sẽ lập bảng các quan hệ giữa các đại lượng không thứ nguyên , , theo σ﹐ từ (1-53), (1-54), (1-55) ở (Phụ lục 1-2). Bảng này tự chúng ta cũng có thể lập bảng trên excel.
Từ phụ lục, nếu biết một trong các đại lượng, tra ra các đại lượng còn lại. Do đó, có thể tính các kích thước hình thang như b, h, R nếu biết bán kính lợi nhất vế thuỷ lực.
Xác định bán kính thủy lực.
Theo lưu lượng cho mặt cắt lợi nhất về thủy lực, ta có:
Agơrôtskin đã tính sẵn quan hệ:
(1-56)
Trong đó hệ số Chezy được tính theo công thức của tác giả và lập thành bảng (Phụ lục 1 -1)
Nếu tính C theo công thức của Maninh hay Phoocơrâyme, thì có thể tính rút trực tiếp ra RLn:
- Theo Maninh: (1-57)
- Theo Phoocơrâyme: (1-58)
Cách vận dụng cụ thể
Bài toán 1: Tìm h khi biết: Q, m, n, i và b.
+ Trước tiên xác định bán kính lợi nhất về thuỷ lực: RLn có thể dùng các công thức (1-57), (1-58) hoặc dùng phụ lục (1-1).
+ Lập tỉ: tra phụ lục (1-2) suy ra được:
+ Tính h theo công thức:
(1-59)
Bài toán 2: Tìm b khi biết: Q, m, n, i và h.
+ Trước tiên xác định RLn như trên
+ Lập tỉ: tra phụ lục (1-2) suy ra được:
+ Tính b theo công thức:
(1-60)
Bài toán 3: Tìm b và h, khi biết: Q, m, n, i và β
+ Xác định RLn như trên.
+ Tính đặc trưng mặt cắt hình thang theo công thức (1-51), tra phụ lục (1-2) suy ra được ,
+ Tính h và b theo công thức: (1-59) và (1-60)
Bài toán 4: Tìm b và h, khi biết: Q, m, n, i và R hoặc v.
+ Xác định RLn như trên.
+ Nếu có R thì lập tỉ số, tra phụ lục (8-3) suy ra được: ,
+ Tính h và b theo công thức: (1-59) và (1-60)
- Nếu biết v: Tính vận tốc theo Chezy, hệ số Chezy xác định theo Manning. Do đó tính bán kính thuỷ lực R theo công thức (1-35), tính ra b và h như trên.
DÒNG CHẢY TRONG ỐNG
Các yếu tố thuỷ lực
Công thức tính diện tích và chu vi mặt cắt hình tròn chảy lưng ống, tuy đơn giản nhưng ít được các tài liệu chứng minh.
Tính diện tích, xét 2 phần: diện tích cung tròn MHG và diện tích tam giác OMN, tức là:
trong đó:
d là đường kính mặt cắt hình tròn;
θ là góc được ghi chú trên hình 3. (rad)
Diện tích cung tròn MHG:
Diện tích phần tam giác OMG:
Vì xét tam giác vuông OMN, ta có:
ta lại có:
Do đó:
Hay:
cosθ =1-2a (1-61)
Đặt:
(1-61a)
Công thức (1-65) và (1-66), giúp chúng ta thiết lập mối quan hệ giữa độ sâu mực nước chảy lưng ống với đường kính ống tròn và góc θ đã đặt, để từ tính diện tích ướt và chu vi ướt.
Diện tích: (1-62)
Đặt: (1-62a)
Chu vi ướt (1-63)
Chiều rộng mặt thoáng B=dsinθ (1-64)
Bán kính thuỷ lực (1-65)
Công thức tính lưu lượng
Tính lưu lượng theo công thức Manning (1-30), thay (1-62) và (1-65), ta được:
(1-66)
(1-67)
Mặt cắt lợi nhất về thuỷ lực
Với i, n và d cho trước, ứng độ sâu mực nước trong ống là bao nhiêu để có lưu lượng lớn nhất khi:
Sau khi lấy đạo hàm hàm số trên, ta được phương trình:
2θ- 5θcos2θ sin2θ =0
Giải phương trình, ta được: θ=1510 hay a=0,94.
Tính vận tốc theo (1-34), thay bán kính thuỷ lực (1-64), ta được:
(1-68)
Với i, n và d cho trước, ứng độ sâu mực nước trong ống là bao nhiêu để có vận tốc lớn nhất khi:
Sau khi lấy đạo hàm hàm số trên, ta được phương trình:
- 2θcos2θ sin2θ =0
Giải phương trình, ta được: θ=1290 hay a=0,81
Các bài thường gặp
Bài toán 1: Bài toán thiết kế, có Q, n và i. Xác định đường kính ống.
Giải.
Từ công thức (1-66), cho thấy Q=f(n, i, d, a), vì vậy bài toán có 2 ẩn số là d và a, nhưng chỉ có một phương trình, nên tuỳ yêu cầu thực tế ta cần lưu lượng lớn thì lấy a=0,94, còn tính theo vân tốc lớn nhất lấy a=0,81.
Khi có a ta kính được θ và kA, tính theo công thức sau:
(1-69)
Bài toán 2: Bài toán kiểm tra, có Q, d, n và i. Xác định độ sau mực nước.
Giải.
Từ (1-67), ta tính được:
(1-70)
Có 2 cách để tìm nghiệm h:
- Cách 1: Phương pháp thử dần (mò nghiệm), tự chọn a tính θ và kA, ta vào biểu thức sau:
(1-71)
Tính đến khi nào h0(θ)≈ h(θ) thì gía trị a đó cần tìm.
- Cách 2: Tra bảng, từ công thức (1-61), (1-62) và (1-71) ta lập bảng tra
Từ công thức (1-70) tính được h0(θ) dựa vào bảng ta tra ra giá tri cần tím a, tính h theo công thức sau:
h=a.d (1-72)
Từ các công thức (1-61a), (1-61), (1-62a) và (1-71), tiến hành lập bảng bằng excel Phụ lục 1-3 để tra, thuân tiện trong việc tính toán bằng máy tính tay. Ta cũng thể dựa vào các công thức trên lập trình tính toán hay dùng phần mềm Mathcad.