Tổng quan nghiên cứu

Nhu cầu mở rộng không gian đô thị và quỹ đất ven biển tại các quốc gia có mật độ dân số cao như Nhật Bản, Hàn Quốc và Singapore đang tạo ra áp lực hạ tầng vô cùng lớn. Phương pháp quai đê lấn biển truyền thống đòi hỏi nguồn vốn đầu tư hàng tỷ USD và thường gây suy thoái tới 60% đến 80% hệ sinh thái thềm lục địa tự nhiên. Trong bối cảnh đó, kết cấu nổi siêu lớn (Very Large Floating Structure - VLFS) dạng phao (Pontoon) với chiều dài từ 500 m đến 5.000 m và chiều rộng lên đến 1.000 m nổi lên như một giải pháp bền vững, thi công module nhanh và giảm thiểu tác động môi trường biển.

Vấn đề kỹ thuật cốt lõi của công trình VLFS là phân tích ứng xử thủy đàn hồi (hydroelasticity) phức tạp khi chịu tác động đồng thời của tải trọng di động (như máy bay cất/hạ cánh hoặc dòng phương tiện giao thông) và lực nén màng trong mặt phẳng tấm (uốn dọc) do sóng biển, dòng chảy hoặc hệ thống neo gây ra. 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 của tác giả Trần Phước An (Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh, năm 2019) tập trung giải quyết bài toán này. Mục tiêu nghiên cứu là phát triển mô hình tính toán số kết hợp phương pháp phần tử biên (Boundary Element Method - BEM) và phương pháp phần tử chuyển động (Moving Element Method - MEM) trên ngôn ngữ lập trình kỹ thuật.

Nghiên cứu khảo sát một kết cấu tấm nổi có kích thước thực tế 200 m × 200 m, đặt trong vùng nước sâu 20 m, chịu tải trọng di động 2,35 kN với dải vận tốc từ 2 m/s đến 14 m/s và lực dọc biến thiên từ 0 đến 90 kN/m. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn cao, giúp giảm hơn 40% khối lượng tính toán ma trận so với phương pháp phần tử hữu hạn truyền thống, đồng thời cung cấp cơ sở định lượng quan trọng để kiểm soát độ võng và an toàn chịu lực cho sân bay nổi, cầu phao và kho dầu 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 được xây dựng trên hệ thống lý thuyết liên ngành vững chắc giữa cơ học kết cấu tấm và thủy động lực học công trình biển:

Lý thuyết tấm mỏng Kirchhoff có xét ảnh hưởng của lực màng (uốn dọc): Tấm nổi được mô hình hóa theo giả thiết pháp tuyến thẳng không đổi độ dài và vuông góc với mặt trung hòa sau khi biến dạng. Phương trình vi phân đạo hàm riêng bậc 4 mô tả dao động uốn của tấm có dạng liên kết giữa độ cứng chống uốn $D$, tải trọng ngang $p$ và các thành phần lực trong mặt phẳng gồm lực dọc $N_x, N_y$ cùng lực trượt $N_{xy}$.

Lý thuyết sóng biển tuyến tính (Linear Potential Wave Theory): Môi trường chất lỏng xung quanh được giả định là lý tưởng, không nén được, không nhớt và chuyển động không xoáy. Hàm thế vận tốc $\Phi(x, y, z, t)$ thỏa mãn phương trình vi phân Laplace trong toàn bộ miền chất lỏng, kết hợp các điều kiện biên động lực học tại đáy biển sâu $H = 20\text{ m}$, biên mặt thoáng tự do và mặt tiếp xúc đáy kết cấu.

Hệ tọa độ tương đối gắn liền với tải trọng chuyển động: Điểm mấu chốt của phương pháp MEM là chuyển đổi hệ trục tọa độ cố định sang hệ tọa độ chuyển động $(r = x - vt, s = y, z)$ với vận tốc không đổi $V$. Việc chuyển đổi này biến bài toán tải trọng thay đổi vị trí theo thời gian thành tải trọng cố định tại gốc tọa độ tương đối, triệt tiêu nhu cầu tái phân chia lưới phần tử qua từng bước thời gian.

Các khái niệm chính trong mô hình bao gồm: cấu trúc nổi siêu lớn dạng phao (Pontoon VLFS), ma trận tương tác thủy động lực, lớp cản nhân tạo hấp thụ sóng phản xạ ở biên vô hạn và phương pháp tích phân từng bước Newmark.

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

