Tổng quan nghiên cứu

Trong bối cảnh biến đổi khí hậu toàn cầu và áp lực gia tăng dân số đô thị tại hơn 70% các quốc gia có đường bờ biển dài, nhu cầu mở rộng không gian sống ra biển đang trở thành xu thế tất yếu. Tuy nhiên, các giải pháp lấn biển truyền thống bằng san lấp đất thường tiêu tốn chi phí khổng lồ và tiềm ẩn rủi ro lún sụt nghiêm trọng tại các vùng nước sâu có nền địa chất yếu. Để giải quyết thách thức này, kết cấu nổi siêu lớn (Very Large Floating Structures - VLFS) dạng pontoon đã được đề xuất như một phương án thay thế tối ưu, thân thiện với hệ sinh thái và có khả năng triệt tiêu hoàn toàn tác động động đất truyền từ đáy biển.

Mặc dù sở hữu nhiều ưu điểm, công trình xây dựng trên nền VLFS phải đối mặt trực tiếp với điều kiện khí hậu hải văn khắc nghiệt. Trong đó, tải trọng gió thường chiếm từ 5% đến 10% tổng hợp lực môi trường tác động nhưng lại có tần suất dao động liên tục và biến đổi ngẫu nhiên, gây ra các rung động nghiêm trọng cho kết cấu thượng tầng. Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình toán học và thuật toán số hoàn chỉnh để phân tích ứng xử động lực học của hệ khung phẳng đặt trên tấm nổi siêu lớn chịu tác dụng của tải trọng gió động trên miền tần số, có xét đến sự tương tác thủy đàn hồi giữa kết cấu và môi trường nước xung quanh.

Phạm vi nghiên cứu tập trung vào mô hình khung phẳng bê tông cốt thép nhiều tầng đặt trên tấm nổi pontoon hình chữ nhật kích thước 50m x 50m, được thực hiện trong khoảng thời gian từ tháng 01/2018 đến tháng 06/2018 tại Đại học Quốc gia Thành phố Hồ Chí Minh. Ý nghĩa khoa học và thực tiễn của đề tài thể hiện qua việc cung cấp công cụ dự báo chính xác biên độ dao động, độ lệch tầng và nội lực chân cột, giúp các kỹ sư công trình biển tối ưu hóa kích thước hình học của VLFS, giảm thiểu hơn 25% rủi ro mất ổn định công trình và tiết kiệm đáng kể chi phí gia cố kết cấu.

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 dựa trên sự tích hợp của ba nền tảng lý thuyết cơ học kết cấu và cơ học thủy khí hiện đại:

Thứ nhất, lý thuyết phần tử hữu hạn (FEM) áp dụng cho hệ khung thanh phẳng 2 nút với 6 bậc tự do (gồm hai chuyển vị thẳng và một góc xoay tại mỗi nút), cho phép thiết lập ma trận độ cứng, ma trận khối lượng và véc tơ tải trọng nút do áp lực gió gây ra thông qua hàm biến đổi Fourier từ miền thời gian sang miền tần số. Đối với kết cấu tấm nổi, nghiên cứu ứng dụng lý thuyết tấm dày Mindlin (Reissner-Mindlin) thay cho lý thuyết tấm mỏng Kirchhoff cổ điển. Lý thuyết này cho phép xét đến ảnh hưởng của biến dạng trượt ngang và quán tính xoay, khắc phục hoàn toàn nhược điểm vi phạm điều kiện biên tĩnh học của tấm mỏng khi chiều dày tấm tăng lên. Mô hình sử dụng phần tử tứ giác đẳng tham số 9 nút (Q9) với hệ số hiệu chỉnh cắt bằng 5/6 để nội suy chính xác trường chuyển vị và độ cong của tấm.

Thứ hai, lý thuyết phần tử biên (BEM) kết hợp hàm thế vận tốc Free-surface Green để mô phỏng trường chất lỏng vô hạn xung quanh và bên dưới tấm nổi. Nước biển được giả thiết là chất lỏng lý tưởng, không nén được, không nhớt, chuyển động không xoáy và thỏa mãn phương trình vi phân Laplace cùng các điều kiện biên mặt thoáng, đáy biển và bức xạ sóng ở vô cực.

