Tổng quan nghiên cứu

Trong bối cảnh đô thị hóa nhanh và biến đổi khí hậu toàn cầu, nhu cầu mở rộng không gian sống và hạ tầng kỹ thuật ra biển trở thành xu thế tất yếu. Các giải pháp truyền thống như quai đê lấn biển đòi hỏi nguồn kinh phí khổng lồ, đặc biệt tại các vùng nước sâu trên 50 mét hoặc nơi có nền đất yếu, đồng thời làm suy giảm nghiêm trọng hệ sinh thái ven bờ. Nhằm giải quyết triệt để bài toán này, kết cấu nổi siêu lớn (Very Large Floating Structures - VLFS) dạng phao nổi (Pontoon) với chiều dài từ 1000 đến 5000 mét, bề rộng hàng trăm mét và chiều dày chỉ từ vài mét đã được nghiên cứu và ứng dụng tại các quốc gia tiên tiến như Nhật Bản, Hoa Kỳ và Hàn Quốc.

Nghiên cứu tập trung giải quyết bài toán phân tích động lực học của kết cấu tấm composite nổi quy mô lớn (kích thước miền khảo sát 4000 m x 3000 m x 5.7 m) trên nền nước sâu 300 mét, chịu tác động đồng thời của tải trọng tập trung di động 3883.54 kN và tương tác thủy đàn hồi (Hydroelasticity). Mục tiêu trọng tâm là thiết lập mô hình tính toán số tiên tiến, làm sáng tỏ ảnh hưởng của vận tốc di chuyển (từ 20 m/s đến 50 m/s), hướng di chuyển của tải trọng (từ 0 đến 45 độ) và tính trực hướng của vật liệu lên phản ứng động lực học của tấm nổi. Kết quả nghiên cứu mang lại ý nghĩa kinh tế và kỹ thuật to lớn, giúp giảm thiểu hơn 35% chi phí đầu tư xây dựng công trình biển và tối ưu hóa độ an toàn cho các công trình hạ tầng nổi như sân bay, bến cảng và khu đô thị nổi ngoài khơi.

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 xây dựng trên nền tảng ba khung lý thuyết chính nhằm mô phỏng chính xác sự tương tác phức tạp giữa kết cấu mỏng nhiều lớp và môi trường chất lỏng:

  1. Lý thuyết tấm mỏng đàn hồi tuyến tính Kirchhoff-Love trực hướng: Mô hình giả thiết các đường vuông góc với mặt trung hòa vẫn thẳng và vuông góc sau biến dạng, cho phép biểu diễn trường chuyển vị qua ba thành phần chuyển vị trên mặt trung hòa. Quan hệ ứng suất - biến dạng được thiết lập qua bốn hằng số đàn hồi độc lập với mô đun đàn hồi theo phương dọc $E_x = 11.9\text{ GPa}$, mô đun đàn hồi theo phương ngang $E_y = 5.95\text{ GPa}$ và hệ số Poisson $\nu = 0.13$, phản ánh chân thực đặc tính vách ngăn gia cường của kết cấu dạng hộp module.
  2. Lý thuyết thế vận tốc chuyển động chất lỏng tuyến tính: Nước biển được xem là chất lỏng lý tưởng (không nén được, không nhớt, chuyển động không xoáy), thỏa mãn phương trình vi phân Laplace trong miền không gian ba chiều với điều kiện động học tại đáy biển sâu 300 mét và mặt tiếp xúc tấm - nước.
  3. Lý thuyết tương tác thủy đàn hồi (Hydroelasticity): Thiết lập phương trình cân bằng liên hợp giữa lực đàn hồi, lực quán tính kết cấu và áp lực thủy động lực học của dòng nước tác động lên mặt đáy tấm.

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

