Tổng quan nghiên cứu

Khu vực bờ biển huyện Gò Công Đông thuộc tỉnh Tiền Giang sở hữu chiều dài bờ biển 32 km kéo dài từ cửa Soài Rạp đến cửa Tiểu, giữ vai trò chiến lược với 2.065 ha rừng phòng hộ ven biển bảo vệ sinh kế và hệ sinh thái thủy sản. Tuy nhiên, tình trạng xói lở và xâm thực bờ biển tại bãi biển Tân Thành diễn ra hết sức nghiêm trọng. Thống kê giai đoạn từ năm 1989 đến năm 2010 cho thấy chiều dày dải rừng ngập mặn suy giảm tới 250 m, tốc độ xói lở trung bình đạt 12 m/năm. Mặc dù địa phương đã nỗ lực trồng thêm 177 ha rừng phòng hộ nhưng thực tế vẫn có hơn 200 ha rừng bị sóng biển cuốn trôi hoàn toàn. Thực trạng này đặt ra yêu cầu cấp bách phải xây dựng hệ thống kè kiên cố bảo vệ mái đê biển trước tác động khắc nghiệt của sóng gió.

Nghiên cứu của tác giả Lê Hồng Hiếu thực hiện năm 2019 hướng tới giải quyết bài toán mô phỏng ứng xử cơ học của kè chắn mái nghiêng sử dụng cấu kiện bê tông tự chèn D27CM dưới tác động trực tiếp của sóng biển Gò Công. Mục tiêu trọng tâm của đề tài là xây dựng mô hình số 2 chiều tương đương, ứng dụng kỹ thuật phân tích tương tác rắn lỏng một chiều để xác định phân bố áp suất sóng thủy động theo thời gian, từ đó tính toán chuyển vị và kiểm tra độ bền ứng suất tương đương von-Mises trên nền đất đàn hồi.

Về mặt ý nghĩa thực tiễn, việc ứng dụng công nghệ mô phỏng số giúp thay thế các thử nghiệm mô hình vật lý thu nhỏ tốn kém, cắt giảm từ 50% đến 70% chi phí thử nghiệm và rút ngắn đáng kể thời gian thiết kế công trình. Kết quả nghiên cứu là luận cứ khoa học quan trọng giúp nâng cao hệ số an toàn đê biển lên trên 1,8, bảo vệ an toàn tuyến đê xung yếu trước mùa gió bão.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu kết hợp hai trụ cột lý thuyết cốt lõi trong cơ học hiện đại gồm Động lực học lưu chất tính toán và Cơ học kết cấu biến dạng thông qua bài toán phân tích tương tác rắn lỏng một chiều. Phương pháp thể tích hữu hạn được áp dụng để giải hệ phương trình bảo toàn khối lượng và phương trình động lượng Navier-Stokes mô tả chuyển động của sóng biển trong kênh hở. Trường dòng chảy được phân loại dựa trên chỉ số Reynolds, với điều kiện dòng chảy rối đạt Re từ 10.000 trở lên khi sóng va đập vào bờ kè.

Ở phần kết cấu công trình, phương pháp phần tử hữu hạn dựa trên nguyên lý công ảo được sử dụng để giải phương trình cân bằng động lực học. Tiêu chuẩn bền von-Mises được lựa chọn để đánh giá trạng thái ứng suất tổng hợp của vật liệu bê tông cốt thép mác M250 có mô đun đàn hồi E bằng 30 GPa, hệ số Poisson 0,2 và khối lượng riêng 2.300 kg/m3.

Khái niệm cấu kiện bê tông tự chèn D27CM giữ vai trò trung tâm trong thiết kế mái kè, với cơ chế ngàm khóa mềm 3 chiều giúp triệt tiêu hoàn toàn chuyển vị xoay và trượt ngang giữa các khối bê tông. Nền đất dưới mái kè được mô hình hóa theo lý thuyết nền đàn hồi với các hệ số nền đặc trưng, kết hợp lớp vải địa kỹ thuật TS 65 và lớp đệm đá dăm 1x2 dày 15 cm nhằm tiêu tán áp lực nước lỗ rỗng.

Phương pháp nghiên cứu

Dữ liệu hình học và kết cấu công trình được trích xuất từ hồ sơ thiết kế kỹ thuật của Sở Nông nghiệp và Phát triển Nông thôn tỉnh Tiền Giang. Dữ liệu khí tượng hải văn lấy theo các khảo sát thực tế tại Gò Công Đông, xác định điều kiện gió mùa Đông Bắc cấp 7 gây ra chiều cao sóng ven bờ trung bình từ 0,7 m đến 0,8 m, đỉnh sóng cực đại vượt mức 1,2 m với chu kỳ sóng dao động 4 đến 6 giây.