Thứ ba, lý thuyết tương tác thủy đàn hồi (Hydroelasticity) trong miền tần số, liên kết trực tiếp phương trình dao động của khung - tấm nổi với áp lực động học của chất lỏng thông qua ma trận khối lượng nước gia thêm và lực kích động thủy động. Hệ cản của toàn bộ kết cấu được thiết lập theo mô hình cản tỷ lệ Rayleigh với tỷ số cản quy chuẩn là 0.05 ứng với hai tần số dao động riêng đầu tiên.

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

Nghiên cứu sử dụng phương pháp mô phỏng số thông qua việc lập trình toàn bộ thuật toán liên hợp FEM-BEM trong môi trường tính toán khoa học Matlab. Nguồn dữ liệu đầu vào bao gồm các thông số cơ lý chuẩn xác: khung phẳng bê tông có mô đun đàn hồi E = 30.000.000 kN/m2, hệ số Poisson 0.2; tấm nổi có mô đun đàn hồi E = 11.900.000 kN/m2, hệ số Poisson 0.15; cùng trường vận tốc gió dao động thực nghiệm được xử lý qua chuỗi phổ năng lượng.

Về quy trình chọn mẫu và chia lưới phân tích, tác giả thực hiện khảo sát tính hội tụ trên ba cấp độ lưới phần tử tấm gồm 16 phần tử (lưới 4x4), 64 phần tử (lưới 8x8) và 256 phần tử (lưới 16x16) tích hợp các panel biên 4 nút. Kết quả cho thấy cấp lưới 64 phần tử mang lại sự cân bằng tối ưu giữa độ chính xác số học (sai số dưới 0.5%) và thời gian tính toán. Cỡ mẫu khảo sát động lực học bao gồm 7 bài toán số toàn diện, biến thiên theo độ sâu mực nước biển từ 10m đến 50m, chiều dày tấm nổi từ 1.0m đến 4.0m, chiều dài tấm từ 50m đến 150m và số tầng khung từ 1 tầng đến 8 tầng (chiều cao mỗi tầng 5.0m, nhịp khung 5.0m).

Lý do lựa chọn phương pháp phân tích liên hợp FEM-BEM là vì BEM chỉ yêu cầu rời rạc hóa bề mặt biên tiếp xúc giữa đáy tấm nổi và chất lỏng, giúp giảm số chiều không gian của bài toán thủy động từ 3D xuống 2D, từ đó cắt giảm hơn 40% khối lượng ẩn số và tài nguyên bộ nhớ so với việc chia lưới toàn bộ miền nước bằng FEM 3D thuần túy, trong khi vẫn đảm bảo sự tương thích chuyển vị hoàn hảo tại các nút liên kết khung.

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

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

Quá trình phân tích số trên 7 bài toán khảo sát đã mang lại các phát hiện khoa học quan trọng về ứng xử động của khung thượng tầng trên nền VLFS:

Thứ nhất, tính chuẩn xác và độ tin cậy của chương trình Matlab tự phát triển được kiểm chứng thành công thông qua việc so sánh với phần mềm thương mại SAP2000. Trong bài toán khung chân ngàm cứng chịu tải trọng gió, độ sai lệch về chuyển vị ngang đỉnh khung và mômen uốn chân cột giữa Matlab và SAP2000 chỉ đạt mức 0.8% đến 2.1%. Khi mở rộng sang bài toán khung đặt trên nền đàn hồi nước tĩnh, sai số tối đa vẫn duy trì ở mức dưới 3.5%, khẳng định tính chính xác tuyệt đối của thuật toán ma trận ghép nối.

Thứ hai, chiều dày tấm nổi (h) có tác động quyết định đến độ cứng tổng thể của toàn hệ thống. Khi tăng chiều dày tấm h từ 1.0m lên 3.0m, độ cứng uốn của tấm tăng gấp nhiều lần, giúp chuyển vị đỉnh khung giảm mạnh 32.4% (từ 0.045m xuống 0.030m tại tần số kích động 0.25 Hz). Đồng thời, giá trị mômen uốn tại chân cột tầng 1 giảm khoảng 18.6%, cho thấy tấm nổi dày hơn giúp phân phối lại nội lực đều hơn và giảm tải trọng cục bộ lên các liên kết khung.

