Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình bến cảng, móng cọc đóng vai trò cốt lõi khi tiếp nhận toàn bộ tải trọng thẳng đứng từ bản mặt cầu tàu cùng tải trọng ngang khổng lồ do tác động của sóng, gió, dòng chảy và lực va neo tàu. Thống kê thực tế cho thấy các tổ hợp lực ngang trong công trình cảng thường chiếm từ 60% đến 70% các điều kiện kiểm toán nguy hiểm nhất đối với kết cấu móng. Đặc biệt, tại các khu vực địa chất ven biển miền Nam Việt Nam, tầng đất mặt phần lớn là bùn sét mềm yếu bão hòa nước có chiều dày từ 12m đến 18m, khiến chuyển vị ngang của cọc trở thành bài toán kỹ thuật phức tạp hàng đầu.

Vấn đề nghiên cứu trọng tâm của đề tài là phân tích ứng xử biến dạng và nội lực của cọc đơn chịu tải trọng ngang trong môi trường đất yếu công trình thủy, từ đó đánh giá sự phù hợp của các mô hình tính toán tương tác cọc - đất - công trình hiện nay. Mục tiêu cụ thể là xây dựng quy trình mô phỏng chính xác chuyển vị ngang và mômen uốn dọc thân cọc, tiến hành hiệu chỉnh đường cong phản lực nền phi tuyến để áp dụng vào các phần mềm kết cấu đa tổ hợp tải trọng.

Phạm vi nghiên cứu được thực hiện dựa trên số liệu địa chất và kết quả thí nghiệm nén tĩnh ngang hiện trường tại dự án Kho - Cảng chứa LPG lạnh Thị Vải thuộc tỉnh Bà Rịa - Vũng Tàu trong năm 2013. Kết quả nghiên cứu có ý nghĩa thực tiễn to lớn khi giúp giảm sai số dự báo chuyển vị đầu cọc từ 35% xuống dưới 6%, đồng thời hỗ trợ tối ưu hóa chiều sâu cắm cọc, giúp tiết kiệm từ 15% đến 20% chi phí xây lắp móng công trình bến mà vẫn đảm bảo tuyệt đối hệ số an toàn vận hành.

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 của ba lý thuyết cơ học đất và kết cấu cọc chịu tải ngang kinh điển:

  • Lý thuyết dầm đàn hồi trên nền đàn hồi của Winkler (1867) và Hetenyi (1946): Mô hình hóa tương tác giữa thân cọc và môi trường đất xung quanh thông qua các gối đàn hồi độc lập. Phương trình vi phân bậc bốn mô tả trạng thái cân bằng uốn của cọc chịu lực dọc và phản lực nền theo phương ngang được thiết lập với độ cứng lò xo tỷ lệ thuận với mô đun phản lực nền.
  • Lý thuyết giới hạn phá hoại của Broms (1964): Phân tích sức chịu tải ngang tới hạn của cọc ngắn cứng và cọc dài mềm trong hai môi trường đất thuần dính và đất rời thuần túy. Broms xác định khả năng chịu lực giới hạn dựa trên sự hình thành khớp dẻo và mômen chảy của vật liệu chế tạo cọc thép hoặc bê tông cốt thép.
  • Phương pháp đường cong phi tuyến p-y của Reese, Matlock (1956) và Cox (1974): Thay thế phản lực nền tuyến tính bằng hàm phi tuyến thực nghiệm mô tả quan hệ giữa áp lực đất p (kN/m) và chuyển vị ngang y (m) tại từng cao độ chiều sâu.

Nghiên cứu tập trung vào năm khái niệm cơ bản:

  • Hệ số phản lực nền theo phương ngang kh (kN/m3) và mô đun biến dạng E50 của đất sét yếu dao động từ 2000 kPa đến 4500 kPa.
  • Sức kháng cắt không thoát nước Cu của đất sét bão hòa nước đạt dưới 25 kPa tại tầng bùn mặt.
  • Mômen kháng uốn giới hạn của vật liệu cọc Myield (kNm).
  • Chuyển vị ngang đỉnh cọc y0 (mm).
  • Chiều dài tự do của cọc L0 (m) nằm ngoài mặt đất tự nhiên chịu uốn tự do trong nước.

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

Nguồn dữ liệu của nghiên cứu bao gồm số liệu khảo sát địa chất công trình chi tiết gồm 3 lớp đất chính tại khu vực cảng Thị Vải và chuỗi số liệu thí nghiệm nén ngang cọc thép đường kính D812.1mm ngoài hiện trường với các cấp tải trọng thử nghiệm từ 50 kN đến 400 kN.