Nghiên cứu sử dụng phương pháp kết hợp phần tử biên và phần tử chuyển động (BEM-MEM) nhằm khắc phục hoàn toàn các nhược điểm của phương pháp phần tử hữu hạn (FEM) truyền thống trong bài toán tải trọng di động trên miền vô hạn:

  • Phương pháp phần tử chuyển động (MEM): Hệ tọa độ được gắn cố định vào tải trọng di động, giúp loại bỏ hoàn toàn việc phải cập nhật lại vectơ tải trọng và lưới phần tử sau mỗi bước thời gian. Mô hình tấm được rời rạc hóa thành lưới 200 x 74 phần tử tấm Hermite bậc cao, đảm bảo tính liên tục của độ võng và góc xoay.
  • Phương pháp phần tử biên (BEM): Giảm số chiều bài toán chất lỏng từ miền khối 3 chiều xuống mặt biên 2 chiều thông qua nghiệm cơ bản hàm Green, với tổng số 4.440.000 panel biên phân bố trên bề mặt tiếp xúc và đáy biển.
  • Thuật toán giải và tích phân số: Sử dụng phương pháp cầu phương Gauss 2x2 để tính toán các ma trận phần tử và phương pháp tích phân từng bước Newmark trong miền thời gian với bước thời gian siêu nhỏ $\Delta t = 0.05\text{ s}$ xuyên suốt thời gian khảo sát 50 giây. Một lớp cản nhân tạo hấp thụ năng lượng tiêu tán sóng phản xạ được bố trí tại biên ngoài nhằm đảm bảo sóng truyền ra vô cực mà không bị dội ngược.

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

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

Quá trình mô phỏng số trên phần mềm Matlab đã mang lại các phát hiện cốt lõi về ứng xử động lực học của tấm composite nổi:

  • Hiện tượng cộng hưởng động và vận tốc tới hạn: Chuyển vị lớn nhất ($w_{max}$) và hệ số ứng xử động thay đổi phi tuyến mạnh theo vận tốc tải trọng. Khi vận tốc tăng từ 20 m/s lên dải vận tốc tới hạn 34 m/s – 37 m/s, hệ số ứng xử động tăng vọt hơn 42%, đạt đỉnh trước khi giảm dần ở các dải vận tốc siêu tới hạn (45 m/s đến 50 m/s).
  • Ảnh hưởng của góc di chuyển tải trọng: Do tính trực hướng của vật liệu ($E_x/E_y = 2$), khi góc hợp bởi hướng tải trọng và trục sợi tăng từ 0 độ lên 45 độ, độ võng cực đại tại tâm vùng lõm tăng thêm xấp xỉ 28%. Hướng di chuyển song song với trục có độ cứng lớn nhất ($E_x$) cho khả năng chịu uốn tối ưu nhất.
  • Độ trễ chuyển vị và biến dạng vùng lõm: Vị trí có độ võng cực đại không trùng với điểm đặt tải trọng mà luôn lệch về phía sau một khoảng trễ. Khi vận tốc tăng từ 20 m/s lên 50 m/s, khoảng cách trễ chuyển vị tăng từ 12.5 mét lên 48.2 mét (tăng gần 285%).
  • Sự biến đổi bước sóng kết cấu: Bước sóng kết cấu lan truyền ở phía sau tải trọng giãn dài tỷ lệ thuận với bình phương vận tốc di chuyển, trong khi bước sóng phía trước bị nén ngắn lại, tạo ra trường sóng bất đối xứng rõ rệt.

Thảo luận kết quả

Kết quả phân tích số cho thấy sự tương thích xuất sắc khi đối chiếu với các nghiên cứu thực nghiệm và bán giải tích trên thế giới, với mức sai số chuyển vị nhỏ hơn 3.5%. Dữ liệu chuyển vị theo thời gian có thể được trực quan hóa qua biểu đồ đường đồng mức 3D mô tả hình dạng phễu lõm chuyển động và bảng so sánh đa thông số.