Mô hình mặt cắt đại diện 2 chiều được xây dựng với quy mô 44 cấu kiện tự chèn D27CM ghép nối liên tục trên mái dốc m = 4, trải dài từ cao trình +3,5 m xuống -0,8 m và liên kết với chân khay bệ chìm sâu đến cao trình -2,45 m. Quá trình chia lưới không cấu trúc được tinh chỉnh cục bộ tại bề mặt tiếp xúc rắn lỏng với kích thước phần tử 0,05 m nhằm đảm bảo độ hội tụ cao.

Lý do lựa chọn phương pháp tương tác rắn lỏng một chiều trên nền tảng phần mềm ANSYS Workbench là do độ cứng của kết cấu bê tông rất lớn, biến dạng đàn hồi nhỏ không làm thay đổi ngược lại trường dòng chảy của sóng. Quy trình thực hiện gồm hai giai đoạn liên hoàn: mô đun Fluent tính toán trường áp suất thủy động theo thời gian từ 0,05 giây đến 2,0 giây, sau đó dữ liệu áp lực bề mặt được chuyển trực tiếp sang mô đun Mechanical Transient Structural để phân tích động học kết cấu trong năm 2019.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phân tích số liệu mô phỏng thủy động lực học cho thấy áp suất sóng biển tác động lên mặt kè biến thiên mạnh mẽ theo các mốc thời gian t từ 0,05 giây đến 2,0 giây. Giá trị áp lực thủy động đạt đỉnh cực đại tại thời điểm t = 1,45 giây đến 1,50 giây khi mũi sóng đánh trực tiếp vào phần giữa mái dốc, tạo ra áp suất va đập lên tới 45 kPa. Tại các vị trí đỉnh tường kè ở cao trình +4,0 m, áp suất sóng giảm mạnh hơn 60% so với khu vực chân mái kè do năng lượng sóng đã tiêu hao đáng kể trong quá trình leo dốc.

Về phản ứng chuyển vị, toàn bộ 44 cấu kiện bê tông tự chèn D27CM làm việc đồng bộ và ổn định. Chuyển vị lớn nhất của kết cấu tập trung tại khu vực tiếp giáp giữa mái nghiêng và chân khay ống buy, với độ võng đàn hồi tối đa ghi nhận chỉ đạt 2,8 mm ở thời điểm 1,5 giây. Nhờ cấu tạo ngàm tự chèn 3 chiều, hiện tượng xô lệch hoặc bật trượt cục bộ giữa các viên bê tông đã giảm hơn 42% so với kết cấu tấm lát bê tông phẳng thông thường.

Trường ứng suất tương đương von-Mises phân bố tập trung chủ yếu tại các điểm tỳ ngàm cấu kiện và khối chân khay bệ chìm, đạt giá trị ứng suất cực đại 12,4 MPa. Giá trị này thấp hơn nhiều so với cường độ chịu nén tính toán 25 MPa của bê tông mác M250, mang lại hệ số an toàn kết cấu đạt 1,85. Hiệu suất tiêu giảm sóng phản xạ của mái nghiêng D27CM đạt trên 35%, giúp hạn chế tối đa nguy cơ xói chân công trình.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp kết cấu kè mái nghiêng Tân Thành duy trì độ ổn định cao dưới tác động của sóng cấp 7 nằm ở cơ chế truyền lực linh hoạt của cấu kiện D27CM. Khi nền đất xuất hiện hiện tượng lún không đều, ngàm khóa mềm cho phép các cấu kiện võng nhẹ theo mặt nền mà không phá vỡ liên kết mảng, triệt tiêu hoàn toàn khe hở tiếp giáp. Đồng thời, tầng đệm đá dăm 1x2 dày 15 cm cùng vải lọc TS 65 phát huy tối đa khả năng thoát nước ngược, ngăn chặn hiện tượng áp lực đẩy nổi phá hoại mảng kè.

So sánh với các nghiên cứu thực nghiệm mô hình vật lý của Nguyễn Viết Tiến và cộng sự năm 2013 hoặc công trình đê trụ rỗng Cà Mau của Viện Thủy công, phương pháp mô phỏng số FSI cho thấy mức độ tương đồng kết quả biến dạng và ứng suất với sai số dưới 5%. Khi trình bày trên đồ thị phân bố áp suất theo chuỗi thời gian 5 mốc quan sát và biểu đồ đường đẳng mức ứng suất von-Mises, kỹ sư có thể nhận diện trực quan vị trí xung yếu nhất tại cao trình -0,8 m để gia cường cục bộ, khắc phục nhược điểm tính toán tĩnh truyền thống vốn bỏ qua lực quán tính quá độ của sóng.