Cỡ mẫu nghiên cứu bao gồm 01 cọc thí nghiệm hiện trường quy chuẩn kết hợp với 12 kịch bản phân tích số mô phỏng cọc chịu tải trọng ngang trong các điều kiện biên và chiều dài tự do khác nhau. Phương pháp chọn mẫu áp dụng theo tiêu chuẩn thử nghiệm tải trọng cọc ASTM D4945 và tiêu chuẩn thiết kế móng cọc TCVN 205:1998, đảm bảo tính đại diện cho điều kiện địa tầng đất yếu ven biển phía Nam.

Phương pháp phân tích được thực hiện qua hai bước chính:

  • Mô phỏng không gian ba chiều bằng phần mềm phần tử hữu hạn Plaxis 3D Foundation để tái hiện chính xác trường ứng suất và biến dạng phi tuyến liên tục của khối đất quanh thân cọc.
  • Mô hình hóa kết cấu thanh trong phần mềm SAP2000 tích hợp hệ thống lò xo phi tuyến p-y phân bố dọc theo chiều dài cọc. Lý do lựa chọn giải pháp này là nhằm kết hợp ưu điểm về độ chính xác cơ học của phương pháp 3D với khả năng tính toán nhanh hàng trăm tổ hợp tải trọng gió, sóng, va tàu trên mô hình kết cấu tổng thể. Toàn bộ quá trình nghiên cứu được triển khai liên tục trong thời gian 6 tháng từ tháng 01/2013 đến tháng 06/2013.

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ố liệu và so sánh giữa mô phỏng với thực nghiệm đã chỉ ra các phát hiện cốt lõi sau:

  • Sai lệch của mô hình p-y tiêu chuẩn: Khi áp dụng trực tiếp đường cong p-y lý thuyết của Matlock cho lớp sét mềm tại Thị Vải, độ cứng của nền đất bị đánh giá cao hơn thực tế, làm cho chuyển vị ngang đầu cọc tính toán nhỏ hơn giá trị đo thực tế ngoài hiện trường khoảng 28% đến 35% ở cấp tải trọng 300 kN.
  • Hiệu chỉnh thành công đường cong p-y: Bằng việc đưa vào hệ số suy giảm độ cứng alpha dao động từ 0.65 đến 0.72 cho tầng sét mềm trên cùng (độ sâu từ 0.0m đến -15.0m), kết quả tính toán chuyển vị và mômen uốn trong SAP2000 đạt độ trùng khớp cao với số liệu quan trắc thực tế, giảm sai số chuyển vị đầu cọc xuống còn dưới 5.2%.
  • Tác động sâu sắc của chiều dài tự do cọc L0: Khi chiều dài tự do của cọc tăng từ 0m lên 6m và 10m (tương ứng với mực nước sâu của bến cảng), chuyển vị ngang đầu cọc tăng lần lượt là 145% và 290% dưới cùng một cấp tải trọng ngang 200 kN. Đồng thời, mômen uốn cực đại Mmax tăng từ 40% đến 85% và vị trí xuất hiện mômen cực đại bị đẩy sâu xuống dưới đáy bùn thêm từ 1.5m đến 2.5m.
  • Khả năng thay thế của mô hình lò xo phi tuyến: So sánh giữa mô hình phần tử hữu hạn Plaxis 3D Foundation và mô hình cọc gắn lò xo p-y hiệu chỉnh trong SAP2000 cho thấy mức độ tương đồng đạt trên 92% về phân bố nội lực mômen, khẳng định tính khả thi của mô hình cọc - lò xo trong thiết kế bến cảng thực tế.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự sai khác giữa lý thuyết Matlock và thực tế là do cấu trúc bùn sét mềm ven biển tại khu vực cửa sông Thị Vải có độ bão hòa nước cao và tính nhạy cảm lớn. Dưới tác dụng của tải trọng ngang tăng dần, vùng đất xung quanh cổ cọc bị xáo động và phá hủy dẻo cục bộ nhanh hơn so với các giả định trong điều kiện thí nghiệm tiêu chuẩn.

Kết quả này hoàn toàn nhất quán với các công bố khoa học quốc tế của Reese và các nghiên cứu địa kỹ thuật tại khu vực Cái Mép - Thị Vải, đồng thời khắc phục triệt để nhược điểm của phương pháp ngàm cứng giả định vốn thường làm sai lệch mômen uốn tới 45%.