Thứ ba, độ sâu mực nước biển (H) ảnh hưởng rõ rệt đến phản ứng thủy đàn hồi ở vùng nước nông nhưng có xu hướng ổn định ở vùng nước sâu. Cụ thể, khi độ sâu nước tăng từ 10m lên 30m, chuyển vị ngang đỉnh khung tăng thêm 14.2% do sự suy giảm hiệu ứng cản phản xạ từ đáy biển. Tuy nhiên, khi độ sâu vượt quá 40m, các giá trị chuyển vị và nội lực chân cột chỉ dao động nhẹ dưới 2.8%, tiệm cận về trạng thái giới hạn của môi trường nước sâu vô hạn.

Thứ tư, hiệu ứng tăng tầng của khung làm gia tăng đáng kể độ mảnh và độ nhạy cảm động lực học. Tỷ số chuyển vị đỉnh trên tổng chiều cao công trình tăng phi tuyến từ 0.12% ở khung 2 tầng lên mức 0.46% ở khung 6 tầng dưới tác động của trường gió động, tiệm cận ngưỡng giới hạn an toàn chuyển vị cho phép theo các tiêu chuẩn thiết kế công trình dân dụng hiện hành.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến các hiện tượng trên bắt nguồn từ cơ chế tương tác thủy đàn hồi phi tuyến giữa ba môi trường: kết cấu thanh, tấm nổi đàn hồi và sóng nước động. Khác với công trình trên đất liền có liên kết ngàm cứng cố định, tấm nổi VLFS đóng vai trò như một gối tựa đàn hồi biến dạng liên tục. Khi gió tác dụng lên khung, lực cắt và mômen tại chân cột truyền xuống tấm làm tấm bị uốn cục bộ, đồng thời kích thích các sóng nước bức xạ lan truyền ra xa, tạo nên lực cản thủy động và khối lượng nước gia thêm tác động ngược trở lại đáy tấm.

Khi so sánh với các công trình nghiên cứu kinh điển trên thế giới, kết quả của luận văn hoàn toàn tương đồng với các kết luận về phản ứng thủy đàn hồi của các công trình nổi được công bố bởi các chuyên gia đầu ngành trong lĩnh vực cơ học công trình biển, đồng thời phát triển sâu hơn khi chỉ ra chi tiết sự phân bố nội lực mômen và lực cắt trong từng phần tử thanh của kết cấu thượng tầng - một khía cạnh mà các nghiên cứu trước đây thường bỏ qua do xem kết cấu thượng tầng như tải trọng phân bố tĩnh.

Về phương diện trực quan hóa, toàn bộ dữ liệu phản ứng động lực học có thể được biểu diễn một cách mạch lạc thông qua hai dạng đồ thị và bảng số liệu chủ đạo:

[Biểu đồ 1: Đường cong đáp ứng tần số (FRF) biểu diễn biên độ chuyển vị đỉnh khung (m) theo dải tần số kích động từ 0.05 Hz đến 2.0 Hz, thể hiện rõ các đỉnh cộng hưởng tại vùng tần số thấp 0.15 - 0.35 Hz]

[Bảng tổng hợp: Ma trận so sánh chuyển vị lệch tầng (drift ratio) và mômen chân cột (kNm) tương ứng với các cấp chiều dày tấm h = 1m, 2m, 3m và độ sâu H = 10m, 20m, 30m, 50m]

Ý nghĩa thực tiễn của phát hiện này chỉ ra rằng việc thiết kế khung trên biển không thể áp dụng đơn thuần các sơ đồ tính tĩnh học trên nền đất cứng, mà bắt buộc phải xét đến dải tần số kích động của gió trùng với các dạng dao động riêng của hệ tương tác tấm - nước nhằm ngăn ngừa hiện tượng cộng hưởng phá hủy.

Đề 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 đề xuất 4 nhóm giải pháp kỹ thuật có tính ứng dụng cao cho công tác thiết kế và thi công công trình trên nền nổi:

Thứ nhất, tối ưu hóa tỷ lệ kích thước hình học giữa tấm nổi VLFS và khung thượng tầng. Đơn vị tư vấn thiết kế cần khống chế tỷ lệ chiều dày tấm trên chiều dài nhịp khung ở mức tối thiểu từ 0.04 đến 0.06, nhằm đảm bảo độ cứng chống uốn cục bộ và giữ cho độ lệch tầng của khung không vượt quá giới hạn 1/500 dưới cấp gió thiết kế cực hạn. Giải pháp này cần được áp dụng ngay trong giai đoạn thiết kế cơ sở với khung thời gian thực hiện 2 tháng cho mỗi dự án.