Đề xuất và khuyến nghị

Áp dụng phương pháp mô phỏng số tương tác rắn lỏng một chiều vào quy trình lập dự án và thẩm định thiết kế các tuyến kè bảo vệ đê biển xung yếu tại khu vực Đồng bằng sông Cửu Long. Mục tiêu cần đạt là rút ngắn 40% thời gian tính toán thủy lực và kiểm soát chính xác tải trọng va đập của sóng. Kế hoạch triển khai thực hiện trong giai đoạn từ năm 2020 đến năm 2022 do Ban Quản lý Dự án Đầu tư Xây dựng Công trình Nông nghiệp và Phát triển Nông thôn các tỉnh ven biển chủ trì.

Chuẩn hóa quy trình sản xuất và thi công lắp ghép khối bê tông tự chèn D27CM mác M250 kết hợp hệ thống vải địa kỹ thuật TS 65 và lớp đệm đá dăm dày 15 cm. Target metric đặt ra là đảm bảo độ sai lệch lắp ghép bề mặt dưới 3 mm, nâng cao độ bền công trình lên trên 30 năm trước môi trường xâm thực mặn. Đơn vị thi công và tư vấn giám sát cần hoàn thiện quy trình kiểm tra chất lượng định kỳ 6 tháng một lần bắt đầu từ năm 2021.

Xây dựng mạng lưới trạm quan trắc tự động theo dõi chuyển vị và áp lực sóng tại 5 phân đoạn kè trọng điểm dọc 32 km bờ biển Gò Công Đông. Mục tiêu là phát hiện sớm 100% hiện tượng sụt lún nền móng hoặc xói ngầm trước mùa gió Đông Bắc hàng năm. Sở Nông nghiệp và Phát triển Nông thôn tỉnh Tiền Giang cần phối hợp cùng các viện nghiên cứu chuyên ngành hoàn thành lắp đặt thiết bị trong giai đoạn 2022-2024.

Tối ưu hóa kết cấu chân kè bệ chìm sử dụng hàng ống buy bát giác đường kính 150 cm kết hợp lăng thể đá đổ ngoài mép nước nhằm nâng cao hệ số ổn định trượt tổng thể lên trên 2,0. Viện Khoa học Thủy lợi Miền Nam cùng các doanh nghiệp tư vấn xây dựng công trình thủy cần hoàn thiện bộ bản vẽ thiết kế điển hình vào quý 4 năm 2023.

Đối tượng nên tham khảo luận văn

Kỹ sư thiết kế công trình thủy lợi và công trình ven biển là nhóm thụ hưởng trực tiếp. Họ có thể khai thác toàn bộ quy trình thiết lập bài toán tương tác rắn lỏng FSI trên phần mềm ANSYS Workbench để tính toán tải trọng sóng động và ứng suất kết cấu cho các dự án kè biển tương tự.

Cán bộ quản lý đê điều và phòng chống thiên tai tại các địa phương ven biển, đặc biệt là khu vực Tiền Giang, Bến Tre và Cà Mau. Luận văn cung cấp bức tranh chi tiết về cơ chế tác động của gió mùa Đông Bắc cấp 7 và diễn biến xói lở, hỗ trợ lập phương án phân bổ ngân sách duy tu và gia cố đê điều hiệu quả.

Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật xây dựng công trình thủy, Cơ học kỹ thuật và Kỹ thuật tài nguyên nước. Tài liệu đóng vai trò là một case study mẫu mực về ứng dụng tích hợp phương pháp thể tích hữu hạn FVM và phương pháp phần tử hữu hạn FEM trong việc giải quyết các bài toán biên phức tạp.

Doanh nghiệp sản xuất cấu kiện bê tông đúc sẵn và nhà thầu xây lắp công trình thủy lợi. Các thông số kích thước hình học, độ bền uốn nén của cấu kiện D27CM và biện pháp thi công chân khay bệ chìm trong luận văn là cơ sở quan trọng để tối ưu hóa dây chuyền đúc khuôn và biện pháp tổ chức thi công ngoài hiện trường.

Câu hỏi thường gặp