Dữ liệu nghiên cứu được trực quan hóa thông qua hệ thống biểu đồ quan hệ chuyển vị - tải trọng ngang (P-y0), biểu đồ phân bố mômen uốn M dọc thân cọc theo từng cấp tải trọng (từ 50 kN đến 400 kN) và bảng so sánh các giá trị chuyển vị tương ứng giữa đo thực tế, phần mềm Plaxis 3D và phần mềm SAP2000 trước và sau khi hiệu chỉnh. Sự kết hợp này mang lại cái nhìn định lượng rõ ràng về động thái làm việc của cọc trong nền đất yếu.

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

Dựa trên các kết luận nghiên cứu, bốn giải pháp kỹ thuật cụ thể được đề xuất nhằm nâng cao chất lượng thiết kế và thi công móng cọc công trình cảng:

  • Áp dụng mô hình cọc - lò xo phi tuyến p-y hiệu chỉnh trong các phần mềm tính toán kết cấu phổ thông (như SAP2000, MIDAS) thay thế hoàn toàn cho mô hình ngàm cứng hoặc ngàm trượt giả định. Mục tiêu kiểm soát sai số nội lực cọc dưới 8%, thực hiện ngay trong giai đoạn thiết kế kỹ thuật (thời gian 1-2 tháng đầu dự án) do các kỹ sư kết cấu cảng và kỹ sư địa kỹ thuật đảm trách.
  • Bắt buộc thực hiện thí nghiệm nén tĩnh ngang cọc hiện trường kết hợp lắp đặt thiết bị đo biến dạng sâu (Inclinometer) dọc thân cọc cho 100% các công trình bến cảng tiếp nhận tàu trọng tải từ 20.000 DWT trở lên. Mục tiêu kiểm soát chuyển vị ngang đỉnh cọc không vượt quá ngưỡng cho phép 25 mm, do chủ đầu tư và nhà thầu thí nghiệm chuyên ngành triển khai trong vòng 30 ngày trước khi thi công cọc đại trà.
  • Xây dựng hệ thống cơ sở dữ liệu đường cong p-y đặc thù cho các tầng bùn sét mềm và cát phân bố tại các vùng kinh tế trọng điểm ven biển phía Nam (Bà Rịa - Vũng Tàu, TP. Hồ Chí Minh, Cần Thơ). Mục tiêu tối ưu hóa chiều dài cắm cọc, tiết kiệm từ 12% đến 18% khối lượng vật liệu cọc, do các viện nghiên cứu chuyên ngành phối hợp với các trường đại học kỹ thuật thực hiện trong lộ trình 12 đến 24 tháng.
  • Cập nhật quy định kiểm toán ảnh hưởng của chiều dài tự do cọc L0 trong tiêu chuẩn thiết kế bến cảng thủy lợi - giao thông quốc gia. Đảm bảo hệ số an toàn chịu uốn của cọc đài cao luôn đạt mức tối thiểu từ 1.5 đến 2.0 dưới tác động va neo tàu cực hạn, do Bộ Giao thông Vận tải và Bộ Xây dựng chủ trì rà soát, ban hành.

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

Tài liệu luận văn mang lại giá trị thực tiễn và học thuật cho bốn nhóm đối tượng chính:

  • Kỹ sư thiết kế công trình thủy và cầu cảng: Tiếp cận phương pháp thiết lập ma trận độ cứng lò xo phi tuyến p-y để mô hình hóa chính xác hệ cọc đài cao trên SAP2000, phục vụ phân tích đồng thời hàng trăm tổ hợp tải trọng va tàu, neo tàu và sóng gió phức tạp.
  • Chuyên gia địa kỹ thuật và tư vấn giám sát: Sử dụng quy trình phân tích và công thức hiệu chỉnh hệ số phản lực nền p-y từ dữ liệu nén tĩnh ngang thực tế để lập báo cáo đánh giá an toàn nền móng công trình ven biển.
  • Chủ đầu tư và ban quản lý dự án cảng biển: Có cơ sở khoa học để thẩm tra phương án kết cấu móng móng cọc đài cao, kiểm soát chặt chẽ biến dạng ngang của bến dưới ngưỡng 20 mm và tối ưu hóa tổng mức đầu tư xây dựng công trình.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Địa kỹ thuật - Công trình biển: Sử dụng làm tài liệu tham khảo chuẩn mực về phương pháp kết hợp mô phỏng phần tử hữu hạn 3D (Plaxis 3D) với bài toán phân tích khung kết cấu không gian.

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

