Tổng quan nghiên cứu

Quá trình gia tăng dân số với tốc độ khoảng 1,1% mỗi năm cùng sự thu hẹp nhanh chóng của quỹ đất đô thị ven biển đang thúc đẩy các quốc gia phát triển hạ tầng vươn ra đại dương. Các công trình kết cấu nổi siêu lớn (Very Large Floating Structures - VLFS) với chiều dài từ 1000m đến 5000m và chiều rộng xấp xỉ 1000m đang trở thành giải pháp đột phá thay thế cho phương pháp lấn biển truyền thống. Nhờ ưu điểm không phá hoại hệ sinh thái đáy biển, không làm gián đoạn dòng hải lưu và tận dụng tối đa lực đẩy nổi của nước, VLFS dạng tấm phẳng (Pontoon hay Mega-Float) được ứng dụng rộng rãi làm sân bay nổi, bến cảng, kho chứa dầu khí và không gian đô thị nổi.

Tuy nhiên, với chiều dày kết cấu chỉ từ vài mét so với chiều dài hàng nghìn mét, VLFS có độ mềm uốn rất cao. Khi chịu tác động của tải trọng di động có vận tốc từ 5 m/s đến 70 m/s như máy bay Boeing 747 hạ cánh, hiện tượng tương tác thủy đàn hồi (hydroelasticity) diễn ra vô cùng phức tạp. Nghiên cứu này được thực hiện tại Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh trong giai đoạn 2018 - 2019 nhằm giải quyết bài toán phân tích động lực học kết cấu tấm nổi có bố trí hệ neo đàn hồi liên kết với đáy biển trong miền thời gian thực.

Phạm vi nghiên cứu tập trung khảo sát các dải thông số kỹ thuật đa dạng: vận tốc tải trọng từ 5 m/s đến 70 m/s, độ cứng hệ neo từ 10000 kN/m đến 60000 kN/m và độ sâu đáy biển từ 10m đến 58m. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi chứng minh hệ neo đàn hồi giúp kiểm soát và giảm độ võng dao động của tấm từ 18% đến hơn 28%, giảm lực cản sóng và đảm bảo an toàn tuyệt đối cho các công trình nổi siêu lớn.

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 tích hợp hai khung lý thuyết cơ bản trong cơ học công trình và cơ học thủy khí:

  • Lý thuyết tấm dày Reissner - Mindlin: Mô hình này xem xét đầy đủ biến dạng uốn và biến dạng trượt cắt trung bình với hệ số hiệu chỉnh cắt k' bằng 5/6, giúp mô tả chính xác ứng xử uốn cục bộ của kết cấu tấm nổi so với lý thuyết tấm mỏng cổ điển Kirchhoff.
  • Lý thuyết sóng tuyến tính và chất lỏng lý tưởng: Miền nước biển được giả thiết là chất lỏng không nén, không nhớt, chuyển động không xoáy và được biểu diễn thông qua hàm thế vận tốc phi thỏa mãn phương trình Laplace.

Mô hình nghiên cứu bao quát 4 khái niệm cốt lõi:

  1. Tương tác thủy đàn hồi (Hydroelastic interaction): Sự tương tác qua lại giữa biến dạng đàn hồi của kết cấu tấm và áp lực động lực học của chất lỏng bên dưới.
  2. Điều kiện biên động học mặt tiếp xúc (Kinematic boundary condition): Vận tốc theo phương đứng của đáy tấm và vận tốc bề mặt chất lỏng tại mọi điểm tiếp xúc luôn bằng nhau, triệt tiêu hoàn toàn khoảng hở giữa tấm và nước.
  3. Hệ neo đàn hồi (Elastic mooring system): Hệ thống dây neo hoặc giếng cọc được mô hình hóa thành các lò xo đàn hồi tuyến tính có độ cứng km, kết nối trực tiếp vào ma trận độ cứng tổng thể của tấm.
  4. Điều kiện biên vô cực (Far-field boundary): Vận tốc sóng triệt tiêu khi tiến ra xa kết cấu, đảm bảo năng lượng bức xạ không bị phản xạ nhân tạo vào miền tính toán.

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

Nguồn dữ liệu và cỡ mẫu mô phỏng số bao gồm 6 bài toán kịch bản đa chiều với hơn 200 phần tử hữu hạn bậc cao, bao phủ các trường hợp thay đổi vận tốc, vị trí đặt neo, số lượng liên kết neo, chiều dài tấm từ 320m đến 360m và độ sâu đáy biển.