Thứ hai, tích hợp thuật toán phân tích liên hợp FEM-BEM vào quy trình thẩm tra kết cấu công trình biển. Các viện nghiên cứu chuyên ngành và cơ quan quản lý kỹ thuật cần áp dụng mô hình tính toán hydroelastic miền tần số để kiểm tra nội lực thực tế tại chân cột, giảm thiểu sai số tính toán dưới 5% so với mô hình giản lược. Kế hoạch triển khai phần mềm tính toán chuyên dụng dự kiến hoàn thành trong vòng 6 tháng.

Thứ three, lắp đặt bổ sung hệ thống giảm chấn khối lượng điều chỉnh (Tuned Mass Damper - TMD) tại tầng cao nhất của các công trình khung từ 4 tầng trở lên đặt trên VLFS. Nhà thầu thi công cơ điện và kết cấu cần phối hợp triển khai nhằm đạt mục tiêu triệt tiêu ít nhất 20% đến 25% biên độ gia tốc dao động ngang do gió giật, thực hiện lắp đặt và cân chỉnh trong thời gian 3 tháng trước khi đưa công trình vào vận hành.

Thứ tư, xây dựng bộ tiêu chuẩn kỹ thuật quốc gia về giới hạn an toàn chuyển vị và tiện nghi động lực học cho công trình dân dụng trên nền nổi. Bộ Xây dựng phối hợp cùng các trường đại học kỹ thuật chủ trì soạn thảo bộ chỉ dẫn thiết kế với ngưỡng vận tốc gió vận hành tối đa từ 28 m/s đến 35 m/s, hoàn thành phê duyệt và ban hành trong lộ trình 12 tháng.

Đố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, mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng sau:

Thứ nhất, các kỹ sư kết cấu và chuyên gia tư vấn thiết kế công trình biển. Tài liệu cung cấp cơ sở lý thuyết vững chắc và các công thức thiết lập ma trận phần tử để mô hình hóa chính xác tương tác phức tạp giữa khung và phao nổi. Trường hợp sử dụng điển hình là tính toán thiết kế các tổ hợp nhà hàng nổi, khách sạn ven biển, sàn công tác và trạm quan trắc hải văn.

Thứ hai, các nhà nghiên cứu và giảng viên chuyên ngành Kỹ thuật Xây dựng công trình thủy, Xây dựng Dân dụng & Công nghiệp. Công trình đóng vai trò như một nghiên cứu tiền đề về phương pháp số liên hợp BEM-FEM giải bài toán thủy đàn hồi 3D. Trường hợp sử dụng là phát triển các đề tài khoa học cấp bộ về đô thị nổi sinh thái, đê chắn sóng thông minh và sân bay nổi trên biển.

Thứ ba, học viên cao học và sinh viên đại học ngành Xây dựng, Công trình biển. Tài liệu là nguồn tham khảo mẫu mực về kỹ thuật lập trình số trên Matlab, ứng dụng phần tử tấm dày Mindlin Q9 và phương pháp biến đổi Fourier trong phân tích động lực học công trình. Trường hợp sử dụng là thực hiện đồ án tốt nghiệp, luận văn thạc sĩ và nghiên cứu khoa học sinh viên.

Thứ tư, các chủ đầu tư, nhà phát triển dự án bất động sản nghỉ dưỡng ven biển và năng lượng tái tạo ngoài khơi. Nghiên cứu giúp nắm bắt các thông số kỹ thuật giới hạn, đánh giá tính khả thi về mặt kết cấu và dự toán chi phí tối ưu cho nền móng nổi. Trường hợp sử dụng là lập báo cáo khả thi tiền dự án cho các trang trại điện mặt trời nổi hoặc khu đô thị lấn biển.

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

Tại sao cần kết hợp cả phương pháp phần tử hữu hạn (FEM) và phần tử biên (BEM) thay vì chỉ dùng một phương pháp?

Phương pháp FEM rất mạnh trong việc mô hình hóa các kết cấu thanh và tấm có hình học, tính chất vật liệu phức tạp, nhưng lại tốn kém tài nguyên tính toán nếu chia lưới toàn bộ miền chất lỏng 3D vô hạn. Ngược lại, BEM chỉ cần rời rạc hóa bề mặt biên tiếp xúc, tự động thỏa mãn các điều kiện biên ở vô cực thông qua hàm Green, giúp giảm đáng kể số bậc tự do và tăng độ chính xác phân tích thủy động lực học.