Tại sao mô hình ngàm cứng giả định không còn phù hợp trong tính toán cọc bến cảng?

Mô hình ngàm cứng giả định tại một độ sâu cố định không phản ánh đúng tương tác phi tuyến giữa đất và cọc. Trong đất sét mềm, cọc bị biến dạng liên tục theo chiều sâu. Việc dùng ngàm cứng làm sai lệch chuyển vị ngang đầu cọc từ 30% đến 50% và đánh giá sai vị trí mômen uốn nguy hiểm nhất.

Phương pháp đường cong p-y có ưu thế gì so với lý thuyết của Broms (1964)?

Lý thuyết của Broms chỉ dự báo được sức chịu tải giới hạn và chuyển vị tại mặt đất cho nền đất đồng nhất, không cung cấp nội lực dọc thân cọc. Ngược lại, phương pháp p-y mô tả chi tiết biểu đồ mômen, lực cắt và chuyển vị tại mọi cao độ trên cọc qua nhiều lớp đất phức tạp.

Hệ số hiệu chỉnh alpha cho lớp sét mềm tại Thị Vải có ý nghĩa gì?

Hệ số alpha (từ 0.65 đến 0.72) giúp hiệu chỉnh giảm độ cứng ban đầu của đường cong p-y tiêu chuẩn Matlock. Nhờ đó, phản lực đất mô phỏng phản ánh đúng trạng thái bão hòa nước và phá hủy dẻo cục bộ của bùn sét ven biển, đưa sai số tính toán về mức dưới 5.2%.

Chiều dài tự do của cọc L0 ảnh hưởng thế nào đến sự làm việc của móng cảng?

Khi chiều dài tự do tăng do mực nước sâu, độ cứng chống uốn tổng thể của cọc giảm mạnh. Ví dụ, khi L0 tăng từ 0m lên 6m, chuyển vị đầu cọc tăng 145%, làm mômen uốn cực đại tăng trên 40% và dịch chuyển điểm nguy hiểm xuống sâu hơn trong lòng đất.

Có thể áp dụng kết quả nghiên cứu này cho các loại cọc bê tông cốt thép hay cọc ván không?

Hoàn toàn có thể áp dụng. Bản chất phương pháp p-y hiệu chỉnh gắn liền với đặc tính cơ lý của tầng đất nền. Khi thay đổi sang cọc bê tông cốt thép ly tâm dự ứng lực hoặc cọc ván thép, kỹ sư chỉ cần cập nhật lại mô men quán tính I và mô đun đàn hồi Ep của cọc.

Kết luận

  • Luận văn đã giải quyết thành công bài toán tương tác phi tuyến phức tạp giữa cọc và nền đất yếu chịu tải trọng ngang trong công trình bến cảng.
  • Xác định chính xác mức độ sai lệch của phương pháp p-y kinh điển và đề xuất hệ số hiệu chỉnh alpha từ 0.65 đến 0.72 tối ưu cho tầng sét mềm khu vực Thị Vải.
  • Khẳng định tính hiệu quả vượt trội của mô hình cọc - lò xo phi tuyến trong SAP2000 với độ tương đồng đạt trên 92% so với mô phỏng 3D Plaxis.
  • Định lượng rõ nét quy luật gia tăng chuyển vị đỉnh cọc tới 290% và độ dịch chuyển mômen uốn cực đại khi chiều dài tự do cọc trong nước tăng lên 10m.
  • Đóng góp giải pháp tính toán kết cấu móng chuẩn xác, hỗ trợ tiết kiệm 15% chi phí đầu tư xây dựng móng cọc cảng biển mà vẫn duy trì hệ số an toàn cao.

Đóng góp chính của đề tài là cầu nối quan trọng giữa lý thuyết địa kỹ thuật phi tuyến chuyên sâu và thực hành thiết kế kết cấu công trình thủy. Trong kế hoạch 6 đến 12 tháng tiếp theo, các nghiên cứu mở rộng cần tập trung phân tích hiệu ứng nhóm cọc và tác động của tải trọng động chu kỳ dài hạn. Các đơn vị tư vấn và kỹ sư thiết kế công trình cảng được khuyến nghị áp dụng ngay quy trình tính toán p-y hiệu chỉnh này vào các dự án hạ tầng bến cảng nước sâu hiện nay.