Quy trình chọn mẫu và rời rạc hóa hình học được thực hiện như sau:

  • Phương pháp phân tử hữu hạn (FEM): Kết cấu tấm nổi được chia lưới bằng phần tử tứ giác 9 nút (Q9) với 3 bậc tự do tại mỗi nút (độ võng w, góc xoay quanh trục x và y). Thuật toán tích phân số sử dụng phép cầu phương Gauss 3 điểm (tổng cộng 9 điểm tích phân Gauss trên mỗi phần tử) nhằm khắc phục triệt để hiện tượng khóa trượt cắt (shear locking).
  • Phương pháp phần tử biên (BEM): Miền chất lỏng có kích thước mở rộng gấp 2 lần diện tích tấm nổi, được rời rạc hóa thành các phần tử biên panel 4 nút trên 4 vùng biên xác định gồm mặt đáy tấm (Hull Bed), mặt thoáng chất lỏng (Free Surface), mặt đáy biển (Sea Bed) và biên vô cực (Far Region).

Lý do lựa chọn phương pháp kết hợp FEM - BEM trong miền thời gian là vì FEM mô phỏng tối ưu độ cứng uốn phức tạp của tấm Mindlin, trong khi BEM chỉ cần chia lưới trên bề mặt biên mà không cần rời rạc hóa toàn bộ thể tích khối nước vô hạn, giúp tiết kiệm hơn 60% tài nguyên tính toán so với các phương pháp số khác. Toàn bộ thuật toán được lập trình tự động trên nền tảng phần mềm Matlab R2015b, triển khai nghiên cứu liên tục từ tháng 7/2018 đến tháng 12/2018 và kiểm chứng đối soát với phần mềm SAP2000 cùng công bố thực nghiệm của các học giả quốc tế.

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

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

Nghiên cứu số đã ghi nhận 4 phát hiện quan trọng về động lực học kết cấu tấm nổi:

  • Vận tốc tải trọng quyết định biên độ biến dạng: Khi vận tốc di động tăng từ 5 m/s lên vận tốc tới hạn 24 m/s, độ võng động tại điểm giữa tấm Z3 tăng vọt hơn 45%, hình thành vùng lõm sâu và tạo ra sóng cản uốn rõ rệt phía trước mũi tải trọng. Tuy nhiên, khi vận tốc vượt ngưỡng 30 m/s đến 70 m/s, chuyển vị cực đại có xu hướng dịch chuyển về phía đuôi tấm Z5 và biên độ uốn cục bộ giảm dần do thời gian tác dụng của tải trọng ngắn hơn chu kỳ dao động riêng của tấm.
  • Hệ neo đàn hồi triệt tiêu dao động vượt trội: Bố trí hệ neo đàn hồi tại các góc tấm với độ cứng 60000 kN/m giúp làm giảm chuyển vị đứng cực đại từ 18% đến 28,5% tại các điểm nút Z1 và Z3 so với kết cấu tấm nổi tự do không có neo.
  • Ảnh hưởng của độ cứng liên kết neo: Khi nâng độ cứng neo km từ 10000 kN/m lên 60000 kN/m, mô men uốn cực đại Mxx_max giảm khoảng 22%, đồng thời chuyển vị động tại thời điểm t = 4,5 giây giảm tuyến tính từ 12% đến 25%, giúp tấm ổn định nhanh chóng sau khi tải trọng đi qua.
  • Tác động của độ sâu đáy biển: Khi độ sâu đáy biển giảm từ 58m xuống 10m (vùng nước nông), hiệu ứng nén ép dòng chất lỏng đáy tấm làm gia tăng biên độ dao động chuyển vị thêm khoảng 35%, đòi hỏi hệ neo phải có khả năng chịu lực cắt và lực nhổ lớn hơn.

Thảo luận kết quả

Cơ chế vật lý của sự suy giảm dao động bắt nguồn từ việc hệ neo đóng vai trò như các gối tựa đàn hồi tiêu tán năng lượng, làm thay đổi dạng dao động riêng và ngăn chặn hiện tượng cộng hưởng thủy đàn hồi.

Khi so sánh với kết quả nghiên cứu của Ismail và mô hình kiểm chứng trên phần mềm SAP2000, sai số chuyển vị tĩnh và động lực học của thuật toán xây dựng chỉ dao động dưới 2,3%. Điều này khẳng định thuật toán BEM - FEM trong miền thời gian có độ chính xác và tin cậy rất cao.

