Tổng quan nghiên cứu

Việt Nam sở hữu đường bờ biển dài hơn 3.260 km cùng tốc độ gia tăng nhu cầu không gian kinh tế biển đạt trên 15% mỗi năm. Các giải pháp san lấp lấn biển truyền thống thường bộc lộ nhiều hạn chế về chi phí khi thi công tại các vùng nước có độ sâu trên 15 mét hoặc khu vực đáy biển có tầng địa chất bùn sét yếu, đồng thời gây ảnh hưởng tiêu cực đến hệ sinh thái ngầm và đường bờ biển lân cận. Để khắc phục thách thức 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 quy mô mặt bằng dài hàng trăm mét đến hơn 1.000 mét nhưng bề dày chỉ từ 1 đến 5 mét đã trở thành giải pháp đột phá cho các công trình như sân bay nổi, bến cảng, trạm kho chứa dầu và bãi đáp trực thăng.

Tuy nhiên, kết cấu tấm nổi nhiều lớp làm việc đồng thời trong hai môi trường: bề mặt phía trên tiếp xúc trực tiếp với bức xạ mặt trời và không khí (nhiệt độ dao động từ 35°C đến trên 50°C), trong khi bề mặt đáy ngâm trong nước biển ổn định ở mức 20°C đến 25°C. Sự chênh lệch nhiệt độ này tạo ra gradient nhiệt lớn theo phương chiều dày, kết hợp cùng tác động của tải trọng di động có vận tốc thay đổi từ 2 m/s đến trên 20 m/s (như máy bay cất hạ cánh hoặc xe cơ giới), gây ra các ứng xử uốn và dao động thủy đàn hồi vô cùng phức tạp.

Luận văn thạc sĩ của tác giả Nguyễn Tấn Lộc, thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh năm 2022, đã tập trung giải quyết bài toán phân tích ảnh hưởng đồng thời của nhiệt độ và tải trọng di động đến ứng xử động lực học của tấm nổi hai lớp trên vùng nước nông. Bằng việc ứng dụng phương pháp phần tử chuyển động (Moving Element Method – MEM), nghiên cứu giúp giảm thiểu sai số dự báo chuyển vị động lực học từ 15% đến 25% so với phương pháp phân tích tĩnh thông thường, tạo tiền đề tối ưu hóa độ dày vật liệu và giảm thiểu từ 12% đến 18% chi phí xây dựng công trình nổ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 được xây dựng trên nền tảng tích hợp của ba hệ thống lý thuyết cơ học và thủy khí động lực học:

Thứ nhất, lý thuyết tấm mỏng Kirchhoff-Love được áp dụng để mô tả trạng thái biến dạng uốn phẳng của hai lớp tấm đẳng hướng song song. Vật liệu tấm tuân theo định luật Hooke với mô đun đàn hồi Young, hệ số Poisson và khối lượng riêng trên một đơn vị diện tích.

Thứ hai, mô hình nền đàn hồi loại Winkler được sử dụng để mô phỏng lớp lõi liên kết giữa tấm trên và tấm dưới với độ cứng liên kết và hệ số cản nhớt. Lớp lõi này truyền tương tác lực cắt và lực nén giữa hai lớp tấm khi kết cấu chịu uốn.

Thứ ba, lý thuyết sóng nước nông dạng tuyến tính (Linear Shallow Water Wave Theory) được ứng dụng để mô tả chuyển động của trường chất lỏng bên dưới tấm với độ sâu nước không đổi. Trường vận tốc của nước được biểu diễn thông qua hàm thế thỏa mãn phương trình vi phân liên tục dạng Laplace.

Các khái niệm then chốt trong nghiên cứu bao gồm: Gradient nhiệt độ phân bố tuyến tính theo chiều dày tấm, mô men uốn nhiệt tương đương, lực dọc trục sinh ra do giãn nở nhiệt trong mặt phẳng tấm, tương tác thủy đàn hồi (Hydroelasticity), và hệ tọa độ di động Euler-Lagrange đồng hành cùng tải trọng.

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

Nguồn dữ liệu và mô hình phân tích được phát triển hoàn toàn dựa trên phương pháp phần tử chuyển động (MEM) kết hợp phương pháp tích phân số Newmark trên nền tảng ngôn ngữ lập trình MATLAB.

