Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình ngoài khơi, các kết cấu như giàn khoan tự nâng và trụ tuabin điện gió thường xuyên phải chịu tải trọng môi trường khổng lồ với lực sóng và gió dao động từ 4 MN đến hơn 25 MN, cùng tải trọng bản thân lên tới 200 MN tại các vùng biển có độ sâu từ 10m đến 30m. Vấn đề cốt lõi đặt ra là các phương pháp phân tích truyền thống thường lý tưởng hóa liên kết chân móng thành ngàm cứng hoặc các gối tựa lò xo đàn hồi tuyến tính đơn giản. Cách tiếp cận này bỏ qua hoàn toàn hiện tượng biến dạng dẻo tích lũy, sự suy giảm độ cứng của nền đất cát và đất sét dưới tác động của tải trọng lặp tuần hoàn, dẫn đến việc đánh giá sai lệch đáng kể phản ứng động học và nội lực thực tế của kết cấu khung bên trên.

Mục tiêu cụ thể của luận văn là xây dựng quy trình và thuật toán tương tác hai chiều hoàn chỉnh giữa mô hình móng phi tuyến đàn dẻo ISIS và mô hình phân tích khung kết cấu không gian 3 chiều SOS_3D. Nghiên cứu tập trung giải quyết bài toán tương tác đất - kết cấu cho giàn khoan 3 chân ngoài khơi sử dụng móng spudcan và móng caisson chịu tác động của hệ tải trọng kết hợp 6 bậc tự do.

Nghiên cứu được hoàn thành tại Trường Đại học Bách Khoa, Đại học Quốc gia Thành phố Hồ Chí Minh. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi giúp tối ưu hóa từ 15% đến 25% độ an toàn kết cấu khung ngoài khơi, đồng thời cắt giảm hơn 30% thời gian tính toán so với các mô hình phần tử hữu hạn 3D khối đất liên tục thông qua việc áp dụng khái niệm phần tử vĩ mô hiệ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 dựa trên sự tích hợp của hai nền tảng lý thuyết cơ học tiên tiến:

Thứ nhất là lý thuyết chảy dẻo nhiều mặt cho nền móng công trình ngoài khơi do Houlsby và Puzrin khởi xướng năm 2000, sau đó được phát triển hoàn thiện bởi Nguyễn Sỹ Lâm năm 2005 trong chương trình ISIS. Mô hình này sử dụng hệ nhiều mặt chảy dẻo lồng nhau cùng quy luật củng cố hỗn hợp đẳng hướng và động học để mô tả chính xác ứng xử trễ đàn - dẻo và sự tiêu tán năng lượng của đất dưới tải trọng chu kỳ.

Thứ hai là lý thuyết phần tử hữu hạn dầm - cột không gian 3 chiều phi tuyến hình học xét đến hiệu ứng P-Delta và ma trận định hướng khớp quay hữu hạn theo nghiên cứu của Kassimali và Abbasnia năm 1991, được tích hợp trong chương trình SOS_3D.

Mô hình nghiên cứu bao gồm bốn khái niệm then chốt:

  • Phần tử vĩ mô (macro-element): Gom toàn bộ ứng xử phức tạp của nền đất vào một điểm tham chiếu tải trọng duy nhất tại đáy móng với 6 thành phần lực và chuyển vị tương ứng.
  • Mặt chảy dẻo chuẩn hóa: Biên giới hạn đàn hồi trong không gian lực 6 chiều với thông số lệch tâm elip a bằng -0.2 cho cát thô silic và hệ số ma sát chuẩn hóa h0 bằng 0.154.
  • Quy luật chảy dẻo không kết hợp: Cho phép mô phỏng chính xác sự biến đổi thể tích, hiện tượng trượt và chuyển vị đứng của móng dưới tác dụng của tải trọng ngang và momen uốn.
  • Ma trận cản Rayleigh: Kết hợp tuyến tính giữa ma trận khối lượng và ma trận độ cứng để lọc bỏ các dao động bậc cao trong phân tích động lực học kết cấu.

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

Nghiên cứu sử dụng nguồn dữ liệu thực nghiệm tiêu chuẩn từ Trung tâm Nền móng Công trình Ngoài khơi thuộc Đại học Tây Úc và các chuỗi thí nghiệm chịu tải phức tạp của Đại học Oxford trên đất sét Speswhite và cát khô Yellow Leighton-Buzzard.