Tại sao nghiên cứu lại lựa chọn mô hình tương tác rắn lỏng 1 chiều thay vì 2 chiều? Kết cấu kè bê tông tự chèn D27CM có độ cứng vật liệu rất cao với mô đun đàn hồi 30 GPa, chuyển vị tối đa đo được chỉ 2,8 mm nên không làm biến dạng đáng kể miền lưu chất của sóng biển. Lựa chọn mô hình 1 chiều giúp tiết kiệm hơn 60% thời gian tính toán trên máy tính mà vẫn duy trì độ chính xác đạt trên 95% so với mô hình thực tế.

Cấu kiện bê tông tự chèn D27CM có ưu điểm gì vượt trội so với tấm bê tông phẳng truyền thống? Khối bê tông D27CM dày 27 cm nặng 95,5 kg được thiết kế ngàm liên kết 3 chiều tạo thành mảng khóa mềm tự thích ứng với độ lún của nền đất. Khi chịu sóng va đập 45 kPa, cấu kiện không bị trượt hay xoay lật, triệt tiêu khe hở mặt tiếp giáp và giảm thiểu trên 80% nguy cơ xói ngầm lớp đất đắp bên dưới.

Đặc tính sóng biển tại Gò Công trong mùa gió Đông Bắc nguy hiểm như thế nào? Theo số liệu quan trắc khí hậu, gió mùa Đông Bắc cấp 7 tác động trực tiếp vào bờ biển hướng Bắc Nam của Gò Công Đông, tạo ra sóng ven bờ cao từ 0,7 m đến 1,2 m. Áp lực va đập của đợt sóng này lớn gấp 3 lần so với sóng mùa gió Tây Nam, là nguyên nhân chính gây mất 250 m rừng phòng hộ trong 21 năm.

Lớp đệm đá dăm 15 cm và vải địa kỹ thuật TS 65 đóng vai trò gì dưới mái kè? Lớp đá dăm 1x2 dày 15 cm kết hợp vải lọc TS 65 đóng vai trò tầng lọc ngược tiêu tán áp lực nước kẽ rỗng sau mỗi đợt sóng rút. Hệ thống này giữ lại toàn bộ hạt đất sét mịn của thân đê, ngăn ngừa hiện tượng rửa trôi cơ học và tăng độ ổn định động học cho toàn bộ mái dốc m = 4.

Làm thế nào để ứng dụng mô hình số này cho các vùng bờ biển có điều kiện địa chất khác? Các kỹ sư chỉ cần thay đổi thông số đầu vào trong mô đun ANSYS gồm vận tốc sóng, chiều cao sóng thiết kế từ 0,5 m đến 1,5 m và hệ số nền đất đặc trưng của từng địa phương như Bến Tre hay Trà Vinh. Quy trình chia lưới 0,05 m và thuật toán truyền dữ liệu FSI hoàn toàn có thể tái sử dụng nguyên vẹn.

Kết luận

  • Xây dựng thành công mô hình số 2 chiều mô phỏng chính xác tương tác rắn lỏng giữa sóng biển và kè bảo vệ mái đê biển Tân Thành, huyện Gò Công Đông.
  • Xác định trường áp suất sóng cực đại đạt 45 kPa tại thời điểm 1,45 giây dưới điều kiện gió mùa Đông Bắc cấp 7 với chiều cao sóng ven bờ từ 0,7 m đến 1,2 m.
  • Chứng minh tính ổn định vượt trội của 44 cấu kiện tự chèn D27CM với chuyển vị đàn hồi cực đại chỉ 2,8 mm và ứng suất von-Mises tối đa 12,4 MPa, đảm bảo hệ số an toàn 1,85.
  • Khẳng định giải pháp kết cấu chân khay bệ chìm đặt ở cao trình từ -0,8 m đến -2,45 m kết hợp hàng ống buy 150 cm tạo điểm tựa vững chắc cho toàn bộ mái dốc m = 4.
  • Đóng góp quy trình tính toán số hóa tiên tiến, giúp tiết kiệm trên 50% chi phí thử nghiệm thực nghiệm và mở ra hướng ứng dụng rộng rãi cho công trình đê kè biển Việt Nam.

Kế hoạch tiếp theo trong giai đoạn 2025-2030 là mở rộng bài toán mô phỏng 3 chiều toàn diện có xét đến kịch bản nước biển dâng do biến đổi khí hậu toàn cầu. Quý bạn đọc, kỹ sư và nhà quản lý hãy áp dụng ngay các giải pháp mô phỏng số chuyên sâu này để nâng cao chất lượng thiết kế và bảo vệ vững chắc hệ thống đê biển trước thiên tai!