Nguyên nhân của hiện tượng tăng vọt chuyển vị tại dải vận tốc 34 m/s – 37 m/s là do vận tốc di chuyển của tải trọng tiến gần đến vận tốc truyền sóng pha nhỏ nhất của hệ tấm - chất lỏng, gây ra hiện tượng cộng hưởng cục bộ. Khoảng trễ chuyển vị tăng nhanh theo vận tốc xuất phát từ quán tính của khối lượng nước bị chiếm chỗ bên dưới đáy tấm; lực nâng thủy động lực học không thể phản ứng tức thời tại điểm đặt lực mà cần độ trễ thời gian để tái phân phối áp lực đáy. Điều này khẳng định rằng đối với kết cấu tấm nổi lớn, việc bỏ qua tính tương tác thủy đàn hồi hoặc giả định nền Winkler/Pasternak đơn giản sẽ dẫn đến sai số nghiêm trọng trong tính toán kết cấu thực tế.

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

Dựa trên các kết quả phân tích định lượng, nghiên cứu đưa ra 4 nhóm giải pháp kỹ thuật cụ thể nhằm ứng dụng an toàn và hiệu quả kết cấu nổi siêu lớn:

  1. Kiểm soát dải vận tốc khai thác an toàn: Các đơn vị vận hành đường băng sân bay nổi hoặc cầu phao cần thiết lập giới hạn tốc độ cất/hạ cánh và lưu thông phương tiện dưới 27 m/s (xấp xỉ 97 km/h) hoặc vượt nhanh qua vùng vận tốc nguy hiểm 34 m/s – 37 m/s, giúp giảm thiểu 35% biên độ dao động uốn cục bộ của sàn kết cấu.
  2. Tối ưu hóa hướng bố trí kết cấu module: Khi thiết kế kết cấu tấm composite hoặc module vỉ thép hộp, các kỹ sư cần định hướng trục có độ cứng uốn lớn nhất ($E_x = 11.9\text{ GPa}$) song song 100% với trục đường băng chính (góc lệch bằng 0 độ), giúp tăng khả năng kháng uốn lên 28% và tiết kiệm đáng kể vật liệu gia cường.
  3. Ứng dụng hệ thống giảm chấn thủy lực biên: Chủ đầu tư và nhà thầu xây dựng cần tích hợp các tấm chắn sóng và hệ giảm chấn tiêu tán năng lượng tại vùng mép kết cấu nổi, mục tiêu hấp thụ tối thiểu 20% năng lượng sóng biển lan truyền trong vòng 12 tháng đầu đưa vào vận hành.
  4. Xây dựng tiêu chuẩn thiết kế kháng tải di động: Các cơ quan quản lý nhà nước như Bộ Xây dựng và Bộ Giao thông Vận tải cần sớm ban hành bộ quy chuẩn kỹ thuật cho kết cấu VLFS, quy định rõ ràng tải trọng thử nghiệm tiêu chuẩn 3883.54 kN và phương pháp kiểm tra tương tác thủy đàn hồi bắt buộc trước năm 2026.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng chuyên ngành:

  • Kỹ sư thiết kế công trình biển và dầu khí: Ứng dụng thuật toán BEM-MEM để phân tích ứng xử động lực học cho các bãi đáp trực thăng nổi, giàn khoan Semi-submersible, bến cảng nổi tiếp nhận tàu trọng tải trên 100.000 DWT.
  • Các nhà nghiên cứu và nghiên cứu sinh cơ học công trình: Khai thác mã nguồn lập trình Matlab và mô hình phần tử chuyển động để mở rộng nghiên cứu bài toán tấm composite chịu tải trọng gia tốc, tải phân bố tuần hoàn hoặc sóng phi tuyến.
  • Cơ quan quy hoạch không gian biển và hạ tầng đô thị: Sử dụng các số liệu kỹ thuật làm cơ sở tiền khả thi cho các dự án xây dựng sân bay nổi Mega-Float, đảo nhân tạo sinh thái quy mô 4000 m x 3000 m nhằm giải quyết bài toán quỹ đất ven biển.
  • Giảng viên và sinh viên sau đại học ngành Kỹ thuật Xây dựng: Làm tài liệu tham khảo chuyên sâu về phương pháp tích phân biên, lý thuyết tấm trực hướng và bài toán dao động phi dừng trong môi trường chất lỏng.

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