Về mô hình rời rạc hóa, miền tính toán hữu hạn của tấm được chia thành các phần tử tấm chữ nhật Hermite 4 nút với tổng cộng 12 bậc tự do (mỗi nút gồm 3 bậc tự do: chuyển vị đứng và hai góc xoay). Miền chất lỏng được rời rạc bằng các phần tử 4 nút với 4 bậc tự do biểu diễn hàm thế vận tốc. Nhằm loại bỏ hoàn toàn hiện tượng sóng phản xạ từ biên nhân tạo của miền tính toán hữu hạn, tác giả đã tích hợp lớp cản ảo nhân tạo (Artificial Absorbing Boundary Layer) sử dụng hàm bước Heaviside bậc 4 theo hai phương toạ độ, giúp hấp thụ toàn bộ năng lượng sóng truyền ra xa.

Lý do lựa chọn phương pháp MEM thay vì phương pháp phần tử hữu hạn (FEM) truyền thống là vì MEM gắn hệ tọa độ trực tiếp lên tải trọng di động. Điều này giúp loại bỏ việc chia lại lưới phần tử theo từng bước thời gian, tiết kiệm hơn 40% chi phí tính toán máy tính. Thuật toán Newmark với tham số gia tốc trung bình đảm bảo tính ổn định số vô điều kiện cho hệ phương trình vi phân phi tuyến. Toàn bộ chương trình tính toán được kiểm chứng qua 5 bộ dữ liệu đối chuẩn độc lập với các công trình quốc tế đã công bố, cho thấy độ lệch chuẩn sai số luôn duy trì ở mức dưới 2.5%. Nghiên cứu được hoàn thành theo đúng tiến độ từ tháng 02/2022 đến tháng 06/2022.

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ố và phân tích các bài toán động lực học đã mang lại 4 phát hiện quan trọng:

Thứ nhất, ảnh hưởng của chênh lệch nhiệt độ: Khi độ chênh lệch nhiệt độ giữa hai bề mặt tấm tăng từ 0°C lên 40°C, chuyển vị đứng lớn nhất của tấm lớp trên tăng thêm từ 18.5% đến 28.3%. Nhiệt độ cao tại mặt trên làm sinh ra mô men uốn nhiệt dương, làm tăng độ võng tĩnh ban đầu và khuếch đại biên độ dao động khi có tải trọng di động đi qua.

Thứ hai, ảnh hưởng của tỷ lệ bề dày hai lớp tấm: Khi tăng chiều dày tấm trên từ giá trị ban đầu lên gấp 2 lần chiều dày tấm dưới, độ cứng chống uốn của lớp trên tăng gấp 8 lần. Kết quả là chuyển vị đứng cục bộ của lớp trên giảm từ 34.2% đến 46.5%, đồng thời giúp triệt tiêu hiện tượng tập trung ứng suất nhiệt tại vùng tiếp giáp nước.

Thứ ba, vai trò của độ cứng lớp lõi liên kết: Khi độ cứng lớp lõi tăng từ 0 lên mức 2 lần tích số trọng lượng riêng chất lỏng và gia tốc trọng trường (tương đương gần 20.000 N/m³), dao động của hai lớp tấm đạt trạng thái đồng pha hoàn toàn. Độ chênh lệch chuyển vị tương đối giữa lớp trên và lớp dưới giảm đi hơn 52.6%, giúp toàn bộ hệ kết cấu làm việc như một khối thống nhất.

Thứ tư, tác động của vận tốc di chuyển tải trọng: Khi vận tốc tải trọng tăng từ 2 m/s lên 4 m/s, đỉnh chuyển vị cực đại dịch chuyển lùi dần về phía sau điểm đặt lực khoảng 1.2 mét đến 2.5 mét do độ trễ truyền sóng trong chất lỏng, đồng thời biên độ chuyển vị đỉnh tăng thêm khoảng 12.4% đến 19.8% do hiệu ứng động học.

Thảo luận kết quả

Cơ chế vật lý của hiện tượng uốn nhiệt bắt nguồn từ sự giãn nở không đồng đều giữa thớ trên và thớ dưới của tấm nổi. Khi mặt trên bị đốt nóng bởi bức xạ mặt trời trong khi mặt dưới được làm mát liên tục bởi nước biển, phân bố nhiệt độ tuyến tính làm xuất hiện mô men uốn nhiệt tương đương $M_T$. Mô men này kết hợp với áp lực thủy tĩnh đẩy nổi và phản lực thủy động của lớp nước nông tạo nên trạng thái ứng suất uốn hai chiều rất phức tạp.

Dữ liệu mô phỏng trong luận văn được minh họa trực quan qua 45 biểu đồ biến thiên chuyển vị và 10 bảng tổng hợp số liệu chi tiết (từ Bảng 3-3 đến Bảng 3-10). Các đồ thị mặt cắt chuyển vị dọc theo trục di chuyển thể hiện rõ dạng sóng uốn lan truyền và sự suy giảm biên độ hoàn toàn tại vùng lớp cản ảo.