Mô hình tính toán số được thiết lập với bộ thông số kỹ thuật tiêu chuẩn của tấm nổi:

  • Chiều dài tấm $L = 200\text{ m}$, chiều rộng $B = 200\text{ m}$, chiều dày ban đầu $h_s = 0,17\text{ m}$.
  • Mô đun đàn hồi của vật liệu $E = 5 \times 10^5\text{ kN/m}^2$, hệ số Poisson $\nu = 0,13$, khối lượng riêng tương đương $m = 0,25625\text{ T/m}^3$.
  • Độ sâu mực nước biển $h = 20\text{ m}$, tải trọng di động $P_0 = 2,35\text{ kN}$.
  • Dải lực dọc khảo sát $P_x, P_y$ chạy từ 0 đến $90\text{ kN/m}$ với bước nhảy $10\text{ kN/m}$.
  • Dải vận tốc khảo sát: $V = 2, 4, 6, 8, 10, 12, 14\text{ m/s}$.

Mô hình sử dụng phần tử tấm tứ giác 4 nút (Q4) với mỗi nút chứa 3 bậc tự do gồm chuyển vị đứng $w$ và hai góc xoay $\theta_x, \theta_y$, tạo nên 12 bậc tự do cho mỗi phần tử. Hàm nội suy Hermite bậc cao được ứng dụng để đảm bảo tính liên tục của độ dốc chuyển vị. Miền tấm được chia lưới thành các phần tử kích thước đồng nhất $2\text{ m} \times 4\text{ m}$. Phương pháp tích phân số cầu phương Gauss được áp dụng để tính toán chính xác ma trận độ cứng, ma trận cản và ma trận khối lượng.

Quy trình giải hệ phương trình vi phân chuyển động phi tuyến kết hợp tấm - chất lỏng được thực hiện bằng phương pháp Newmark dạng chuyển vị và gia tốc trong miền thời gian từ $t = 0\text{ s}$ đến $t = 20\text{ s}$. Phương pháp BEM-MEM được lựa chọn vì BEM giảm số chiều không gian của miền chất lỏng vô hạn xuống bề mặt biên, trong khi MEM loại bỏ hoàn toàn sai số tích lũy khi cập nhật vị trí tải trọng của phương pháp FEM thông thường.

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 8 bài toán khảo sát toàn diện đã đưa ra các kết quả cơ học định lượng rõ nét:

Thứ nhất, phân hóa phản ứng động theo ngưỡng vận tốc tới hạn ($C_{\min} \approx 6 - 8\text{ m/s}$): Khi vận tốc tải trọng nhỏ hơn vận tốc sóng tới hạn ($V = 2, 4, 6\text{ m/s}$), hình dạng chuyển vị của tấm nổi có tính đối xứng cục bộ quanh điểm đặt tải. Ngược lại, khi vận tốc vượt ngưỡng tới hạn ($V = 8, 10, 12, 14\text{ m/s}$), ứng xử thủy đàn hồi chuyển sang trạng thái bất đối xứng rõ rệt, hình thành sóng lan truyền kéo dài phía sau vệt di chuyển của tải trọng.

Thứ hai, tác động khuếch đại chuyển vị của lực nén màng: Khi lực dọc theo phương chuyển động $P_x$ tăng từ $0\text{ kN/m}$ lên $90\text{ kN/m}$, độ võng lớn nhất của tấm nổi tăng đều đặn từ 12% đến 18% ở vùng vận tốc dưới tới hạn. Hiện tượng này bắt nguồn từ sự suy giảm độ cứng uốn hiệu dụng của tấm dưới tác động của ứng suất nén trong mặt phẳng.

Thứ ba, sự vượt trội của ảnh hưởng lực dọc dọc trục ($P_x$) so với lực dọc ngang trục ($P_y$): Lực dọc $P_x$ cùng phương với quỹ đạo chuyển động của tải trọng tạo ra độ gia tăng chuyển vị cực đại lớn hơn từ 15% đến 20% so với trường hợp tấm chịu lực dọc ngang $P_y$ có cùng cường độ $90\text{ kN/m}$.

Thứ tư, hiệu ứng giảm võng khi tăng chiều dày kết cấu: Khi tăng chiều dày tấm $h_s$ từ $0,17\text{ m}$ lên các mức dày hơn trong bài toán khảo sát số 8, độ võng động cực đại của tấm giảm hơn 45%, chứng minh vai trò chi phối của độ cứng uốn hình học đối với khả năng chịu tải di động.

Thảo luận kết quả

Cơ chế vật lý của hiện tượng bất đối xứng ở vận tốc cao ($V > C_{\min}$) được giải thích qua sự tương tác năng lượng giữa tải trọng di chuyển và sóng mặt nước. Khi tải trọng di chuyển nhanh hơn tốc độ truyền sóng tự do trên bề mặt tấm nổi, năng lượng dao động không kịp phân tán về phía trước mà bị nén lại và bức xạ thành chuỗi sóng thủy đàn hồi về phía sau, tương tự như hiện tượng sóng xung kích Mach trong khí động lực học.