Ưu điểm vượt trội của phương pháp BEM-MEM so với phương pháp FEM truyền thống là gì?
Phương pháp BEM-MEM gắn hệ tọa độ chuyển động trực tiếp theo tải trọng, giúp lưới phần tử không cần chia lại và vectơ tải không cần cập nhật qua từng bước thời gian. Điều này giúp giảm hơn 70% dung lượng bộ nhớ máy tính và tiết kiệm 60% thời gian chạy mô phỏng đối với bài toán miền vô hạn kích thước 4000 m x 3000 m.

Tại sao cần phải xét đến tính trực hướng của tấm nổi composite?
Trong thực tế, kết cấu nổi siêu lớn được cấu tạo từ các module hàn ghép với hệ vách ngăn dọc và ngang bất đối xứng, dẫn đến độ cứng uốn hai phương chênh lệch nhau gấp 2 lần (11.9 GPa so với 5.95 GPa). Việc mô hình hóa trực hướng giúp phản ánh chính xác độ võng tấm, vốn có thể lệch tới 28% khi hướng tải thay đổi.

Vận tốc tới hạn của tải trọng di động trên tấm nổi là bao nhiêu và nguy hiểm như thế nào?
Vận tốc tới hạn được xác định trong khoảng 34 m/s đến 37 m/s (khoảng 122 - 133 km/h). Tại dải tốc độ này, tải trọng di chuyển cùng tốc độ với sóng uốn kết cấu, gây ra hiện tượng cộng hưởng làm hệ số phản ứng động tăng vọt 42%, dễ dẫn đến nứt vỡ cục bộ hoặc mất ổn định mỏi.

Khoảng trễ chuyển vị cực đại có ảnh hưởng gì đến phương tiện di chuyển bên trên?
Độ trễ chuyển vị (từ 12.5 m đến 48.2 m) khiến vùng lõm lớn nhất luôn nằm phía sau bánh xe hoặc vệt đáp máy bay. Điều này tạo ra độ dốc uốn liên tục ngay trước mũi phương tiện, làm gia tăng đáng kể lực cản lăn và kéo dài thời gian chạy đà cất cánh của máy bay trên sân bay nổi.

Mô hình chất lỏng trong nghiên cứu đã xử lý điều kiện biên vô cực như thế nào?
Nghiên cứu áp dụng lớp cản nhân tạo hấp thụ năng lượng tiêu tán dựa trên hàm bước Heaviside tại các biên $r = L_0$ và $s = L_0$. Giải pháp này triệt tiêu hoàn toàn 100% sóng phản xạ giả tạo dội ngược từ biên vào miền tính toán, đảm bảo tính chuẩn xác cho trường sóng tự do.

Kết luận

  • Nghiên cứu đã xây dựng thành công thuật toán và chương trình tính toán số BEM-MEM trên nền tảng Matlab, giải quyết trọn vẹn bài toán động lực học tương tác thủy đàn hồi của tấm composite nổi chịu tải trọng di động.
  • Xác định chính xác vùng vận tốc tới hạn nguy hiểm từ 34 m/s đến 37 m/s, nơi hệ số ứng xử động lực học đạt giá trị cực đại tăng hơn 42% so với tải tĩnh.
  • Chứng minh định lượng tính trực hướng của vật liệu ($E_x/E_y = 2$) và góc di chuyển tải trọng làm thay đổi tới 28% biên độ võng cục bộ của kết cấu nổi.
  • Làm sáng tỏ quy luật mở rộng bước sóng kết cấu phía sau tải trọng và sự gia tăng khoảng trễ chuyển vị cực đại lên tới 48.2 mét ở vận tốc 50 m/s.
  • Đóng góp cơ sở khoa học tin cậy cho việc tính toán, thiết kế và thẩm định an toàn các công trình kết cấu nổi siêu lớn (VLFS) tại Việt Nam trong giai đoạn 2025–2030; khuyến khích các đơn vị tư vấn và nghiên cứu áp dụng ngay giải pháp BEM-MEM để tối ưu hóa quy trình thiết kế công trình biển hiện đại.