Vấn đề: đầy một nửa không có nghĩa mọi đại lượng đều bằng một nửa
Trong một cống tròn chảy có mặt thoáng, hình dạng phần nước thay đổi theo độ sâu. Diện tích ướt, chu vi ướt và bề rộng mặt nước thay đổi theo những quy luật khác nhau. Bởi vậy, lấy lưu lượng đầy ống nhân với phần trăm chiều sâu chỉ đúng ở một vài trường hợp đặc biệt, không phải quy tắc chung.
Bài này xét tiết diện tròn lý tưởng, không có bùn lắng, với dòng ổn định đều và hệ số nhám không đổi theo độ sâu. Mục đích là xây một phép tính hình học–thủy lực để kiểm tra mô hình. Kết quả chưa bao gồm tổn thất cửa vào, cửa ra, hố ga hoặc nước dềnh từ hạ lưu.
Từ đường kính đến góc ở tâm
Ký hiệu D là đường kính trong cống, y là chiều sâu nước tính từ đáy trong, với 0<y≤D. HEC-RAS định nghĩa kích thước đứng lớn nhất của tiết diện tròn bằng đường kính, còn kích thước ngang tương ứng cũng bằng đường kính. [1] Đối với tính thủy lực, phải dùng kích thước lòng dẫn thực; đường kính ngoài không phải diện tích dành cho nước.
Gọi θ là góc ở tâm ôm phần chu vi ướt, tính bằng radian, nằm từ 0 đến 2π. Từ tam giác nối tâm đường tròn đến mép mặt nước, có:
cos(θ/2) = 1−2y/D
Do đó θ=2 arccos(1−2y/D). Nếu máy tính đặt chế độ độ thay vì radian, không được lấy trực tiếp con số góc ấy đưa vào công thức diện tích bên dưới. Hàm sin và bản thân θ phải dùng cùng quy ước radian để phép trừ θ−sinθ có ý nghĩa.
Phần diện tích nước bằng diện tích hình quạt trừ tam giác có dấu, cho:
A = D²(θ−sinθ)/8; P = Dθ/2
P chỉ là cung tròn tiếp xúc nước. Mặt thoáng không thuộc chu vi ướt. Bán kính thủy lực R=A/P, còn bề rộng mặt nước T=2√[y(D−y)]. T và P phục vụ các mục đích khác nhau, không được thay cho nhau trong công thức Manning.
Kiểm tra các trường hợp đặc biệt
Khi y=D/2, θ=π, A=πD²/8 và P=πD/2. Khi giới hạn y tiến tới D, θ tiến tới 2π, A tiến tới πD²/4 và P tiến tới πD. Trong cả hai trường hợp, R=D/4. Đây là một kiểm tra rất hữu ích: ở nửa đầy và đầy, diện tích khác nhau hai lần nhưng bán kính thủy lực bằng nhau.
Nếu bảng tính cho R ở nửa đầy bằng D/8, có thể đã cộng mặt thoáng vào P hoặc chia sai hệ số. Nếu A ở đầy không trở về diện tích hình tròn, cần kiểm tra đơn vị góc và hệ số D²/8 trước khi tính lưu lượng. Các kiểm tra hình học giúp tìm lỗi mà không cần dữ liệu công trường.
Liên hệ với lưu lượng
Trong hệ SI, với cùng hệ số nhám n và độ dốc năng lượng S:
Q(y) = A(y)·R(y)^(2/3)·√S/n
Q tính m³/s nếu A bằng m² và R bằng m. Tài liệu TxDOT phân biệt dòng có mặt thoáng với dòng có áp và sử dụng độ dốc năng lượng trong công thức sức cản. [2] Ở bài toán dòng đều đang xét, S bằng độ dốc đáy; ngoài điều kiện đó, cần xác định S từ mô hình thủy lực phù hợp.
Định nghĩa Qfull là giá trị Manning ở tiết diện đầy với cùng D, n và S. Khi y=D/2, R không đổi so với đầy và A bằng một nửa, nên Q/Qfull=0,5. Kết quả đặc biệt này không thể suy rộng thành Q/Qfull=y/D ở mọi độ sâu, vì R thay đổi phi tuyến ở các mức nước khác.
Ví dụ giả định: cống đường kính trong 1 m
Chọn D=1,00 m, n=0,013 và S=0,003. Các số liệu là giả định để kiểm tra công thức, không chỉ định loại cống hoặc tiêu chí thiết kế. Ở tiết diện đầy, A=0,785398 m², R=0,25 m và Qfull≈1,313209 m³/s.
Ở y=0,50 m, A=0,392699 m², R=0,25 m, cho Q≈0,656604 m³/s. Đúng như lập luận trên, lưu lượng bằng một nửa Qfull. Vận tốc trung bình ở hai trường hợp bằng nhau trong mô hình vì V=R^(2/3)√S/n và R bằng nhau.
Ở y=0,75 m, θ=4π/3≈4,188790 rad. Diện tích A≈0,631852 m², chu vi ướt P≈2,094395 m và R≈0,301687 m. Lưu lượng Q≈1,197486 m³/s, bằng khoảng 91,19% Qfull, chứ không phải 75%.
Nếu chỉ nhân 1,313209 với 0,75, kết quả khoảng 0,984907 m³/s. Sai khác gần 0,21258 m³/s xuất phát từ việc bỏ biến thiên R. Đây là sai số mô hình quy đổi, không phải sai số do làm tròn vài chữ số ở hàm lượng giác.
Vì sao gần đầy có thể cho Q lớn hơn giá trị đầy?
Ở y=0,90 m, tính được A≈0,744523 m², P≈2,498092 m, R≈0,298037 m. Q≈1,399615 m³/s, lớn hơn Qfull khoảng 6,58%. Ở y=0,95 m, Q≈1,411062 m³/s, khoảng 7,45% cao hơn Qfull.
Điều này không vi phạm bảo toàn khối lượng. Phần diện tích ở gần đỉnh được thêm vào khi đầy ống khá nhỏ, trong khi chu vi ướt tăng nhanh và R giảm về D/4. Trong mô hình giữ n và S không đổi, sự thay đổi sức cản tương đối có thể làm Q giảm khi mức nước tiến từ gần đầy tới đầy.
Không nên dùng kết quả này để mặc định cho cống vận hành sát đỉnh nhằm tối đa lưu lượng. Mô hình dòng đều có mặt thoáng không mô tả đầy đủ chuyển tiếp có áp, khí trong cống, nước dềnh và dao động mực nước. Giá trị Qfull ở một S cố định cũng không phải một giới hạn tuyệt đối cho mọi tình huống có áp; khi chênh năng lượng thay đổi, quan hệ lưu lượng thay đổi theo.
Độ nhạy theo đường kính và độ nhám
Nếu giữ tỷ lệ đầy y/D, A tỷ lệ D², P tỷ lệ D và R tỷ lệ D. Từ đó suy ra Q tỷ lệ D^(8/3) khi n và S cố định. Tăng D 10% làm Q nhân 1,1^(8/3), khoảng 1,2894, tức tăng gần 28,9%. Đây là so sánh ở cùng tỷ lệ đầy, không phải cùng độ sâu y tuyệt đối.
Nếu n tăng từ 0,013 lên 0,015 trong cùng hình học và S, Q giảm theo tỷ lệ 0,013/0,015≈0,86667, tức giảm 13,33%. Khi bùn lắng xuất hiện, chỉ tăng n có thể chưa đủ: phần tròn bị lấp làm thay đổi cả A và P. Không thể luôn thay cống thành một đường tròn có đường kính nhỏ hơn bằng D trừ chiều dày bùn, vì hình dạng còn lại không còn tròn đồng dạng.
Sai lệch đường kính thực, lớp lót hoặc vật cản tại mối nối có thể đáng kể khi cống đã gần mức nước khảo sát. Do đó, dữ liệu phục vụ đánh giá một tuyến đang sử dụng cần mô tả lòng dẫn thực, không chỉ đường kính danh nghĩa trong danh mục vật tư.
Từ một mặt cắt đến cả tuyến cống
USACE nêu khả năng thông qua của cống phụ thuộc cửa vào, thân cống và cửa ra. [3] Bài toán mặt cắt chỉ khảo sát phần thân trong điều kiện lý tưởng. Một cửa vào bị hạn chế có thể làm mực nước thượng lưu tăng dù thân cống còn diện tích; nước hạ lưu cao có thể làm tuyến chịu áp khác với tính dòng đều.
Đối với tuyến nhiều hố ga, cần xét cao độ đáy từng đoạn, tổn thất tại nối, dòng hợp lưu và điều kiện xả. Một độ sâu tính riêng từ Q tại mỗi đoạn không tự tạo thành đường mặt nước liên tục qua toàn mạng. Nếu mục tiêu là ngập sân theo thời gian, phải bổ sung thủy đồ vào và dung tích trữ tại các nút.
Hồ sơ có thể kiểm tra nên lưu D, y hoặc tỷ lệ đầy, θ, A, P, R, n và S. Kiểm tra nửa đầy–đầy giúp phát hiện lỗi hình học; đối chiếu cửa vào–cửa ra giúp phát hiện giới hạn của mô hình. Nhờ vậy, người đọc hiểu được kết quả nào thuộc mặt cắt và dữ liệu nào còn thiếu để đánh giá khả năng thoát nước của cả hệ thống.
Nguồn tham khảo
[1] HEC-RAS User’s Manual — Pipe Network Geometry — U.S. Army Corps of Engineers, Hydrologic Engineering Center. Truy cập ngày 09/11/2026.
[2] Hydraulic Design Manual — Section 2: Flow in Conduits — Texas Department of Transportation. Truy cập ngày 09/11/2026.
[3] HEC-RAS Hydraulic Reference Manual — Introduction to Culvert Terminology — U.S. Army Corps of Engineers, Hydrologic Engineering Center. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Trong một cống tròn chảy có mặt thoáng, hình dạng phần nước thay đổi theo độ sâu. Diện tích ướt, chu vi ướt và bề rộng mặt nước thay đổi theo những quy luật khác nhau. Bởi vậy, lấy lưu lượng đầy ống nhân với phần trăm chiều sâu chỉ đúng ở một vài trường hợp đặc biệt, không phải quy tắc chung.
Bài này xét tiết diện tròn lý tưởng, không có bùn lắng, với dòng ổn định đều và hệ số nhám không đổi theo độ sâu. Mục đích là xây một phép tính hình học–thủy lực để kiểm tra mô hình. Kết quả chưa bao gồm tổn thất cửa vào, cửa ra, hố ga hoặc nước dềnh từ hạ lưu.
Từ đường kính đến góc ở tâm
Ký hiệu D là đường kính trong cống, y là chiều sâu nước tính từ đáy trong, với 0<y≤D. HEC-RAS định nghĩa kích thước đứng lớn nhất của tiết diện tròn bằng đường kính, còn kích thước ngang tương ứng cũng bằng đường kính. [1] Đối với tính thủy lực, phải dùng kích thước lòng dẫn thực; đường kính ngoài không phải diện tích dành cho nước.
Gọi θ là góc ở tâm ôm phần chu vi ướt, tính bằng radian, nằm từ 0 đến 2π. Từ tam giác nối tâm đường tròn đến mép mặt nước, có:
cos(θ/2) = 1−2y/D
Do đó θ=2 arccos(1−2y/D). Nếu máy tính đặt chế độ độ thay vì radian, không được lấy trực tiếp con số góc ấy đưa vào công thức diện tích bên dưới. Hàm sin và bản thân θ phải dùng cùng quy ước radian để phép trừ θ−sinθ có ý nghĩa.
Phần diện tích nước bằng diện tích hình quạt trừ tam giác có dấu, cho:
A = D²(θ−sinθ)/8; P = Dθ/2
P chỉ là cung tròn tiếp xúc nước. Mặt thoáng không thuộc chu vi ướt. Bán kính thủy lực R=A/P, còn bề rộng mặt nước T=2√[y(D−y)]. T và P phục vụ các mục đích khác nhau, không được thay cho nhau trong công thức Manning.
Kiểm tra các trường hợp đặc biệt
Khi y=D/2, θ=π, A=πD²/8 và P=πD/2. Khi giới hạn y tiến tới D, θ tiến tới 2π, A tiến tới πD²/4 và P tiến tới πD. Trong cả hai trường hợp, R=D/4. Đây là một kiểm tra rất hữu ích: ở nửa đầy và đầy, diện tích khác nhau hai lần nhưng bán kính thủy lực bằng nhau.
Nếu bảng tính cho R ở nửa đầy bằng D/8, có thể đã cộng mặt thoáng vào P hoặc chia sai hệ số. Nếu A ở đầy không trở về diện tích hình tròn, cần kiểm tra đơn vị góc và hệ số D²/8 trước khi tính lưu lượng. Các kiểm tra hình học giúp tìm lỗi mà không cần dữ liệu công trường.
Liên hệ với lưu lượng
Trong hệ SI, với cùng hệ số nhám n và độ dốc năng lượng S:
Q(y) = A(y)·R(y)^(2/3)·√S/n
Q tính m³/s nếu A bằng m² và R bằng m. Tài liệu TxDOT phân biệt dòng có mặt thoáng với dòng có áp và sử dụng độ dốc năng lượng trong công thức sức cản. [2] Ở bài toán dòng đều đang xét, S bằng độ dốc đáy; ngoài điều kiện đó, cần xác định S từ mô hình thủy lực phù hợp.
Định nghĩa Qfull là giá trị Manning ở tiết diện đầy với cùng D, n và S. Khi y=D/2, R không đổi so với đầy và A bằng một nửa, nên Q/Qfull=0,5. Kết quả đặc biệt này không thể suy rộng thành Q/Qfull=y/D ở mọi độ sâu, vì R thay đổi phi tuyến ở các mức nước khác.
Ví dụ giả định: cống đường kính trong 1 m
Chọn D=1,00 m, n=0,013 và S=0,003. Các số liệu là giả định để kiểm tra công thức, không chỉ định loại cống hoặc tiêu chí thiết kế. Ở tiết diện đầy, A=0,785398 m², R=0,25 m và Qfull≈1,313209 m³/s.
Ở y=0,50 m, A=0,392699 m², R=0,25 m, cho Q≈0,656604 m³/s. Đúng như lập luận trên, lưu lượng bằng một nửa Qfull. Vận tốc trung bình ở hai trường hợp bằng nhau trong mô hình vì V=R^(2/3)√S/n và R bằng nhau.
Ở y=0,75 m, θ=4π/3≈4,188790 rad. Diện tích A≈0,631852 m², chu vi ướt P≈2,094395 m và R≈0,301687 m. Lưu lượng Q≈1,197486 m³/s, bằng khoảng 91,19% Qfull, chứ không phải 75%.
Nếu chỉ nhân 1,313209 với 0,75, kết quả khoảng 0,984907 m³/s. Sai khác gần 0,21258 m³/s xuất phát từ việc bỏ biến thiên R. Đây là sai số mô hình quy đổi, không phải sai số do làm tròn vài chữ số ở hàm lượng giác.
Vì sao gần đầy có thể cho Q lớn hơn giá trị đầy?
Ở y=0,90 m, tính được A≈0,744523 m², P≈2,498092 m, R≈0,298037 m. Q≈1,399615 m³/s, lớn hơn Qfull khoảng 6,58%. Ở y=0,95 m, Q≈1,411062 m³/s, khoảng 7,45% cao hơn Qfull.
Điều này không vi phạm bảo toàn khối lượng. Phần diện tích ở gần đỉnh được thêm vào khi đầy ống khá nhỏ, trong khi chu vi ướt tăng nhanh và R giảm về D/4. Trong mô hình giữ n và S không đổi, sự thay đổi sức cản tương đối có thể làm Q giảm khi mức nước tiến từ gần đầy tới đầy.
Không nên dùng kết quả này để mặc định cho cống vận hành sát đỉnh nhằm tối đa lưu lượng. Mô hình dòng đều có mặt thoáng không mô tả đầy đủ chuyển tiếp có áp, khí trong cống, nước dềnh và dao động mực nước. Giá trị Qfull ở một S cố định cũng không phải một giới hạn tuyệt đối cho mọi tình huống có áp; khi chênh năng lượng thay đổi, quan hệ lưu lượng thay đổi theo.
Độ nhạy theo đường kính và độ nhám
Nếu giữ tỷ lệ đầy y/D, A tỷ lệ D², P tỷ lệ D và R tỷ lệ D. Từ đó suy ra Q tỷ lệ D^(8/3) khi n và S cố định. Tăng D 10% làm Q nhân 1,1^(8/3), khoảng 1,2894, tức tăng gần 28,9%. Đây là so sánh ở cùng tỷ lệ đầy, không phải cùng độ sâu y tuyệt đối.
Nếu n tăng từ 0,013 lên 0,015 trong cùng hình học và S, Q giảm theo tỷ lệ 0,013/0,015≈0,86667, tức giảm 13,33%. Khi bùn lắng xuất hiện, chỉ tăng n có thể chưa đủ: phần tròn bị lấp làm thay đổi cả A và P. Không thể luôn thay cống thành một đường tròn có đường kính nhỏ hơn bằng D trừ chiều dày bùn, vì hình dạng còn lại không còn tròn đồng dạng.
Sai lệch đường kính thực, lớp lót hoặc vật cản tại mối nối có thể đáng kể khi cống đã gần mức nước khảo sát. Do đó, dữ liệu phục vụ đánh giá một tuyến đang sử dụng cần mô tả lòng dẫn thực, không chỉ đường kính danh nghĩa trong danh mục vật tư.
Từ một mặt cắt đến cả tuyến cống
USACE nêu khả năng thông qua của cống phụ thuộc cửa vào, thân cống và cửa ra. [3] Bài toán mặt cắt chỉ khảo sát phần thân trong điều kiện lý tưởng. Một cửa vào bị hạn chế có thể làm mực nước thượng lưu tăng dù thân cống còn diện tích; nước hạ lưu cao có thể làm tuyến chịu áp khác với tính dòng đều.
Đối với tuyến nhiều hố ga, cần xét cao độ đáy từng đoạn, tổn thất tại nối, dòng hợp lưu và điều kiện xả. Một độ sâu tính riêng từ Q tại mỗi đoạn không tự tạo thành đường mặt nước liên tục qua toàn mạng. Nếu mục tiêu là ngập sân theo thời gian, phải bổ sung thủy đồ vào và dung tích trữ tại các nút.
Hồ sơ có thể kiểm tra nên lưu D, y hoặc tỷ lệ đầy, θ, A, P, R, n và S. Kiểm tra nửa đầy–đầy giúp phát hiện lỗi hình học; đối chiếu cửa vào–cửa ra giúp phát hiện giới hạn của mô hình. Nhờ vậy, người đọc hiểu được kết quả nào thuộc mặt cắt và dữ liệu nào còn thiếu để đánh giá khả năng thoát nước của cả hệ thống.
Nguồn tham khảo
[1] HEC-RAS User’s Manual — Pipe Network Geometry — U.S. Army Corps of Engineers, Hydrologic Engineering Center. Truy cập ngày 09/11/2026.
[2] Hydraulic Design Manual — Section 2: Flow in Conduits — Texas Department of Transportation. Truy cập ngày 09/11/2026.
[3] HEC-RAS Hydraulic Reference Manual — Introduction to Culvert Terminology — U.S. Army Corps of Engineers, Hydrologic Engineering Center. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