So sánh với các nghiên cứu sử dụng phương pháp BEM-FEM truyền thống, thuật toán BEM-MEM trong luận văn duy trì được tính hội tụ ổn định ngay cả khi tải trọng di chuyển đến gần các cạnh biên $L = 200\text{ m}$. Sai số số học do nội suy vị trí tải giữa các bước thời gian được triệt tiêu hoàn toàn.

Dữ liệu tính toán trong luận văn được minh họa trực quan qua 46 đồ thị và 8 bảng số liệu chi tiết. Các biểu đồ mặt cắt chuyển vị dọc theo trục $x$ và trục $y$ làm nổi bật vị trí đáy phễu võng lớn nhất lệch khỏi tâm tải trọng khi $V \ge 8\text{ m/s}$. Các bảng so sánh chuyển vị cực đại theo từng cấp lực nén $0, 10, 20, \dots, 90\text{ kN/m}$ cung cấp cái nhìn định lượng chuẩn xác phục vụ thiết kế thực tế.

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

Dựa trên các kết quả mô phỏng số, tác giả đề xuất 4 nhóm giải pháp kỹ thuật cụ thể:

  1. Tích hợp module thuật toán BEM-MEM vào các phần mềm mô phỏng kết cấu công trình biển chuyên dụng. Nhóm kỹ sư phát triển phần mềm và các viện nghiên cứu kết cấu thủy công cần triển khai trong vòng 6 đến 12 tháng, hướng tới mục tiêu giảm ít nhất 35% thời gian xử lý bài toán thủy đàn hồi 3D quy mô lớn.

  2. Thiết lập giới hạn vận tốc khai thác an toàn cho phương tiện và máy bay trên mặt sàn tấm nổi. Đơn vị quản lý vận hành sân bay nổi hoặc cầu phao nổi cần quy định tốc độ di chuyển của phương tiện không vượt quá 70% vận tốc tới hạn $C_{\min}$ (tương đương duy trì $V \le 5\text{ m/s}$ trong điều kiện mô phỏng tiêu chuẩn) nhằm kiểm soát độ võng động trong giới hạn an toàn cho phép $L/500$.

  3. Tối ưu hóa hệ thống neo định vị để triệt tiêu lực nén màng $P_x$ dọc theo hướng cất hạ cánh chính. Các kỹ sư tư vấn thiết kế công trình ngoài khơi cần bố trí dây neo phân bổ ứng suất căng đều, đảm bảo ứng suất nén dư trong tấm không vượt quá $30\text{ kN/m}$ trong giai đoạn lập dự án 12 tháng.

  4. Xây dựng bộ tiêu chuẩn kỹ thuật thiết kế tấm nổi siêu lớn chịu tải trọng động lực học. Bộ Xây dựng phối hợp cùng Bộ Giao thông Vận tải và các trường đại học khối kỹ thuật cần xây dựng quy chuẩn đánh giá ứng xử uốn dọc cho công trình nổi trong lộ trình 24 đến 36 tháng tới.

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

Nội dung nghiên cứu mang giá trị học thuật và ứng dụng chuyên sâu cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế công trình biển và công trình thủy: Nắm bắt phương pháp tính toán tương tác thủy đàn hồi chính xác để lựa chọn chiều dày tấm, vật liệu composite hoặc bê tông nhẹ và bố trí hệ thống neo chống uốn dọc cho sân bay nổi, bến cảng container nổi.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Xây dựng - Thủy lợi: Khai thác khung lý thuyết Kirchhoff cải biên, phương pháp phần tử chuyển động MEM và code tính toán Matlab làm tài liệu tham khảo để mở rộng sang các mô hình tấm Mindlin, tấm vật liệu có cơ tính biến thiên (FGM).
  • Đơn vị quy hoạch hạ tầng giao thông và đô thị biển: Tiếp cận giải pháp công trình nổi siêu lớn VLFS như một phương án thay thế lấn biển thân thiện môi trường, phục vụ quy hoạch đảo nhân tạo và cầu vượt biển nổi.
  • Chuyên gia phân tích số và phát triển phần mềm CAE: Tham khảo kỹ thuật xây dựng ma trận phần tử tứ giác 4 nút Hermite kết hợp tích phân Newmark và lớp cản hấp thụ sóng để lập trình các gói giải thuật động lực học tải trọng di động.

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

Phương pháp BEM-MEM có điểm gì ưu việt hơn so với phương pháp FEM truyền thống khi phân tích tấm nổi chịu tải di động? Phương pháp MEM thiết lập hệ tọa độ tương đối di chuyển đồng tốc với tải trọng, biến tải trọng di động thành tải trọng tĩnh tại gốc tọa độ. Điều này giúp loại bỏ hoàn toàn việc liên tục phân chia lại lưới phần tử và cập nhật vectơ tải sau mỗi bước thời gian $t$, tiết kiệm trên 40% chi phí tính toán và tránh hiện tượng mất ổn định số tại biên kết cấu $200\text{ m}$.

