Tổng quan nghiên cứu

Không gian biển đang trở thành trọng tâm phát triển chiến lược toàn cầu khi đại dương chiếm hơn 70% diện tích bề mặt Trái Đất và gần 50% dân số tại các quốc gia công nghiệp hóa đang sinh sống trong phạm vi 1 km cách bờ biển. Nhu cầu mở rộng quỹ đất ven biển đối mặt với bài toán kinh tế và môi trường phức tạp. Các biện pháp quai đê lấn biển truyền thống chỉ phát huy hiệu quả kinh tế ở những vùng nước nông có độ sâu dưới 20 m. Khi độ sâu vượt quá giới hạn này hoặc địa chất đáy biển yếu, công trình nổi siêu lớn dạng phao (Very Large Floating Structures - VLFS Pontoon) dài hơn 60 m nổi lên như một giải pháp đột phá, tiêu biểu như mô hình thử nghiệm Mega-Float dài 1.000 m tại Vịnh Tokyo hay các căn cứ chứa dầu nổi Shirashima có dung tích 5,6 triệu kiloliters.

Vấn đề cốt lõi trong thiết kế kết cấu nổi siêu lớn là sự tương tác động lực học phức tạp giữa tấm kết cấu và môi trường nước xung quanh (hydroelasticity) dưới tác động của tải trọng di động như máy bay hạ cánh hoặc phương tiện vận tải. Do tỉ số giữa chiều dày vài mét và chiều dài hàng trăm mét rất nhỏ, công trình có ứng xử mềm uốn rõ rệt thay vì chuyển động vật rắn thuần túy. Luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng Công trình Dân dụng và Công nghiệp tại Đại học Bách Khoa - Đại học Quốc gia TP.HCM giai đoạn 2017 - 2018 tập trung xây dựng mô hình số kết hợp giữa phương pháp phần tử hữu hạn và phương pháp phần tử biên. Mục tiêu là phân tích chính xác ứng xử uốn động học của tấm nổi trực hướng kích thước 60 m x 60 m x 0,17 m, giúp giảm thiểu sai số thiết kế xuống dưới 1% và tối ưu hóa độ cứng kết cấu theo các phương chuyển động.

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 ứng dụng lý thuyết uốn tấm mỏng đàn hồi Kirchhoff mở rộng cho vật liệu trực hướng kết hợp cơ học chất lỏng lý tưởng. Khung lý thuyết bao gồm 3 trụ cột toán - cơ học chính:

  1. Lý thuyết tấm trực hướng đàn hồi tuyến tính: Quan hệ ứng suất - biến dạng được biểu diễn qua 4 hằng số đàn hồi độc lập gồm mô đun đàn hồi theo 2 phương chính $E_x, E_y$, hệ số Poisson $\nu_{xy}, \nu_{yx}$ và mô đun trượt $G_{xy}$. Phương trình vi phân cân bằng uốn bậc 4 mô tả chuyển vị thẳng đứng $w(x,y)$ chịu ảnh hưởng của các độ cứng uốn $D_x, D_y$ và độ cứng xoắn hữu hiệu $2H$.
  2. Lý thuyết thế năng động học chất lỏng lý tưởng: Nước biển được giả định là chất lỏng không nhớt, không nén được và chuyển động không xoáy. Trường thế vận tốc $\phi(x,y,z,t)$ thỏa mãn phương trình đạo hàm riêng thuần nhất Laplace $\nabla^2 \phi = 0$. Áp suất thủy động tác dụng lên đáy tấm được tích phân trực tiếp qua công thức tích phân Cauchy - Lagrange tuyến tính hóa.
  3. Lý thuyết tương tác động lực học chất lỏng - kết cấu (Hydroelasticity): Mô tả sự truyền sóng tương hỗ tại bề mặt tiếp xúc, trong đó vận tốc pháp tuyến của tấm và chất lỏng triệt để đồng nhất tại mọi vị trí tiếp giáp $\partial\phi/\partial n = \partial w/\partial t$.

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