Lý thuyết tấm dày Mindlin mang lại ưu điểm gì vượt trội so với lý thuyết tấm mỏng Kirchhoff trong mô phỏng VLFS?

Lý thuyết tấm mỏng Kirchhoff bỏ qua hoàn toàn biến dạng trượt ngang và quán tính xoay, dẫn đến sai số lớn khi dự báo mômen uốn và lực cắt tại các vị trí biên và chân cột khung. Lý thuyết tấm dày Mindlin với phần tử tứ giác 9 nút đẳng tham số Q9 bổ sung hai bậc tự do xoay độc lập, giúp phản ánh chính xác trạng thái ứng suất thực tế của các tấm nổi có tỷ lệ chiều dày lớn trong môi trường biển.

Tải trọng gió động tác động lên công trình trên nền nổi khác biệt như thế nào so với trên đất liền?

Trên đất liền, tải trọng gió tác dụng lên công trình có liên kết chân cột ngàm cố định, phản ứng dao động chủ yếu phụ thuộc vào độ cứng của bản thân khung. Trên nền nổi VLFS, gió không chỉ gây rung lắc khung mà còn làm tấm nổi bị nghiêng và uốn, kích động sóng nước tạo ra lực phản kháng thủy đàn hồi, dẫn đến hiện tượng dao động ghép đa bậc tự do với biên độ chuyển vị phức tạp hơn nhiều.

Độ sâu của mực nước biển ảnh hưởng ra sao đến phản ứng của khung thượng tầng?

Ở các vùng nước nông dưới 20m, sóng bức xạ đáy phản xạ mạnh trở lại mặt dưới của tấm nổi, làm tăng áp lực động của nước và gia tăng chuyển vị đỉnh khung lên khoảng 14.2%. Tuy nhiên, khi độ sâu nước đạt trên 40m, ảnh hưởng của phản xạ đáy biển suy giảm hoàn toàn, đáp ứng động lực học của khung tiệm cận về giá trị không đổi của bài toán vùng nước sâu vô hạn.

Thuật toán và mã nguồn Matlab trong luận văn có thể tích hợp với các phần mềm thương mại không?

Mã nguồn Matlab được xây dựng theo cấu trúc module hóa chuẩn mực, hoàn toàn có thể xuất và nhập dữ liệu ma trận độ cứng, ma trận khối lượng với các phần mềm phân tích kết cấu thương mại như SAP2000, ANSYS hoặc ETABS. Điều này cho phép các kỹ sư dễ dàng chuyển đổi mô hình khung phức tạp từ phần mềm thương mại sang thuật toán BEM-Matlab để giải bài toán tương tác chất lỏng một cách tự động.

Kết luận

  • Luận văn đã thiết lập thành công mô hình giải tích số hoàn chỉnh phân tích tương tác thủy đàn hồi giữa khung kết cấu và tấm nổi siêu lớn (VLFS) chịu tải trọng gió động trên miền tần số.
  • Chương trình tính toán phát triển trên nền tảng Matlab đạt độ chính xác cao với sai số kiểm chứng dưới 3.5% so với phần mềm thương mại SAP2000.
  • Nghiên cứu chỉ ra quy luật định lượng về ảnh hưởng của chiều dày tấm nổi, độ sâu mực nước và số tầng khung đến chuyển vị đỉnh và nội lực chân cột, trong đó tăng chiều dày tấm h lên 3.0m giúp giảm 32.4% chuyển vị dao động đỉnh.
  • Lý thuyết tấm dày Mindlin Q9 kết hợp phương pháp phần tử biên BEM chứng minh tính ưu việt vượt trội về độ chính xác cơ học và hiệu năng tính toán so với các phương pháp số truyền thống.
  • Đóng góp khoa học cốt lõi của đề tài là cung cấp công cụ dự báo đáng tin cậy cho việc thiết kế an toàn các công trình hạ tầng đô thị trên mặt nước.

Lộ trình phát triển tiếp theo trong 12 đến 24 tháng tới bao gồm việc mở rộng thuật toán sang mô hình khung không gian 3D chịu tác động đồng thời của sóng phi tuyến và gió xoáy ngẫu nhiên. Hãy tham khảo và ứng dụng ngay các kết quả nghiên cứu này để nâng cao chất lượng và độ an toàn cho các dự án công trình biển của bạn.