So với các nghiên cứu trước đây chỉ khảo sát tấm đơn lớp chịu tải di động mà bỏ qua yếu tố nhiệt độ, nghiên cứu này đã chứng minh rằng nếu không tính đến gradient nhiệt độ 30°C - 40°C, kết cấu nổi sẽ bị thiết kế thiếu an toàn từ 20% đến 30% về khả năng chịu võng. Ngoài ra, việc sử dụng phương pháp MEM trên nền nước nông tuyến tính giúp rút ngắn 65% thời gian xử lý ma trận tương tác rắn - lỏng so với phương pháp kết hợp phần tử hữu hạn và phần tử biên (FEM-BEM) truyền thống.

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

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

Thứ nhất, tối ưu hóa tỷ lệ chiều dày kết cấu: Các đơn vị thiết kế cần lựa chọn tỷ lệ chiều dày giữa lớp trên và lớp dưới ở mức từ 1.5 đến 2.0. Giải pháp này giúp giảm độ võng uốn cực đại tới 40% mà không làm tăng tổng khối lượng vật liệu quá 15%, cần áp dụng ngay trong giai đoạn thiết kế cơ sở (thời gian thực hiện từ 1 đến 3 tháng) do các kỹ sư kết cấu công trình biển chủ trì.

Thứ hai, lựa chọn vật liệu lõi đàn hồi giảm chấn: Khuyến nghị các nhà sản xuất vật liệu ứng dụng lớp đệm đàn hồi composite hoặc cao su kỹ thuật có độ cứng lớn hơn 20.000 N/m³ kết hợp hệ số cản nhớt cao. Mục tiêu nhằm giảm độ rung chấn truyền qua giữa hai lớp tấm xuống dưới 10%, triển khai chế tạo thử nghiệm trong vòng 6 tháng bởi các nhà thầu cung ứng vật liệu chuyên dụng.

Thứ ba, quy định giới hạn vận tốc di chuyển phương tiện: Cơ quan quản lý vận hành cảng và sân bay nổi cần ban hành quy chuẩn kiểm soát tốc độ phương tiện giao thông trên mặt sàn nổi ở mức dưới 15 m/s đối với xe tải nặng và kiểm soát gia tốc hạ cánh của máy bay nhằm duy trì hệ số xung kích động lực dưới 1.25, áp dụng liên tục trong suốt vòng đời khai thác công trình.

Thứ tư, thương mại hóa và mở rộng bộ công cụ tính toán MEM: Các viện nghiên cứu và trường đại học chuyên ngành xây dựng công trình thủy cần tiếp tục tài trợ kinh phí trong 12 tháng tới để phát triển thuật toán MATLAB hiện tại thành phần mềm thiết kế chuyên dụng có giao diện đồ họa, bổ sung các mô-đun phân tích sóng ngẫu nhiên phi tuyến và bức xạ nhiệt mặt trời thay đổi theo giờ.

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

Luận văn là tài liệu tham khảo chuyên sâu và hữu ích cho 4 nhóm đối tượng chính sau:

Thứ nhất, các kỹ sư thiết kế kết cấu công trình biển và ven bờ (Offshore & Coastal Engineers): Tài liệu cung cấp toàn bộ công thức giải tích ma trận và biểu đồ ứng xử chuyển vị để ứng dụng trực tiếp vào tính toán thiết kế sân bay nổi, cầu phao khẩu độ lớn, bến cảng nổi và bãi đáp trực thăng ngoài khơi.

Thứ hai, học viên cao học, nghiên cứu sinh và giảng viên ngành Kỹ thuật Xây dựng, Công trình Thủy: Luận văn là tài liệu học thuật mẫu mực về phương pháp phần tử chuyển động (MEM), phương pháp tích phân Newmark giải bài toán tương tác thủy đàn hồi (FSI), cùng các đoạn mã nguồn MATLAB hoàn chỉnh tại phần phụ lục.

Thứ ba, các chủ đầu tư, ban quản lý dự án phát triển kinh tế biển và lấn biển: Giúp các nhà quản lý có thêm cơ sở khoa học để đánh giá tính khả thi kinh tế - kỹ thuật của giải pháp kết cấu nổi Pontoon thay thế phương án san lấp lấn biển tại các vùng nước sâu từ 10 mét đến 30 mét có nền địa chất yếu.

Thứ tư, các cơ quan kiểm định chất lượng và tư vấn thẩm tra công trình giao thông thủy: Sử dụng các bảng tra số liệu và ngưỡng chuyển vị giới hạn trong luận văn làm tài liệu đối chiếu khi thẩm định mức độ an toàn chịu lực của các kết cấu nổi nhiều lớp dưới tác động đồng thời của khí hậu và tải trọng động.

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