Lực nén màng (lực uốn dọc) ảnh hưởng thế nào đến độ võng của kết cấu nổi siêu lớn? Khi lực nén màng $P_x, P_y$ tăng từ $0\text{ kN/m}$ lên $90\text{ kN/m}$, độ võng uốn của tấm tăng tuyến tính từ 12% đến 18%. Lực nén tạo ra hiệu ứng làm mềm hình học, làm giảm độ cứng chống uốn tổng thể của tấm mỏng Kirchhoff, khiến kết cấu dễ bị lún sâu hơn dưới cùng một giá trị tải trọng động $2,35\text{ kN}$.

Ngưỡng vận tốc tới hạn của tải trọng di động trên tấm nổi mang ý nghĩa gì trong thực tế? Vận tốc tới hạn $C_{\min}$ (khoảng $6 - 8\text{ m/s}$ trong bài toán khảo sát) là ranh giới thay đổi hoàn toàn đặc trưng dao động của tấm. Ở dải vận tốc $V > C_{\min}$, chuyển vị không còn đối xứng mà hình thành các dải sóng thủy đàn hồi kéo dài phía sau vệt di chuyển, gây rung lắc mạnh và đòi hỏi kiểm soát an toàn nghiêm ngặt khi máy bay hạ cánh.

Tại sao phần tử tứ giác 4 nút Q4 với hàm nội suy Hermite lại được lựa chọn trong nghiên cứu? Hàm nội suy Hermite bậc cao chứa các đạo hàm riêng tương ứng với các góc xoay tại nút, đảm bảo tính liên tục của độ võng $w$ và độ dốc mặt trung hòa giữa các phần tử $2\text{ m} \times 4\text{ m}$. Nhờ đó, ma trận độ cứng phần tử phản ánh chính xác thế năng biến dạng uốn và mô men xoắn theo lý thuyết tấm mỏng Kirchhoff.

Mô hình tính toán trong luận văn có thể mở rộng cho những dạng bài toán kết cấu nào khác? Thuật toán BEM-MEM có thể mở rộng linh hoạt cho bài toán tấm dày Mindlin, tấm composite nhiều lớp trên nền đàn hồi Pasternak hoặc nền đàn nhớt. Ngoài ra, mô hình hoàn toàn có khả năng phân tích chuỗi nhiều tải trọng di động phức tạp như đoàn tàu cao tốc hoặc đoàn xe tải nặng di chuyển liên tục trên cầu nổi.

Kết luận

  • Xây dựng thành công thuật toán kết hợp BEM-MEM trên hệ tọa độ chuyển động, giải quyết triệt để bài toán tương tác thủy đàn hồi giữa tấm mỏng Kirchhoff và chất lỏng lý tưởng dưới tải trọng di động $2,35\text{ kN}$.
  • Đánh giá định lượng tác động của lực nén màng $P_x, P_y$ từ $0\text{ kN/m}$ đến $90\text{ kN/m}$, chỉ ra rằng ứng suất nén làm gia tăng độ võng cực đại của tấm lên tới 18%, trong đó lực dọc trục $P_x$ có mức độ ảnh hưởng lớn hơn 15% đến 20% so với lực ngang $P_y$.
  • Xác định rõ sự chuyển đổi trạng thái cơ học từ đối xứng sang bất đối xứng hình thành chuỗi sóng kéo theo khi vận tốc tải trọng vượt ngưỡng tới hạn $V > 6\text{ m/s}$ trong miền thời gian $0 \to 20\text{ s}$.
  • Chứng minh tính ưu việt của việc tăng chiều dày tấm $h_s$ giúp triệt giảm hơn 45% độ võng động, cung cấp thông số đầu vào thiết yếu cho bài toán tối ưu hóa hình học kết cấu phao nổi 200 m × 200 m.
  • Khẳng định độ tin cậy và khả năng ứng dụng thực tế cao của chương trình tính toán số viết trên Matlab, mở ra hướng đi bền vững cho các dự án công trình nổi siêu lớn thay thế giải pháp lấn biển truyền thống.

Trong lộ trình nghiên cứu tiếp theo từ 12 đến 24 tháng tới, mô hình tính toán cần được tích hợp thêm các yếu tố phi tuyến của sóng biển ngẫu nhiên và ảnh hưởng của hệ thống neo đàn hồi. Các kỹ sư và nhà nghiên cứu công trình biển quan tâm có thể khai thác trực tiếp cơ sở dữ liệu và thuật toán số này để nâng cao độ an toàn cho các công trình hạ tầng nổi ngoài khơi thế hệ mới.