Phương pháp phân tích chủ đạo là phương pháp phần tử hữu hạn phi tuyến kết hợp thuật toán lặp Newton-Raphson từng bước gia tải. Trong mỗi bước lặp, ma trận độ cứng tiếp tuyến từ mô-đun móng ISIS được cập nhật trực tiếp vào ma trận độ cứng tổng thể của mô-đun kết cấu SOS_3D thông qua mã nguồn được lập trình bằng ngôn ngữ FORTRAN.

Cỡ mẫu nghiên cứu được thiết lập cụ thể trên mô hình giàn khoan tự nâng 3 chân điển hình với móng spudcan hình nón tương đương có đường kính tiếp xúc đáy biển 15m. Phương pháp chọn mẫu điển hình này được áp dụng vì giàn khoan tự nâng 3 chân chiếm hơn 70% tổng số lượng công trình thăm dò dầu khí tại các vùng biển nông và trung bình từ 10m đến 30m. Phương pháp phân tích phần tử vĩ mô được lựa chọn vì vượt trội hơn hẳn mô hình phần tử hữu hạn khối đất liên tục về tốc độ xử lý, loại bỏ việc chia lưới hàng vạn phần tử đất phức tạp mà vẫn giữ được độ chính xác kỹ thuật vượt trội. Quá trình tính toán mô phỏng và kiểm chứng thuật toán được thực hiện liên tục từ tháng 07/2011 đến tháng 02/2013.

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

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

Phân tích số cho thấy sự khác biệt rõ rệt giữa các mô hình tính toán, thể hiện qua 3 phát hiện kỹ thuật cốt lõi:

Thứ nhất, việc xét đến tính mềm và biến dạng dẻo của móng làm tái phân phối đáng kể nội lực trong kết cấu khung. So với mô hình ngàm cứng lý tưởng, mô hình có xét tương tác móng làm giảm momen uốn tại vị trí chân cột từ 22% đến 35%. Tuy nhiên, chuyển vị ngang tại đỉnh giàn khoan lại tăng lên từ 18% đến 28% do liên kết chân móng bị xoay và trượt dẻo.

Thứ hai, mô hình 5 mặt chảy dẻo thể hiện tính ưu việt vượt trội trong việc mô tả hiện tượng suy giảm độ cứng liên tục và hình thành các vòng trễ năng lượng khép kín dưới tải trọng đổi chiều. Ngược lại, mô hình một mặt chảy dẻo xem toàn bộ miền bên trong mặt giới hạn là đàn hồi thuần túy, dẫn đến việc bỏ qua các biến dạng dẻo tích lũy nhỏ và gây sai lệch độ lún dư tích lũy lên tới hơn 40% sau nhiều chu kỳ tải.

Thứ ba, sức chịu tải ngang và khả năng kháng momen của móng spudcan đạt giá trị cực đại khi tải trọng đứng duy trì ở mức khoảng 50% sức chịu tải đứng giới hạn lớn nhất. Khi tải trọng đứng vượt quá ngưỡng này hoặc giảm xuống quá thấp, biên độ mặt chảy dẻo bị co hẹp, làm tăng nguy cơ trượt phẳng hoặc bật nhổ móng.

Thảo luận kết quả

Nguyên nhân chính của sự suy giảm momen chân cột là do đất nền dưới đáy móng spudcan bị phá hủy dẻo cục bộ và chuyển vị xoay, hoạt động như một khớp dẻo bán cứng giúp giải phóng năng lượng tập trung. Kết quả này hoàn toàn phù hợp với các quan sát thực nghiệm của Martin năm 1994 và Cassidy năm 1999 tại hiện trường Biển Bắc, đồng thời khẳng định rằng các thiết kế giả định ngàm cứng đang đánh giá quá cao momen uốn chân cột nhưng lại đánh giá quá thấp độ võng tổng thể của giàn.

Trong báo cáo kỹ thuật, dữ liệu tương tác được trực quan hóa sinh động thông qua các đồ thị quan hệ lực ngang - chuyển vị ngang và momen uốn - góc xoay. Các biểu đồ vòng trễ năng lượng hình thoi thể hiện rõ khả năng dập tắt dao động của nền đất. Bảng so sánh nội lực giữa 3 chân móng cho thấy chân móng phía đón gió chịu tải trọng nhổ và lực cắt gia tăng mạnh, trong khi hai chân móng phía khuất gió chịu tải trọng nén tăng thêm khoảng 30%, chứng minh tính bất đối xứng rõ rệt trong làm việc không gian của kết cấu.

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

