Tổng quan luận án
Công trình kiểm soát nguồn nước tại các vùng đồng bằng ven biển và cửa sông giữ vai trò trọng yếu trong việc điều phối dòng chảy, ngăn mặn, giữ ngọt, tiêu thoát lũ và bảo vệ môi trường sinh thái phục vụ dân sinh và sản xuất nông - thủy sản. Tại vùng đồng bằng ven biển Việt Nam, đặc biệt là khu vực Đồng bằng sông Cửu Long, địa hình tương đối bằng phẳng, cao độ tự nhiên thấp, lòng sông rộng và chịu tác động trực tiếp của chế độ bán nhật triều hoặc nhật triều với biên độ lớn. Tình trạng xâm nhập mặn trong mùa kiệt ngày càng phức tạp do tác động kép của biến đổi khí hậu, nước biển dâng và sự vận hành của các bậc thang thủy điện ở thượng nguồn sông Mê Công (với tổng công suất khoảng 23.200 MW tại Trung Quốc và 2.275 MW tại Lào). Nhằm đáp ứng yêu cầu kiểm soát nước, hàng loạt hệ thống cống ngăn triều khẩu độ lớn đã và đang được triển khai như cụm công trình thủy lợi Cái Lớn – Cái Bé, hệ thống cống kiểm soát triều chống ngập Thành phố Hồ Chí Minh (Mương Chuối, Cây Khô, Phú Xuân, Tân Thuận, Phú Định), đập dâng sông Trà Khúc, cùng các dự án tại Bến Tre, Nam Sông Hậu và kênh Nguyễn Tấn Thành.
Các công trình kiểm soát nước vùng cửa sông ven biển đòi hỏi khẩu độ thông nước lớn, chế độ vận hành đóng mở không thường xuyên, chịu tác động của mực nước thay đổi hai chiều theo chu kỳ triều. Kết cấu cửa van phẳng kéo đứng nhịp lớn dạng dàn không gian ống thép (mặt cắt ngang hình thang hoặc vòm) học tập từ kinh nghiệm quốc tế (như công trình sông Ems tại Đức) đã được áp dụng tại Việt Nam. Tuy nhiên, thực tế thiết kế trong nước hiện nay chủ yếu dựa vào phương pháp "thử - đúng dần" mang tính kinh nghiệm, dẫn đến hệ số dự trữ độ bền giữa các bộ phận thiếu đồng đều, thanh bụng xiên có hệ số sử dụng vật liệu thấp hơn nhiều so với thanh đứng, làm tăng trọng lượng bản thân của cửa, kéo theo gia tăng chi phí chế tạo, vận chuyển, kích thước trụ pin và công suất thiết bị nâng hạ. Đồng thời, các tính toán thiết kế trước đây chưa đánh giá đầy đủ bài toán tương tác động lực học giữa khối nước thượng - hạ lưu với kết cấu van (FSI) khi chịu tải trọng động đất, cũng như ứng xử mỏi tại các mối nối chân dàn chịu chu kỳ áp lực nước lặp.
Mục tiêu nghiên cứu của luận án được xác định gồm 4 nội dung cụ thể:
- Nghiên cứu trạng thái ứng suất – biến dạng của kết cấu cửa van phẳng kéo đứng nhịp lớn trong điều kiện làm việc hai chiều.
- Nghiên cứu ảnh hưởng của thay đổi hình thức kết cấu thanh bụng xiên của dàn chính đến độ bền và độ cứng của cửa van phẳng kéo đứng nhịp lớn khi cửa van làm việc ở trạng thái đóng.
- Nghiên cứu ứng xử của kết cấu cửa van phẳng kéo đứng nhịp lớn khi chịu tác dụng của tải trọng động đất có xét đến tác dụng tương hỗ của khối nước trước và sau cửa van đang ở trạng thái đóng.
- Nghiên cứu tuổi thọ mỏi của kết cấu cửa van phẳng kéo đứng nhịp lớn khi chịu tác dụng của dao động mực nước trước và sau cửa van.
Đối tượng và phạm vi nghiên cứu:
- Đối tượng nghiên cứu: Kết cấu chịu lực dạng dàn vòm ống thép không gian của cửa van phẳng kéo đứng nhịp lớn trong công trình kiểm soát nước vùng triều tại Việt Nam.
- Phạm vi nghiên cứu: Hình dạng và kích thước hệ dàn ống thép của cửa van phẳng kéo đứng có chiều rộng thông nước trên 30 m ($B > 30\text{ m}$), biên độ chênh lệch cột nước thượng và hạ lưu nhỏ từ $1,0\text{ m} \div 2,5\text{ m}$, cửa van làm việc ở trạng thái đóng trên ngưỡng; tải trọng động khảo sát là gia tốc động đất theo thời gian; tuổi thọ mỏi được xác định thông qua trạng thái tập trung ứng suất và lan truyền vết nứt tại mối nối chân dàn chính liên kết với dầm biên.
Tổng quan tài liệu và vị trí của luận án
Luận án đã hệ thống hóa các dạng kết cấu cửa van nhịp lớn trên thế giới và trong nước, bao gồm 9 loại hình cơ bản: cửa van kéo đứng, cửa van nâng đứng, cửa van cung đơn trục ngang, cửa van viên phân xoay, cửa van bản sập, cửa van cung đôi trục đứng, cửa van cánh cửa, cửa van kéo ngang và cửa van ống mềm. Các công trình tiêu biểu sử dụng cửa van nhịp lớn được ghi nhận gồm có: Eastern Scheldt (Hà Lan, 62 khoang $B = 42\text{ m}$), St. Petersburg (Nga, cửa van phao $B = 200\text{ m}$, nâng đứng $B = 110\text{ m}$), Venice (Ý, cửa bản sập $3 \times 400\text{ m}$ và $1 \times 360\text{ m}$), sông Ems (Đức, cửa kéo đứng $1 \times 50\text{ m}$ và $4 \times 63,5\text{ m}$), Hartel (Hà Lan, $B = 98\text{ m}$), Nakdong (Hàn Quốc, $B = 95\text{ m}$), IJssel (Hà Lan, $B = 81,2\text{ m}$), Rhine (Hà Lan, $B = 81\text{ m}$). Tại Việt Nam, cửa van phẳng kéo đứng nhịp lớn chiếm ưu thế với các cống Cái Lớn (11 khoang $B = 40\text{ m}$), Mương Chuối (4 khoang $B = 40\text{ m}$), Cái Bé (2 khoang $B = 35\text{ m}$), Đập sông Trà Khúc ($14 \times 37,8\text{ m}$ và $5 \times 37,8\text{ m}$), Đập hạ lưu sông Dinh ($6 \times 37,8\text{ m}$), cống Rạch Móp, cống Nguyễn Tấn Thành ($B = 40\text{ m}$).
| Tên công trình – Quốc gia |
Năm xây dựng |
Nhịp cửa $B$ (m) |
Chiều cao $H$ (m) |
Cột nước $H_c$ (m) |
Trọng lượng (tấn) |
Hình thức kết cấu |
| Hartel – Hà Lan |
1996 |
98,0 |
9,30 |
7,90 |
640 |
Bản cong, dàn song song, vòm 2 chiều |
| Nakdong – Hàn Quốc |
2014 |
95,0 |
8,50 |
7,80 |
515 |
Bản cong, dàn nghiêng, vòm 2 chiều |
| IJssel – Hà Lan |
1976 |
81,2 |
11,50 |
4,50 |
590 |
Bản phẳng, dàn song song, vòm 1 chiều |
| Rhine – Hà Lan |
1981 |
81,0 |
11,55 |
4,55 |
450 |
Bản phẳng, dàn song song, vòm 1 chiều |
| Ems – Đức |
2002 |
64,2 |
$3 \times 8,5$; $1 \times 10,5$ |
6,60 |
$3 \times 400$; $1 \times 480$ |
Bản phẳng, dàn nghiêng, vòm 1 chiều |
Về tối ưu hóa kết cấu dàn thép: Lịch sử nghiên cứu ghi nhận công trình của Schmit (1960) kết hợp phần tử hữu hạn và quy hoạch phi tuyến; Venkayya (1970) với phương pháp tiêu chí tối ưu; Schmit và cộng sự (1974) phát triển kỹ thuật xấp xỉ; Goldberg và Samtani (1986) ứng dụng thành công thuật toán di truyền (GA) giải bài toán dàn 10 thanh. Tại Việt Nam và trên thế giới gần đây, Nguyễn Trọng Hà (2014) sử dụng phương pháp tìm kiếm ngẫu nhiên theo độ tin cậy; Phạm Hoàng Anh (2016) áp dụng thuật toán tiến hóa vi phân cải tiến (ANDE) tối ưu dàn phẳng với ràng buộc tần số; Phạm Văn Hưng và Đoàn Văn Duẩn (2017) áp dụng nguyên lý cực trị Gauss và nội lực hàm ẩn; Kim (2018) tối ưu phi tuyến dàn không gian; Lê Hoài Nam (2020) dùng thuật toán tối ưu bầy đàn (PSO) cho cầu dàn thép; Adil Baykasoğlu và Cengiz Baykasoğlu (2021) dùng thuật toán WSAR; MahdiAzizi và Uwe Aickelin (2022) nghiên cứu thuật toán CGO; Wang và Zhang (2002), Lluis Gil và Antoni Andreu (2001), Alemseged Gebrehiwot Weldeyesus (2020) nghiên cứu tối ưu hóa hình học, vị trí nút và cấu trúc liên kết dàn không gian dưới các điều kiện ràng buộc chuyển vị và ổn định tổng thể.
Về dao động và tương tác động lực học khối nước – kết cấu (FSI): Kolkman (1980) phân tích hiện tượng dao động cửa van Krimpen ($B = 80\text{ m}$) do dòng chảy đáy; Seong Haeng Lee (2014) thử nghiệm mô hình vật lý và mô hình số cửa van Nakdong ($B = 95\text{ m}$) dạng thấu kính; Xue Huifang và cộng sự (2012) dùng ANSYS chứng minh tần số dao động riêng của cửa van giảm tới 28,5% ở dạng dao động bậc 1 khi có nước; YAN Genhua và cộng sự (2020) chỉ ra tần số dao động riêng giảm từ 22% đến 40% do khối nước trước cửa; Khúc Hồng Vân (2018) dùng ADINA mô phỏng tương tác 3D sóng bề mặt thượng lưu tác dụng lên cửa van theo lý thuyết áp lực thủy động Westergaard.
Về phá hoại mỏi kết cấu cửa van: Shi Zhezhu (2008) chỉ ra nút liên kết chân dàn chính chịu tải trọng lặp có tuổi thọ mỏi 28 năm ($13,93 \times 10^6$ chu trình) khi đặt gân tăng cứng vuông góc trục thanh, và tăng lên 110 năm ($52,35 \times 10^6$ chu trình) khi bố trí gân tăng cứng song song trục thanh; Feng Hui Dong (2017) dùng ANSYS xác định hệ số tập trung ứng suất (SCF) của nút dàn ống thép tròn nối hàn với dầm hộp chữ nhật.
Khoảng trống nghiên cứu: Chưa có công trình nào nghiên cứu một cách hệ thống về việc cải tiến hình thức thanh bụng dàn chính làm việc hai chiều phù hợp với chế độ thủy triều Việt Nam; chưa có mô hình tích hợp thuật toán tối ưu đa mục tiêu (MOGA) cho dàn vòm ống thép không gian cửa van lớn; thiếu các đánh giá động lực học FSI đồng thời với hai khối nước (thượng lưu và hạ lưu) khi chịu gia tốc động đất; và chưa kiểm tra độ bền mỏi chi tiết chân dàn liên kết dầm biên theo mô hình con (submodeling).
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án sử dụng các hệ thống lý thuyết và công cụ phân tích số chuyên ngành công trình thủy và cơ học kết cấu:
-
Lý thuyết tối ưu hóa kết cấu:
Bài toán tối ưu hóa đa mục tiêu (Multi-Objective Genetic Algorithm - MOGA) được thiết lập nhằm tìm cực tiểu hàm mục tiêu trọng lượng kết cấu ($W_T$) hoặc thể tích vật liệu ($V$) dưới các ràng buộc trạng thái giới hạn sử dụng (chuyển vị cho phép $[d]$ hoặc độ võng $f \le [f]$) và trạng thái giới hạn cường độ (ứng suất tương đương $\sigma \le [\sigma]$). Các biến thiết kế (Design Variables - DVs) bao gồm đường kính ngoài và chiều dày thành ống của các nhóm thanh cánh, thanh bụng, thanh dàn ngang và thanh chịu trọng lượng.
-
Lý thuyết động lực học kết cấu và tương tác khối nước – kết cấu (FSI):
Phương trình vi phân chuyển động của hệ kết cấu nhiều bậc tự do có dạng:
$$[M]{\ddot{u}(t)} + [C]{\dot{u}(t)} + [K]{u(t)} = {P(t)}$$
Trong đó $[M]$, $[C]$, $[K]$ lần lượt là ma trận khối lượng, ma trận cản và ma trận độ cứng tổng thể; ${\ddot{u}(t)}$, ${\dot{u}(t)}$, ${u(t)}$ là véc tơ gia tốc, vận tốc và chuyển vị nút; ${P(t)}$ là véc tơ tải trọng động ngoài (gia tốc nền động đất). Phương trình được giải bằng phương pháp tích phân trực tiếp HHT (The Hilbert-Hughes-Taylor Method). Tương tác FSI giữa môi trường chất lỏng lý tưởng và kết cấu thép được mô hình hóa thông qua ma trận khối lượng quy đổi và điều kiện biên liên tục áp suất - chuyển vị tại mặt tiếp xúc cửa van và hai khối nước thượng - hạ lưu.
-
Lý thuyết phá hoại mỏi và tập trung ứng suất:
Quá trình phân tích mỏi dựa trên quy tắc cộng dồn tổn thất tuyến tính Miner, đường cong mỏi $S-N$ của vật liệu thép kết cấu $S355JR$, tính toán dải biên độ ứng suất ($\Delta \sigma$) và hệ số tập trung ứng suất (SCF). Mô hình nứt và phá hoại mỏi được chia thành giai đoạn hình thành vết nứt ban đầu ($N_{CI}$), giai đoạn phát triển nứt ổn định ($N_{SCG}$) và giai đoạn nứt không ổn định dẫn đến phá hủy ($N_{UCG}$).
-
Phương pháp thực tế và công cụ mô phỏng số:
- Công cụ tính toán: Phần mềm phần tử hữu hạn ANSYS sử dụng ngôn ngữ thiết kế tham số APDL (ANSYS Parametric Design Language). Tác giả đã lập trình hai mô đun tính toán chuyên dụng gồm
LIFTGATE (tính toán phân tích nội lực, chuyển vị, phổ ứng suất) và LGATE_OP (tối ưu hóa kích thước và tiết diện thanh dàn).
- Các loại phần tử PTHH sử dụng trong mô hình:
BEAM188: Phần tử dầm không gian 3 chiều (3D) có 2 nút, 6 hoặc 7 bậc tự do tại mỗi nút, có xét đến biến dạng trượt Timoshenko, dùng để mô phỏng toàn bộ hệ thanh dàn ống thép, dầm đứng và dầm biên.
SHELL63: Phần tử vỏ/mặt phẳng 3 chiều (3D) 4 nút, mỗi nút có 6 bậc tự do, dùng để mô phỏng bản mặt chắn nước bằng thép tấm.
FLUID30: Phần tử khối chất lỏng 3 chiều (3D) 8 nút, dùng để mô phỏng môi trường nước thượng lưu và hạ lưu tương tác với kết cấu van.
- Nguồn số liệu: Số liệu địa hình, thủy văn, địa chất, băng gia tốc động đất và các thông số thiết kế hình học lấy từ hồ sơ thiết kế cơ sở và thiết kế kỹ thuật của các dự án kiểm soát nước tại Đồng bằng sông Cửu Long (Cái Lớn – Cái Bé, Mương Chuối, Nguyễn Tấn Thành).
Nội dung chính theo từng chương
Chương 1: Tổng quan về cửa van nhịp lớn và những vấn đề cần nghiên cứu đặt ra với luận án
Chương 1 tổng hợp các khái niệm kỹ thuật và tiêu chí xác định cửa van nhịp lớn ($B/H \ge 4$, $B > 30\text{ m}$). Tác giả phân loại ưu nhược điểm của 9 hình thức cửa van lớn trên thế giới và phạm vi làm việc của từng loại. Phân tích kết cấu cửa van phẳng kéo đứng của công trình sông Ems (Đức) và cống Cái Lớn (Việt Nam), chỉ rõ đặc điểm: cửa van gồm 2 dàn chính dạng vòm đặt nghiêng góc $10^\circ \div 11^\circ$, bản mặt phẳng gắn hệ dầm phụ dọc hình thang và dầm đứng chữ T; dàn chịu trọng lượng bố trí ở phía hạ lưu.
Chương này chỉ ra hạn chế của sơ đồ dàn truyền thống: hệ thanh bụng xiên được bố trí để chịu lực cắt một chiều; khi cửa van làm việc trong môi trường chịu áp lực nước đảo chiều theo triều, thanh bụng xiên bị đổi dấu ứng suất (kéo thành nén hoặc ngược lại), dẫn đến hiệu suất sử dụng vật liệu của thanh bụng xiên rất thấp so với thanh đứng. Từ đó, xác định sự cần thiết phải nghiên cứu cải tiến hình thức thanh dàn và ứng dụng mô hình số tối ưu hóa.
Chương 2: Cơ sở lý thuyết tính toán kết cấu cửa van bằng mô phỏng số
Chương 2 trình bày toàn bộ cơ sở toán học cho 4 bài toán phân tích kết cấu:
- Bài toán tối ưu hóa hình dạng và kích thước dàn thép không gian dựa trên thuật toán di truyền đa mục tiêu (MOGA), thiết lập các phương trình quan hệ giữa hàm trọng lượng $W_T$ với biến thiết kế kích thước mặt cắt ngang $A_g$, mômen quán tính $I_g$, mômen chống uốn $W_g$.
- Bài toán động lực học kết cấu chịu kích động gia tốc theo thời gian, thiết lập hệ phương trình vi phân chuyển động tích phân từng bước bằng thuật toán HHT.
- Bài toán tương tác FSI giữa kết cấu van với 2 khối nước (thượng lưu và hạ lưu), sử dụng phương pháp ma trận khối lượng chất lỏng bổ sung và phần tử
FLUID30.
- Bài toán mỏi theo quy trình chu trình ứng suất kết hợp đường cong phá hoại $S-N$ của thép $S355JR$.
Chương cũng trình bày cấu trúc thuật toán và giao diện nhập liệu tham số của hai chương trình tự xây dựng trên nền tảng ANSYS/APDL là LIFTGATE và LGATE_OP.
Chương 3: Nghiên cứu thay đổi hình thức và tính toán tối ưu kết cấu dàn cửa van phẳng kéo đứng nhịp lớn
Chương 3 đề xuất cải tiến hình thức kết cấu dàn chính từ dạng truyền thống (gồm cả thanh đứng và hệ thanh bụng xiên ziczac) sang hình thức kết cấu mới: loại bỏ hệ thanh bụng xiên chịu lực một chiều, chỉ giữ lại hệ thanh đứng kết hợp thanh cánh thượng cong và thanh cánh hạ cong dạng vòm 2 chiều, hoặc bố trí lại thanh xiên đối xứng làm việc hai chiều.
Tác giả tiến hành khảo sát sự thay đổi nội lực (lực dọc $N$), ứng suất ($\sigma$) và chuyển vị ($f_z$ theo phương dòng chảy, $f_y$ theo phương thẳng đứng) khi thay đổi các thông số hình học:
- Chiều cao dàn chính tại giữa nhịp $B$: Khi tăng chiều cao $B$, chuyển vị ngang $f_z$ giảm nhanh, nhưng làm tăng chiều dài các thanh đứng và thanh dàn ngang; giá trị tỷ lệ $B/L \approx 0,12 \div 0,14$ cho hiệu quả chịu lực tối ưu.
- Khoảng cách giữa 2 thanh cánh hạ dàn chính $H$: Ảnh hưởng trực tiếp đến mômen chống uốn tổng thể và độ cứng chống vặn của cửa van.
- Độ cong và vị trí tương đối đỉnh thanh cánh thượng: Tạo dạng vòm ngược phân phối lại lực dọc trục, giúp giảm ứng suất kéo lớn nhất trong thanh cánh.
- Góc nghiêng của hai mặt phẳng dàn chính: Lựa chọn góc nghiêng hợp lý $10^\circ \div 11^\circ$ giúp cân bằng trọng tâm cửa van trùng với mặt phẳng tác dụng của xi lanh nâng hạ, hạn chế mômen lệch tâm gây kẹt bánh xe bên trong khe van.
Sử dụng chương trình LGATE_OP và thuật toán MOGA, tác giả thực hiện tối ưu hóa kích thước mặt cắt ngang của 6 nhóm thanh dàn ống thép (thanh cánh hạ, thanh cánh thượng, thanh đứng dàn chính, thanh bụng xiên, thanh dàn ngang, thanh dàn chịu trọng lượng). Kết quả so sánh giữa phương án truyền thống và phương án cải tiến đề xuất cho thấy: độ võng tương đối của kết cấu cửa van giảm từ $1/985$ xuống $1/1042$ (đều thỏa mãn điều kiện $[f] \le L/800$), trong khi tổng khối lượng kết cấu thép của hệ dàn giảm được từ 12% đến 18,5% so với thiết kế ban đầu.
Chương 4: Áp dụng kết quả nghiên cứu tính toán cho công trình thực tế
Chương 4 áp dụng toàn bộ quy trình tính toán và hình thức kết cấu đề xuất cho công trình Cống âu Nguyễn Tấn Thành (thuộc tỉnh Tiền Giang, kiểm soát nguồn nước kênh Nguyễn Tấn Thành và sông Tiền).
Thông số thiết kế của công trình:
- Bề rộng thông nước khoang cống: $B = 40,0\text{ m}$; chiều cao cửa van: $H_v = 8,0\text{ m}$.
- Vật liệu kết cấu: Thép hợp kim thấp cường độ cao $S355JR$ có giới hạn chảy $f_y = 355\text{ MPa}$, ứng suất cho phép $[\sigma] = 236,7\text{ MPa}$, mô đun đàn hồi $E = 2,1 \times 10^5\text{ MPa}$, hệ số Poisson $\nu = 0,3$.
- Các tổ hợp tải trọng tính toán:
- Tổ hợp 1 (mực nước thiết kế 1 chiều): MNTL $= +2,20\text{ m}$, MNHL $= -1,00\text{ m}$ (chênh lệch cột nước $\Delta H = 3,20\text{ m}$).
- Tổ hợp 2 (mực nước thiết kế đảo chiều): MNTL $= -0,50\text{ m}$, MNHL $= +2,00\text{ m}$ (chênh lệch cột nước $\Delta H = 2,50\text{ m}$).
- Tổ hợp đặc biệt: Tải trọng tĩnh kết hợp tải trọng động đất tác dụng theo phương ngang (băng gia tốc theo thời gian với đỉnh gia tốc cực đại $PGA = 0,05g$).
Mô hình PTHH gồm toàn bộ hệ thanh dàn (BEAM188), bản mặt chắn nước (SHELL63) kết hợp trường chất lỏng 3D (FLUID30) cho cả hai phía thượng và hạ lưu. Kết quả tính toán cho thấy:
- Chuyển vị lớn nhất theo phương dòng chảy $f_z = 24,6\text{ mm}$ (ứng với độ võng tương đối $f/L = 1/1626 < [1/800]$).
- Ứng suất chính lớn nhất trong hệ thanh dàn $\sigma_{\max} = 168,4\text{ MPa} < [\sigma] = 236,7\text{ MPa}$.
- Phân tích động lực học FSI: Khối nước ở cả 2 phía làm giảm tần số dao động riêng bậc 1 của kết cấu từ $2,84\text{ Hz}$ (trường hợp không nước) xuống $1,92\text{ Hz}$ (có xét FSI hai phía), đồng thời biên độ chuyển vị động lớn nhất trong thời gian động đất tăng khoảng 14,2% so với tính toán tĩnh quy đổi.
- Kiểm tra phá hoại mỏi tại vị trí nguy hiểm nhất (mối nối chân dàn chính liên kết với dầm biên, khảo sát chi tiết tại nút số 3339 bằng kỹ thuật submodeling): Ứng suất mỏi tương đương cực đại đạt $68,5\text{ MPa}$, tổn thất mỏi tích lũy sau 50 năm vận hành $D = 0,185 < 1,0$, khẳng định kết cấu nút đảm bảo an toàn chịu mỏi dài hạn.
Kết quả và những đóng góp mới
- Đóng góp về mặt lý luận và phương pháp:
- Làm sáng tỏ quy luật phân bố ứng suất và biến dạng của hệ dàn vòm ống thép không gian cửa van phẳng kéo đứng khi làm việc chịu áp lực nước hai chiều.
- Xây dựng thành công quy trình và mô hình tính toán số hoàn chỉnh tích hợp trong phần mềm ANSYS/APDL thông qua hai chương trình
LIFTGATE và LGATE_OP, cho phép tự động hóa quá trình phân tích tĩnh, tối ưu hóa kích thước mặt cắt theo thuật toán di truyền MOGA, phân tích động lực học tương tác chất lỏng – kết cấu (FSI 3D hai phía) chịu động đất và phân tích mỏi theo mô hình con (submodeling).
- Chứng minh được ảnh hưởng cản dịnh và khối lượng bổ sung của hai khối nước trước và sau cửa van đối với việc suy giảm tần số dao động riêng và gia tăng chuyển vị động của cửa van khi chịu kích động động đất.
- Đóng góp về mặt thực tiễn và giải pháp kỹ thuật:
- Đề xuất cải tiến hình thức hệ thanh bụng dàn chính cửa van kéo đứng nhịp lớn: loại bỏ hoặc sắp xếp lại hệ thanh bụng xiên đơn hướng truyền thống sang hệ thanh làm việc đối xứng hai chiều, giúp tăng hiệu quả sử dụng vật liệu của các thanh bụng và thanh đứng.
- Đưa ra bộ kích thước hình học hợp lý cho dàn chính: tỷ lệ chiều cao dàn trên nhịp $B/L \approx 0,12 \div 0,14$, góc nghiêng hai mặt phẳng dàn chính $10^\circ \div 11^\circ$, liên kết chân dàn sử dụng gân tăng cứng bố trí song song với trục thanh.
- Áp dụng trực tiếp vào phương án thiết kế cửa van cống âu Nguyễn Tấn Thành ($B = 40\text{ m}$), giúp giảm trọng lượng thép kết cấu hệ dàn trên 15% trong khi vẫn đảm bảo tuyệt đối các điều kiện về độ bền ứng suất, độ cứng chuyển vị và độ bền mỏi công trình.
Kiến nghị giải pháp: Cần áp dụng mô hình phân tích FSI 3D có xét cả 2 khối nước khi tính toán kháng chấn cho các cửa van cống vùng triều có nhịp $B \ge 30\text{ m}$; cấu tạo các mối nối hàn chân dàn chịu lực lớn bắt buộc phải sử dụng các bản gân gia cường định hướng dọc theo trục truyền lực của thanh ống thép.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của luận án:
- Nghiên cứu chỉ tập trung khảo sát trạng thái cửa van đóng hoàn toàn trên ngưỡng chịu chênh lệch mực nước tĩnh và tải trọng động đất; chưa xét đến giai đoạn mở cửa một phần chịu tác động thủy động lực học của dòng chảy xiết dưới đáy cửa hoặc hiện tượng dòng chảy tràn qua đỉnh cửa.
- Tải trọng sóng mới chỉ xét trong điều kiện sóng triều chu kỳ dài, chưa mô phỏng chi tiết phổ sóng ngẫu nhiên do gió bão tác động trực tiếp lên bản mặt cửa.
- Giả thiết khối chất lỏng trong phân tích FSI là môi trường đàn hồi tuyến tính không nhớt, bỏ qua hiện tượng nhớt của dòng bùn cát vùng cửa sông.
Hướng nghiên cứu tiếp theo:
- Nghiên cứu dao động cưỡng bức và hiện tượng cộng hưởng thủy lực (Vortex-Induced Vibration - VIV) của cửa van phẳng kéo đứng nhịp lớn ở các độ mở khác nhau dưới tác động của dòng chảy ngập và dòng chảy đáy.
- Nghiên cứu bài toán tối ưu hóa kết cấu có xét đến yếu tố ăn mòn kim loại theo thời gian trong môi trường nước lợ/mặn vùng cửa biển.
- Thực nghiệm đo đạc hiện trường các tham số dao động thực tế và ứng suất mỏi của cửa van tại các công trình đã đưa vào vận hành (như cống Cái Lớn, Nguyễn Tấn Thành) để kiểm chứng và hiệu chuẩn mô hình mô phỏng số.
Giá trị tham khảo
Luận án là tài liệu tham khảo chuyên sâu và hữu ích đối với các nhóm đối tượng:
- Các kỹ sư thiết kế công trình thủy lợi, thủy công và kết cấu công trình biển: Có thể sử dụng trực tiếp quy trình tính toán, cấu trúc tham số APDL của chương trình
LIFTGATE/LGATE_OP và các tỷ lệ hình học hợp lý ($B/L$, góc nghiêng dàn, cấu tạo mối nối chân dàn) để thiết kế các cửa van nhịp lớn từ $30\text{ m}$ đến $100\text{ m}$.
- Nghiên cứu sinh, học viên cao học và giảng viên ngành Kỹ thuật xây dựng công trình thủy, Cơ học kết cấu: Tham khảo phương pháp thiết lập bài toán tối ưu hóa kết cấu phi tuyến MOGA, kỹ thuật mô phỏng tương tác động lực học FSI 3D và phương pháp đánh giá tuổi thọ mỏi của kết cấu thép lớn.
- Các cơ quan quản lý đầu tư, ban quản lý dự án thủy lợi vùng ven biển: Có cơ sở khoa học để thẩm tra, đánh giá tính hợp lý về mặt kỹ thuật - kinh tế, khối lượng thép và tuổi thọ khai thác của các dự án kiểm soát mặn, ngăn triều.
Câu hỏi thường gặp
1. Thế nào là "cửa van nhịp lớn" theo định nghĩa và phạm vi nghiên cứu của luận án?
Trong luận án, cửa van nhịp lớn được định nghĩa là loại cửa van có kích thước bề rộng chắn nước $B$ lớn hơn nhiều lần chiều cao $H$ (thường $B/H \ge 4$) và có bề rộng nhịp thông nước $B > 30\text{ m}$. Ở quy mô này, nếu sử dụng kết cấu dầm đặc thông thường sẽ không đảm bảo tính kinh tế do thừa độ bền khi khống chế độ võng; vì vậy, kết cấu bắt buộc phải chuyển sang dạng dầm bụng rỗng hoặc hệ dàn vòm ống thép không gian.
2. Tại sao hình thức thanh bụng dàn chính truyền thống (kiểu sông Ems hay Cái Lớn) cần được cải tiến khi áp dụng cho công trình vùng triều tại Việt Nam?
Kết cấu kiểu sông Ems ban đầu được thiết kế cho công trình chịu áp lực nước chủ yếu một chiều. Khi áp dụng cho vùng triều Việt Nam có chế độ nước đảo chiều liên tục giữa thượng lưu và hạ lưu, hệ thanh bụng xiên ziczac bị đổi dấu ứng suất liên tục nhưng khả năng chịu lực cắt phân bố không đồng đều, dẫn đến hệ số sử dụng vật liệu của thanh bụng xiên rất thấp so với thanh đứng. Việc cải tiến bố trí lại thanh bụng (loại bỏ thanh xiên đơn hướng, chuyển sang làm việc hai chiều) giúp phân phối lại nội lực đồng đều và giảm từ 12% đến 18,5% khối lượng dàn thép.
3. Phương pháp phân tích tương tác khối nước – kết cấu (FSI) 3D hai phía của luận án có điểm gì khác biệt so với việc sử dụng công thức Westergaard truyền thống?
Công thức Westergaard truyền thống chỉ quy đổi áp lực nước tăng thêm thành khối lượng nước tương đương gắn vào một mặt bản van trong bài toán chuyển động một phía. Đối với cửa van cống vùng triều khi đóng, cửa van chìm đồng thời trong hai khối nước (thượng lưu và hạ lưu) có chiều sâu và tần số dao động riêng khác nhau. Mô hình FSI 3D sử dụng phần tử FLUID30 mô phỏng đầy đủ tác động chèn ép, khối lượng tham gia và lực cản của cả hai khối chất lỏng, phản ánh chính xác sự suy giảm tần số riêng của hệ (giảm từ 22% đến 40%) và sự gia tăng biên độ chuyển vị động khi có động đất.
4. Kết quả nghiên cứu mỏi đưa ra khuyến nghị gì đối với cấu tạo chi tiết mối nối chân dàn chính liên kết với dầm biên?
Nghiên cứu cho thấy vị trí tập trung ứng suất mỏi cao nhất nằm tại các đường hàn chuyển tiếp giữa chân ống thép thanh dàn chính với dầm biên. Luận án chỉ ra rằng các bản gân tăng cường tại nút liên kết bắt buộc phải được bố trí song song với trục của thanh ống thép chịu lực; giải pháp này giúp phân tán dòng truyền ứng suất dọc trục, giảm đột biến hệ số tập trung ứng suất SCF và nâng tuổi thọ mỏi của nút liên kết từ 28 năm lên trên 100 năm.
5. Công trình thực tế nào đã được lựa chọn để kiểm chứng kết quả nghiên cứu và kết quả tính toán cụ thể ra sao?
Công trình cống âu Nguyễn Tấn Thành ($B = 40\text{ m}$, $H_v = 8,0\text{ m}$, vật liệu thép $S355JR$) được chọn để áp dụng. Kết quả tính toán kiểm tra: chuyển vị ngang lớn nhất $f_z = 24,6\text{ mm}$ (độ võng $1/1626 \ll 1/800$), ứng suất lớn nhất $\sigma_{\max} = 168,4\text{ MPa} < [\sigma] = 236,7\text{ MPa}$, tần số dao động riêng FSI bậc 1 đạt $1,92\text{ Hz}$, và hệ số tổn thất mỏi sau 50 năm chỉ đạt $0,185$, đáp ứng đầy đủ các yêu cầu kỹ thuật và tiết kiệm đáng kể khối lượng kết cấu thép.
Kết luận
Luận án tiến sĩ của tác giả Trần Xuân Hải đã giải quyết trọn vẹn bài toán cơ sở khoa học và thực tiễn trong việc cải tiến hình thức kết cấu hệ dàn vòm ống thép cho cửa van phẳng kéo đứng nhịp lớn ($B > 30\text{ m}$) phục vụ các công trình kiểm soát nguồn nước vùng triều. Thông qua việc tích hợp thuật toán tối ưu hóa di truyền đa mục tiêu (MOGA), lý thuyết tương tác động lực học khối nước – kết cấu (FSI 3D hai phía) và mô hình hóa mỏi cục bộ trong môi trường ANSYS/APDL, công trình đã đề xuất giải pháp hình học hợp lý giúp giảm khối lượng thép từ 12% đến 18,5% mà vẫn đảm bảo an toàn tuyệt đối về độ bền, độ võng và tuổi thọ công trình. Kết quả áp dụng cho công trình cống âu Nguyễn Tấn Thành khẳng định tính đúng đắn, độ tin cậy và khả năng ứng dụng trực tiếp của phương pháp vào công tác tư vấn thiết kế các công trình thủy lợi, ngăn mặn quy mô lớn tại Việt Nam.