Các kết quả phân tích số được trực quan hóa thông qua hệ thống biểu đồ lịch sử thời gian chuyển vị tại các nút Z1, Z3, Z5 trong khoảng thời gian từ 0 đến 10 giây và bảng ma trận so sánh nội lực uốn ứng với 5 cấp vận tốc (10 m/s, 16 m/s, 21 m/s, 24 m/s, 34 m/s). Biểu đồ phân bố lực cản theo phương X và Y cho thấy rõ ràng lực cản sóng tác dụng lên tấm giảm đi rõ rệt khi hệ neo được bố trí cân đối ở cả hai biên dọc và ngang của công trình.

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

Dựa trên các kết quả phân tích động lực học, nghiên cứu đưa ra 4 khuyến nghị kỹ thuật cụ thể:

  • Tối ưu hóa vị trí và mật độ bố trí neo: Các đơn vị tư vấn thiết kế công trình biển cần ưu tiên bố trí các cụm neo đàn hồi dày hơn tại vùng tiếp đất của đường băng nổi (nơi tải trọng máy bay tạo ra xung kích uốn lớn nhất), với khoảng cách phân bố dọc theo chiều dài tấm từ 300m đến 1000m nhằm mục tiêu giảm ít nhất 25% độ võng cục bộ, triển khai áp dụng ngay trong giai đoạn 2024 - 2026.
  • Lựa chọn vật liệu neo có độ cứng biến thiên: Các nhà thầu thi công nên áp dụng hệ neo cáp sợi tổng hợp kết hợp gối cao su giảm chấn có dải độ cứng đàn hồi từ 20000 kN/m đến 80000 kN/m để tăng khả năng tự điều chỉnh ứng suất, hướng tới mục tiêu nâng cao tuổi thọ công trình lên trên 50 năm dưới tác động xâm thực của môi trường biển.
  • Tích hợp hệ thống cảm biến quan trắc thời gian thực: Ban quản lý vận hành dự án cần lắp đặt hệ thống cảm biến gia tốc và đo biến dạng với tần số lấy mẫu trên 100 Hz tại các điểm xung yếu như mép biên và trung tâm tấm nổi nhằm giám sát liên tục phản ứng thủy đàn hồi trong suốt vòng đời khai thác 30 năm của công trình.
  • Xây dựng tiêu chuẩn thiết kế kết cấu nổi siêu lớn: Bộ Xây dựng và Bộ Giao thông Vận tải cần sớm chủ trì biên soạn quy chuẩn kỹ thuật quốc gia về tính toán thiết kế kết cấu VLFS có xét đến liên kết đàn hồi với đáy biển, hoàn thành lộ trình ban hành khung tiêu chuẩn trong vòng 24 tháng tới.

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

Luận văn là tài liệu tham khảo giá trị cao cho 4 nhóm đối tượng trọng tâm:

  • Kỹ sư kết cấu công trình thủy và công trình biển: Tiếp cận quy trình tính toán liên kết FEM - BEM trong miền thời gian, ứng dụng trực tiếp để thiết kế hệ neo và kiểm tra ổn định động cho các dự án đường băng sân bay nổi, bến cảng Mega-Float ngoài khơi.
  • Nhà nghiên cứu và học viên cao học chuyên ngành Xây dựng công trình: Khai thác mô hình toán học phần tử Q9 tấm Mindlin kết hợp thế vận tốc chất lỏng để mở rộng đề tài phân tích kết cấu nổi chịu tải trọng sóng thần, tải trọng gió bão hoặc va chạm tàu thuyền.
  • Chủ đầu tư và nhà quản lý dự án hạ tầng lấn biển: Nắm bắt các thông số kỹ thuật về độ cứng neo và chiều sâu nước biển từ 10m đến 60m để lập dự toán, tối ưu hóa chi phí đầu tư đê chắn sóng và giải pháp móng neo cho các tổ hợp đô thị nổi.
  • Cơ quan quản lý nhà nước và đơn vị thẩm tra thiết kế: Sử dụng các số liệu và biểu đồ lực cản trong luận văn làm cơ sở dữ liệu đối chuẩn phục vụ công tác thẩm định an toàn chịu lực và phê duyệt quy hoạch công trình nổi ngoài khơi.

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