Nhằm nâng cao chất lượng thiết kế và an toàn vận hành cho các công trình biển, nghiên cứu đưa ra 4 khuyến nghị hành động thiết thực:

  1. Áp dụng quy trình phân tích tương tác phi tuyến khung - móng vào công tác thiết kế thực tế cho 100% các dự án giàn khoan tự nâng và tháp tuabin gió móng caisson từ nay đến năm 2026. Chủ thể thực hiện là các kỹ sư trưởng và công ty tư vấn thiết kế công trình dầu khí biển, với mục tiêu cắt giảm từ 8% đến 12% khối lượng thép kết cấu chân giàn nhờ tận dụng sự phân phối lại nội lực.
  2. Tích hợp thuật toán phần tử vĩ mô ISIS vào các bộ phần mềm thương mại phân tích kết cấu biển phổ biến trong giai đoạn 12 đến 18 tháng tới. Nhóm nghiên cứu tại các trường đại học kỹ thuật và viện dầu khí chịu trách nhiệm chuyển giao công nghệ, nhằm rút ngắn thời gian tính toán động đất và sóng bão xuống từ 3 đến 5 lần so với phương pháp chia lưới đất thông thường.
  3. Thực hiện đo đạc và quan trắc lún thực địa định kỳ trên ít nhất 5 giàn khoan đang vận hành tại thềm lục địa Việt Nam trong vòng 2 năm. Chủ đầu tư và các đơn vị vận hành mỏ cần tiến hành thu thập số liệu góc xoay chân móng sau mỗi mùa mưa bão để hiệu chỉnh các tham số độ dẻo của đất cát và đất sét tại chỗ.
  4. Cập nhật và bổ sung tiêu chuẩn kỹ thuật thiết kế công trình ngoài khơi quốc gia trước năm 2027 do Bộ Xây dựng phối hợp Bộ Công Thương ban hành, trong đó quy định bắt buộc phải kiểm tra ứng xử dẻo tích lũy của móng spudcan dưới tác động của bão biển cấp 12 trở lên.

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

Nội dung luận văn mang lại giá trị chuyên sâu và công cụ thực hành cụ thể cho 4 nhóm đối tượng:

  • Kỹ sư thiết kế kết cấu công trình biển: Nắm vững phương pháp mô hình hóa chính xác liên kết chân móng, áp dụng trực tiếp để tối ưu hóa kích thước tiết diện dầm cột cho hơn 10 loại hình giàn khoan và giàn khai thác cố định.
  • Kỹ sư địa kỹ thuật ngoài khơi: Sử dụng khung lý thuyết chảy dẻo 6 bậc tự do và quy luật củng cố hỗn hợp để dự báo chính xác sức chịu tải, độ lún tích lũy và chuyển vị ngang của móng spudcan và suction caisson trong nền đất yếu.
  • Giảng viên và nghiên cứu sinh chuyên ngành Xây dựng Dân dụng - Công nghiệp và Công trình Biển: Khai thác mã nguồn thuật toán FORTRAN, hệ phương trình ma trận độ cứng phi tuyến và các lưu đồ tính toán làm tài liệu tham khảo cho các đề tài nghiên cứu chuyên sâu về cơ học đất phi tuyến.
  • Quản lý dự án và chủ đầu tư năng lượng ngoài khơi: Hiểu rõ cơ chế làm việc thực tế của móng công trình biển để đưa ra quyết định lựa chọn phương án móng tối ưu, giúp tiết kiệm từ 10% đến 15% tổng chi phí đầu tư xây lắp nền móng.

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

Mô hình phần tử vĩ mô có ưu điểm vượt trội gì so với phương pháp phần tử hữu hạn khối đất thông thường? Mô hình phần tử vĩ mô đóng gói toàn bộ ứng xử phi tuyến của đất nền vào một điểm đại diện duy nhất tại đáy móng với 6 bậc tự do. Cách làm này giúp giảm hơn 90% số bậc tự do của toàn hệ thống, rút ngắn thời gian xử lý từ vài ngày xuống còn vài phút mà vẫn đảm bảo độ chính xác trên 95% trong việc mô phỏng phản lực và chuyển vị chân móng.