Dưới đây là 5 câu hỏi thường gặp liên quan đến nội dung và ứng dụng thực tiễn của luận văn:

Tại sao cần phải xét đến ảnh hưởng của nhiệt độ khi thiết kế tấm nổi nhiều lớp? Bởi vì kết cấu nổi nằm giữa hai môi trường có nhiệt dung khác biệt lớn: mặt trên chịu nắng nóng không khí có thể lên tới 50°C, trong khi mặt dưới ngâm trong nước biển 20°C. Gradient nhiệt độ này sinh ra mô men uốn nhiệt tương đương $M_T$, làm gia tăng chuyển vị đỉnh của tấm thêm từ 18% đến 28%, dễ dẫn đến nứt gãy vật liệu nếu bị bỏ qua.

Phương pháp phần tử chuyển động (MEM) có ưu thế gì vượt trội so với FEM truyền thống? Phương pháp MEM sử dụng hệ tọa độ động di chuyển đồng hành cùng ngoại lực, biến bài toán di động trong không gian cố định thành bài toán tĩnh tương đối. Nhờ đó, MEM loại bỏ hoàn toàn nhu cầu tái chia lưới phần tử theo từng bước thời gian của FEM, giúp giảm hơn 40% thời gian tính toán mô phỏng mà vẫn đảm bảo sai số dưới 2.5%.

Lớp cản ảo nhân tạo trong mô hình số có chức năng gì? Trong mô phỏng số, kích thước lưới tấm phải bị giới hạn trong một miền hữu hạn. Lớp cản ảo nhân tạo sử dụng hàm suy giảm bậc 4 đặt tại các biên nhằm hấp thụ triệt để 100% năng lượng sóng kết cấu lan truyền ra ngoài, ngăn chặn hiện tượng sóng dội ngược làm méo mó kết quả chuyển vị tại tâm đặt tải trọng.

Tỷ lệ chiều dày giữa tấm trên và tấm dưới ảnh hưởng như thế nào đến độ bền công trình? Khi tăng bề dày lớp trên lên gấp 2 lần lớp dưới ($h_t = 2h_b$), độ cứng trụ chống uốn của lớp trên tăng gấp 8 lần, giúp giảm độ võng của tấm trên tới 46.5%. Tấm trên dày hơn đóng vai trò như một lớp áo giáp chịu lực chính và phân tán nhiệt, bảo vệ lớp đáy dưới nước luôn ổn định.

Kết quả nghiên cứu này có thể ứng dụng vào những dự án thực tế nào tại Việt Nam? Nghiên cứu có thể áp dụng ngay cho các dự án xây dựng bãi đáp trực thăng nổi tại quần đảo Trường Sa, cầu phao vượt biển tại vịnh Cam Ranh, bến đỗ du thuyền tại Nha Trang - Đà Nẵng, hoặc các trạm kho chứa xăng dầu nổi tại vùng ven biển Đồng bằng sông Cửu Long nơi có nền đất bùn rất yếu.

Kết luận

Năm đóng góp cốt lõi và kết luận nổi bật của luận văn bao gồm:

  • Xây dựng thành công mô hình phần tử chuyển động (MEM) giải quyết trọn vẹn bài toán tương tác thủy đàn hồi của tấm 2 lớp nổi trên vùng nước nông chịu tác động đồng thời của tải di động và gradient nhiệt độ.
  • Xác định định lượng sự gia tăng chuyển vị đỉnh từ 18.5% đến 28.3% khi chênh lệch nhiệt độ giữa hai bề mặt tấm tăng từ 0°C lên 40°C.
  • Chứng minh giải pháp tăng chiều dày lớp trên ($h_t = 2h_b$) và nâng cao độ cứng lớp lõi ($k_c \ge 20.000\text{ N/m}^3$) giúp giảm tới 46.5% độ võng động lực học của kết cấu.
  • Phát triển và kiểm chứng thành công bộ chương trình tính toán số trên MATLAB với độ tin cậy cao, độ lệch chuẩn so với dữ liệu quốc tế dưới 2.5%.
  • Cung cấp cơ sở khoa học và dữ liệu thiết kế quan trọng cho ngành xây dựng công trình biển Việt Nam trong việc phát triển các siêu kết cấu nổi VLFS thích ứng với biến đổi khí hậu.

Kế hoạch phát triển tiếp theo giai đoạn 2023 - 2025 là mở rộng mô hình tính toán cho tấm vật liệu FGM (Functionally Graded Materials) chịu tác động của phổ sóng biển ngẫu nhiên. Các kỹ sư và chuyên gia công trình biển hãy áp dụng ngay phương pháp tính toán MEM để tối ưu hóa độ an toàn và chi phí cho các dự án kết cấu nổi tương lai.