Nghiên cứu áp dụng phương pháp số tổ hợp ghép nối trực tiếp (hybrid numerical scheme) giữa phương pháp phần tử hữu hạn (FEM) và phương pháp phần tử biên (BEM):

  • Phương pháp phần tử hữu hạn cho kết cấu: Bản nổi được rời rạc hóa bằng phần tử tấm chữ nhật 4 nút với 12 bậc tự do (mỗi nút gồm chuyển vị thẳng đứng $w$ và 2 góc xoay $\theta_x, \theta_y$). Phương trình dao động tổng thể dạng ma trận: $[M]{\ddot{w}} + [C]{\dot{w}} + [K]{w} = {F(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 kết cấu; ${P(t)}$ là véc-tơ áp lực thủy động từ chất lỏng.
  • Phương pháp phần tử biên cho miền chất lỏng: Phương trình tích phân biên Green 3D được thiết lập với hàm nghiệm cơ bản hàm Green cho miền bán vô hạn. Phương pháp BEM giúp giảm số chiều của bài toán từ 3D xuống 2D trên bề mặt tiếp xúc, giảm thiểu khối lượng bộ nhớ tính toán và xử lý triệt để điều kiện bức xạ sóng ở vô cực.
  • Quy trình rời rạc và chọn mẫu kiểm chứng: Toàn bộ thuật toán được lập trình trên môi trường Matlab. Mô hình tấm 60 m x 60 m được khảo sát qua 4 cấp lưới chia: 10x10 (100 phần tử), 20x20 (400 phần tử), 30x30 (900 phần tử) và 40x40 (1.600 phần tử). Lưới 30x30 với 961 nút và 2.883 bậc tự do được lựa chọn tối ưu, đảm bảo độ hội tụ cao và thời gian giải trong miền thời gian thực nghiệm 2,0 giây đạt độ chính xác tối ưu.

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 và khảo sát số trên mô hình tấm nổi trực hướng 60 m x 60 m x 0,17 m chịu tải trọng tập trung di động $P = 250\text{ kg}$ với 5 trường hợp biến thiên mô đun đàn hồi $E_1/E_2$ từ 0,5 đến 2,0 đã ghi nhận 4 phát hiện quan trọng:

  1. Sự tồn tại của dải vận tốc tới hạn: Chuyển vị võng lớn nhất của tấm nổi tăng mạnh khi vận tốc di chuyển của xe tăng từ 2,8 m/s (10 km/h) và đạt đỉnh cực đại trong dải vận tốc tới hạn từ 6,9 m/s đến 8,3 m/s (khoảng 25 đến 30 km/h). Khi vượt qua ngưỡng này và tăng lên 22,2 m/s (80 km/h), chuyển vị giảm nhanh chóng và duy trì trạng thái ổn định.
  2. Hiệu ứng gia cường của tính trực hướng: Khi nâng tỉ số mô đun đàn hồi dọc theo phương chuyển động từ $E_1 = 0,5E_2$ lên $E_1 = 2,0E_2$, chuyển vị lớn nhất của tấm giảm mạnh từ 12,465 mm xuống 6,696 mm (giảm 1,86 lần, tương đương giảm 46,28%). Ở dải vận tốc thấp 2,8 m/s, tăng độ cứng $E_1$ lên gấp 2 lần giúp giảm 18,5% chuyển vị, và tăng lên gấp 4 lần giúp giảm đến 54,79% chuyển vị võng.
  3. Quy luật phân bố tần số dao động tự nhiên: Ở mode dao động uốn thứ nhất, tần số riêng của tấm tăng từ 0,0432 Hz ($E_1 = 0,5E_2$) lên 0,0630 Hz ($E_1 = 2,0E_2$). Đáng chú ý, ở mode dao động thứ 6, tần số riêng tăng vọt từ 0,24529 Hz lên 0,35779 Hz, tương ứng mức tăng 45,9%.
  4. Biến dạng hình học dạng sóng động: Tải trọng di động tạo ra các vùng lõm cục bộ và gợn sóng biến dạng lan truyền trên mặt tấm. Tăng độ cứng theo phương di chuyển làm số lượng đỉnh sóng dọc trục giảm rõ rệt, tập trung năng lượng đàn hồi và triệt tiêu dao động lan truyền sang hai biên tự do.

Thảo luận kết quả

Các kết quả tính toán số phản ánh đúng bản chất vật lý của hiện tượng cộng hưởng thủy đàn hồi (hydroelastic resonance). Hiện tượng chuyển vị đạt cực đại tại dải vận tốc 6,9 - 8,3 m/s xảy ra do vận tốc di chuyển của tải trọng trùng với vận tốc pha của sóng mặt nước cảm ứng bên dưới đáy kết cấu. Khi tải trọng vượt qua vận tốc tới hạn, kết cấu không kịp phản ứng với xung lực động, dẫn đến biên độ dao động suy giảm.

Dữ liệu mô phỏng được trực quan hóa qua biểu đồ chuyển vị không gian 3 chiều và các đường trắc dọc $w(x)$ tại thời điểm $t = 2\text{ s}$. Sai số tính toán giữa chương trình Matlab tự lập trình và phần mềm thương mại SAP2000 trong bài toán tĩnh và dao động riêng chỉ dao động từ 0,04% đến 0,16% trên lưới 30x30. Đồng thời, mô hình kiểm chứng với bài toán tấm nổi đẳng hướng kích thước 9,75 m x 1,95 m x 0,053 m cho kết quả trùng khớp hoàn toàn với số liệu công bố quốc tế của tác giả Ismail năm 2004, khẳng định độ tin cậy tuyệt đối của phương pháp BEM-FEM.

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

Dựa trên các phát hiện từ phân tích thủy đàn hồi, 4 nhóm giải pháp kỹ thuật cụ thể được khuyến nghị áp dụng cho các dự án kết cấu nổi siêu lớn trong giai đoạn 2026 - 2030:

  • Tối ưu hóa bố trí sườn gia cường trực hướng: Các đơn vị thiết kế kết cấu cần định hướng bố trí hệ sườn dọc và cáp dự ứng lực theo phương di chuyển chính của phương tiện để đạt tỉ số $E_1/E_2 \ge 1,5$, giúp giảm trên 30% độ võng động mà không cần tăng chiều dày toàn khối của bản sàn bê tông.
  • Kiểm soát vận tốc khai thác giao thông: Ban quản lý vận hành công trình nổi như đường băng sân bay nổi hoặc cầu nổi vượt biển cần đặt biển giới hạn tốc độ, nghiêm cấm phương tiện di chuyển ổn định trong dải vận tốc nguy hiểm từ 25 km/h đến 30 km/h (6,9 - 8,3 m/s) nhằm triệt tiêu hoàn toàn rủi ro cộng hưởng thủy đàn hồi.
  • Tiêu chuẩn hóa công cụ tính toán ghép nối BEM-FEM: Các viện nghiên cứu chuyên ngành xây dựng công trình thủy cần tích hợp thuật toán BEM-FEM vào quy trình thẩm tra thiết kế nền tảng nổi, thay thế các mô hình dầm Winkler đơn giản vốn đánh giá sai lệch biên độ sóng uốn đến 25%.
  • Mở rộng mô hình hóa sóng biển ngẫu nhiên phi tuyến: Nhóm tác giả và các kỹ sư phát triển dự án cần tiếp tục nghiên cứu nâng cấp thuật toán sang bài toán tấm chịu tác động đồng thời của tải trọng di chuyển và phổ sóng biển đại dương ngẫu nhiên JONSWAP trong chu kỳ bão 50 năm.

Đố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à giá trị cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế kết cấu công trình biển và ven bờ: Nắm vững phương pháp phân bố độ cứng trực hướng tối ưu, áp dụng tính toán trực tiếp cho hệ mặt sàn các công trình bến nổi, cầu phao module và sân bay nổi.
  • Nhà nghiên cứu và chuyên gia cơ học tính toán: Khai thác thuật toán tổ hợp ma trận BEM-FEM trong miền thời gian, ứng dụng giải quyết các bài toán tương tác chất lỏng - kết cấu phức tạp mà các phần mềm thương mại khó xử lý.
  • Chủ đầu tư và nhà thầu phát triển hạ tầng nổi: Sử dụng dữ liệu khoa học để đánh giá tính khả thi kỹ thuật, tối ưu hóa khối lượng vật liệu xây dựng và lập phương án đảm bảo an toàn khai thác cho các dự án đô thị nổi ven biển.
  • Học viên cao học và sinh viên chuyên ngành Công trình Thủy - Xây dựng: Tiếp cận bài giảng mẫu mực về phương pháp phần tử biên, lý thuyết uốn tấm Kirchhoff và quy trình lập trình mô phỏng số trên Matlab.

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

Tại sao kết cấu nổi siêu lớn VLFS không thể coi là một vật thể rắn tuyệt đối khi tính toán? Do tỉ số giữa chiều dày khoảng 0,17 m đến vài mét so với chiều dài mặt bằng 60 m đến hàng nghìn mét là cực kỳ nhỏ, độ cứng uốn tổng thể của công trình tương đối mềm. Khi chịu tải trọng tập trung, phản ứng của tấm là biến dạng uốn đàn hồi cục bộ và tạo sóng thủy đàn hồi thay vì chuyển vị tịnh tiến cứng toàn khối.

Vì sao việc kết hợp FEM và BEM lại vượt trội hơn khi chỉ sử dụng một phương pháp FEM duy nhất? Nếu dùng FEM cho cả miền nước vô hạn, số lượng phần tử chất lỏng 3D sẽ lên đến hàng triệu bậc tự do, gây nghẽn bộ nhớ. Phương pháp BEM chỉ cần rời rạc hóa trên bề mặt tiếp xúc 2D, giải quyết hoàn hảo điều kiện bức xạ sóng ở vô cực và giảm trên 70% thời gian tính toán của hệ thống.

Độ cứng uốn trực hướng ảnh hưởng như thế nào đến độ võng của tấm nổi? Gia tăng độ cứng theo phương chuyển động $E_1$ giúp phân tán năng lượng uốn dọc trục hiệu quả. Khi tỉ số $E_1/E_2$ tăng từ 0,5 lên 2,0, độ võng cực đại giảm từ 12,465 mm xuống 6,696 mm, tương ứng mức giảm 1,86 lần, giúp công trình duy trì trạng thái cân bằng động học ổn định.

Vận tốc tới hạn của tải trọng di động trên tấm nổi là gì và mức độ nguy hiểm ra sao? Vận tốc tới hạn là vận tốc di chuyển khiến tần số kích thích ngoại lực trùng với tần số truyền sóng trong môi trường nước - kết cấu. Trong nghiên cứu này, vận tốc tới hạn nằm trong khoảng 6,9 m/s đến 8,3 m/s, khiến chuyển vị võng đạt giá trị lớn nhất và dễ gây nứt cục bộ bề mặt kết cấu.

Mô hình tính toán của luận văn đã được kiểm chứng độ tin cậy bằng những phương pháp nào? Chương trình tính toán được kiểm chứng kép: so sánh chuyển vị với kết quả thực nghiệm số của tác giả Ismail năm 2004 trên tấm 9,75 m x 1,95 m và đối chiếu tần số dao động riêng trên lưới 30x30 với phần mềm SAP2000, ghi nhận sai số tuyệt đối cực nhỏ dưới 0,16%.

Kết luận

  • Đóng góp học thuật xuất sắc: Luận văn đã xây dựng thành công thuật toán ghép nối chính xác giữa phương pháp phần tử biên (BEM) và phần tử hữu hạn (FEM) để giải quyết trọn vẹn bài toán thủy đàn hồi trong miền thời gian.
  • Làm sáng tỏ quy luật động học: Xác định chính xác dải vận tốc nguy hiểm 6,9 - 8,3 m/s và chứng minh việc tăng độ cứng trực hướng $E_1/E_2$ lên 2,0 lần giúp triệt tiêu tới 46,28% độ võng công trình.
  • Tối ưu hóa tài nguyên tính toán: Thiết lập hệ lưới rời rạc 30x30 đạt độ hội tụ hoàn hảo với sai số kiểm chứng dưới 0,16% so với chuẩn phần mềm quốc tế SAP2000.
  • Lộ trình triển khai kế tiếp: Hoàn thiện gói phần mềm phân tích chuyên dụng trong giai đoạn 2026 - 2028 và tích hợp tác động phi tuyến của sóng biển nhiều hướng phục vụ các dự án đảo nổi nhân tạo.
  • Khuyến nghị hành động: Các kỹ sư công trình biển và nhà quy hoạch hạ tầng cần áp dụng ngay khung phân tích trực hướng BEM-FEM để nâng cao hệ số an toàn và tối ưu hóa chi phí đầu tư cho các kết cấu nổi siêu lớn hiện đại.