Tại sao mô hình nhiều mặt chảy dẻo lại chính xác hơn mô hình một mặt chảy dẻo truyền thống? Mô hình một mặt chảy dẻo giả định toàn bộ không gian bên trong mặt giới hạn hoàn toàn đàn hồi tuyến tính, do đó không thể mô tả được độ suy giảm độ cứng và sự tích lũy biến dạng dẻo dưới tải trọng đổi chiều. Mô hình 5 mặt chảy dẻo phân chia miền ứng suất thành nhiều tầng giới hạn lồng nhau, giúp tái hiện chính xác vòng trễ năng lượng và độ lún dẻo tích lũy qua từng chu kỳ sóng.

Quy trình kết nối giữa chương trình ISIS và SOS_3D được triển khai như thế nào? Hai chương trình được liên kết chặt chẽ thông qua giao thức truyền dữ liệu lặp từng bước được viết bằng ngôn ngữ FORTRAN. Trong mỗi bước tăng tải, chuyển vị tại các nút chân cột từ SOS_3D được truyền sang ISIS để tính toán trạng thái dẻo của móng, sau đó ISIS cập nhật ma trận độ cứng tiếp tuyến ngược trở lại ma trận độ cứng tổng thể của SOS_3D cho đến khi hội tụ nghiệm Newton-Raphson.

Nghiên cứu có thể áp dụng cho loại móng công trình biển nào khác ngoài móng spudcan? Toàn bộ khung lý thuyết và thuật toán xây dựng trong luận văn được thiết lập tổng quát cho bài toán 6 bậc tự do, do đó hoàn toàn áp dụng hiệu quả cho móng giếng chìm hạ bằng phương pháp hút nước (suction caisson) và móng nông hình tròn phẳng dùng cho trụ tuabin gió ngoài khơi trên cả nền cát và nền sét.

Hiệu ứng phi tuyến hình học P-Delta ảnh hưởng như thế nào đến độ ổn định của giàn khoan? Khi giàn khoan bị nghiêng do tải trọng gió và sóng, tải trọng đứng rất lớn tác dụng lệch tâm sẽ sinh ra thêm momen uốn phụ P-Delta dọc thân cột, làm tăng momen chân cột thêm từ 15% đến 25%. Sự gia tăng này đẩy nhanh điểm trạng thái ứng suất tại đáy móng chạm tới mặt chảy dẻo, làm giảm khả năng chịu tải tổng thể của công trình.

Kết luận

  • Luận văn đã xây dựng thành công thủ tục tương tác hoàn chỉnh giữa mô hình móng phi tuyến ISIS và mô hình kết cấu khung không gian SOS_3D bằng ngôn ngữ FORTRAN.
  • Ứng dụng thành công lý thuyết chảy dẻo nhiều mặt với hệ 5 mặt giới hạn để mô phỏng chính xác hiện tượng suy giảm độ cứng và vòng trễ năng lượng của nền đất dưới tải trọng tuần hoàn.
  • Chứng minh rõ ràng quy luật phân phối lại nội lực khung: giảm momen uốn chân cột tới 35% nhưng làm tăng chuyển vị ngang đỉnh giàn thêm 28% so với mô hình ngàm cứng truyền thống.
  • Thiết lập công cụ tính toán hiệu quả cao cho móng spudcan và móng caisson chịu tác động của hệ lực không gian phức tạp với đầy đủ 6 bậc tự do.
  • Đóng góp luận cứ khoa học quan trọng giúp nâng cao tính an toàn và tiết kiệm đáng kể chi phí vật liệu trong thiết kế công trình dầu khí và năng lượng tái tạo ngoài khơi.

Trong vòng 12 tháng tới, hướng nghiên cứu tiếp theo sẽ tập trung mở rộng mô hình cho bài toán phân tích động lực học ngẫu nhiên theo miền thời gian dưới tác động đồng thời của động đất và sóng bão cực đại. Các đơn vị tư vấn và kỹ sư công trình biển hãy áp dụng ngay phương pháp tương tác khung - móng tiên tiến này vào các dự án thực tế để tối ưu hóa thiết kế và đảm bảo an toàn tuyệt đối cho công trình biển.