Kết cấu nổi siêu lớn VLFS dạng Pontoon có gì khác biệt so với kết cấu phao nổi thông thường? VLFS dạng Pontoon sở hữu chiều dài vượt trội từ 1000m đến 5000m nhưng chiều dày chỉ từ 2m đến 5m, khiến độ mềm uốn đàn hồi vượt trội hoàn toàn so với chuyển động vật thể cứng. Trong thực tế, phản ứng của VLFS chủ yếu là biến dạng thủy đàn hồi (hydroelastic) phức tạp dưới tác động của tải trọng di động, thay vì chỉ dập dềnh đơn thuần như các phao nổi truyền thống.

Tại sao nghiên cứu lại kết hợp phương pháp FEM và BEM thay vì chỉ dùng duy nhất FEM? Sự kết hợp FEM cho kết cấu tấm và BEM cho miền chất lỏng mang lại hiệu quả tối ưu. Phương pháp BEM chỉ cần chia lưới trên các mặt biên của khối nước mà không cần rời rạc hóa toàn bộ thể tích vô hạn, giúp tiết kiệm hơn 60% dung lượng tính toán nhưng vẫn xác định chính xác trường thế vận tốc phi của chất lỏng lý tưởng.

Hệ neo đàn hồi có vai trò như thế nào trong việc giảm dao động của tấm nổi? Hệ neo đàn hồi hoạt động như các lò xo tiêu tán năng lượng, làm thay đổi tần số dao động riêng và ghìm giữ chuyển vị đứng của kết cấu. Kết quả mô phỏng cho thấy khi bố trí hệ neo với độ cứng 60000 kN/m, chuyển vị đứng cực đại tại trung tâm tấm giảm từ 18% đến 28,5% so với mô hình không neo.

Vận tốc của tải trọng di động ảnh hưởng như thế nào đến độ uốn của tấm? Khi vận tốc tải trọng tăng từ 10 m/s đến ngưỡng tới hạn khoảng 24 m/s, vùng lõm biến dạng mở rộng và chuyển vị động tăng hơn 40% do hiệu ứng sóng cản. Khi vận tốc vượt qua 30 m/s đến 70 m/s, độ võng tại vị trí tải trọng có xu hướng giảm do tải di chuyển nhanh hơn vận tốc truyền sóng uốn trong tấm.

Độ tin cậy của thuật toán tính toán trong luận văn được kiểm chứng ra sao? Chương trình tính toán lập trình bằng Matlab R2015b đã được kiểm chứng đối sánh độc lập với các công bố khoa học quốc tế của Ismail và mô hình phần tử hữu hạn chuẩn trên phần mềm SAP2000. Sai số chuyển vị tĩnh giữa các phương pháp chỉ dưới mức 2,3%, chứng minh độ chính xác và tính tin cậy tuyệt đối của mô hình.

Kết luận

  • Thiết lập thành công mô hình toán học tích hợp phương pháp phần tử hữu hạn FEM (phần tử Q9) và phương pháp phần tử biên BEM 3D trong miền thời gian để giải quyết bài toán tương tác thủy đàn hồi.
  • Chứng minh hệ neo đàn hồi giúp triệt tiêu dao động hiệu quả, làm giảm chuyển vị đứng cực đại của tấm từ 18% đến 28,5% và giảm mô men uốn hơn 20%.
  • Xác định chính xác vận tốc tới hạn của tải trọng di động ở mức xấp xỉ 24 m/s, tạo cơ sở dữ liệu quan trọng để thiết kế an toàn cho đường băng sân bay nổi.
  • Làm rõ quy luật gia tăng 35% biên độ chuyển vị của tấm khi độ sâu đáy biển giảm từ 58m xuống 10m ở vùng nước nông.
  • Cung cấp chương trình tính toán số tin cậy trên nền tảng Matlab R2015b với sai số kiểm chứng đối chuẩn dưới 2,3%.

Đóng góp chính của luận văn là đã hoàn thiện phương pháp luận và công cụ tính toán số chính xác cao để đánh giá ứng xử động lực học của kết cấu nổi siêu lớn có xét liên kết đàn hồi với đáy biển. Trong giai đoạn 2025 - 2027, hướng nghiên cứu tiếp theo sẽ tập trung mở rộng mô hình cho kết cấu tấm nổi đa module liên kết bán cứng và chịu tác động của sóng ngẫu nhiên phi tuyến. Các kỹ sư và nhà quản lý hãy áp dụng ngay phương pháp phân tích FEM - BEM kết hợp hệ neo đàn hồi để nâng cao tính an toàn và tối ưu hóa chi phí cho các công trình hạ tầng biển tương lai.