Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình ngầm và địa kỹ thuật hiện đại, móng cọc bê tông cốt thép đóng vai trò là giải pháp kết cấu chịu lực nền móng cốt lõi nhằm chuyển tải trọng siêu trường từ thượng tầng xuống các tầng địa chất ổn định. Tuy nhiên, việc xác định chính xác sức chịu tải của cọc đơn luôn là thách thức lớn đối với các kỹ sư do độ bất định cao của đất nền và sự chênh lệch đáng kể giữa lý thuyết tính toán với thực tế thi công. Tại nhiều dự án hạ tầng quy mô lớn, sai số giữa các công thức giải tích cổ điển và kết quả thí nghiệm nén tĩnh hiện trường thường dao động trong khoảng 20% đến 45%, dẫn đến lãng phí chi phí vật liệu hoặc tiềm ẩn rủi ro lún lệch nghiêm trọng.

Nghiên cứu tập trung giải quyết bài toán phân tích sức chịu tải của cọc và làm rõ mức độ ảnh hưởng của ma sát thành bên đến khả năng chịu lực tổng thể của cọc bê tông cốt thép. Mục tiêu trọng tâm là thiết lập mối tương quan chuẩn xác giữa phương pháp tính toán giải tích, kết quả thí nghiệm nén tĩnh hiện trường và phương pháp mô phỏng số bằng phần tử hữu hạn. Phạm vi nghiên cứu được thực nghiệm trực tiếp tại Dự án Kho chứa lạnh LPG Thị Vải thuộc Khu công nghiệp Cái Mép, huyện Tân Thành, tỉnh Bà Rịa – Vũng Tàu. Công trình có quy mô diện tích 5,3 ha với công suất thiết kế 60.000 tấn LPG lạnh, đặt trên nền địa chất ven biển phân tầng phức tạp. Đối tượng khảo sát là cọc bê tông cốt thép đúc sẵn kích thước tiết diện 350x350 mm với chiều dài hạ cọc 34 m ngàm qua 5 lớp đất đặc trưng. Ý nghĩa của đề tài thể hiện qua việc nâng cao độ tin cậy trong dự báo sức chịu tải ở giai đoạn thiết kế cơ sở, giúp tối ưu hóa từ 10% đến 15% chiều dài cọc và giảm thiểu chi phí thử tải tĩnh tốn kém trong giai đoạn chuẩn bị dự án.

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ơ học đất và nền móng chuyên sâu với các hệ thống lý thuyết kinh điển:

  • Lý thuyết cân bằng giới hạn phá hoại dẻo của Terzaghi kết hợp mô hình phá hoại trượt đất nền của Meyerhof và Vesic. Các mô hình này thiết lập công thức tổng quát xác định sức kháng mũi dựa trên chỉ số độ cứng và góc ma sát trong của đất nền.
  • Lý thuyết tương tác đất - cọc theo phương pháp alpha của Tomlinson đối với đất dính và phương pháp beta của Burland đối với đất rời, cho phép tích phân lực dính đơn vị và ứng suất pháp hiệu quả dọc theo chiều dài cọc.
  • Phương pháp phân tích CPTu hiện đại của Eslami và Fellenius, ứng dụng chỉ số sức kháng mũi hiệu quả để phân loại đất và xác định ma sát thành bên.
  • Mô hình đàn dẻo Mohr-Coulomb trong cơ học môi trường liên tục, mô tả quan hệ ứng suất - biến dạng phi tuyến của đất xung quanh thân cọc và dưới mũi cọc thông qua 5 thông số đầu vào chính: mô đun đàn hồi biến dạng $E$, hệ số Poisson $\nu$, lực dính $c$, góc ma sát trong $\varphi$, và góc giãn nở $\psi$.

Các khái niệm trọng tâm bao gồm sức chịu tải cực hạn $Q_u$, sức kháng ma sát bên $Q_s$, sức kháng mũi $Q_p$, và hệ số an toàn tính toán $FS$ dao động từ 1,5 đến 3,0 theo tiêu chuẩn xây dựng TCXDVN 205:1998 và TCVN 10304:2014. Theo quy luật cơ học tiếp xúc, sức kháng bên đạt giá trị cực hạn rất nhanh ở độ lún đầu cọc chỉ từ 3 đến 5 mm, trong khi sức kháng mũi đòi hỏi chuyển vị lớn hơn nhiều để huy động toàn bộ cường độ kháng cắt.

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

Nguồn dữ liệu của nghiên cứu được thu thập từ hồ sơ khảo sát địa kỹ thuật chi tiết của hố khoan A-BH01 sâu 36,5 m kết hợp dữ liệu thí nghiệm xuyên tiêu chuẩn SPT được lấy mẫu định kỳ ở khoảng cách 1,5 đến 2,0 m dọc theo chiều sâu. Cỡ mẫu nghiên cứu bao gồm toàn bộ chuỗi số liệu địa chất của 5 phân lớp đất tại vị trí tim cọc thử nghiệm, kết hợp số liệu quan trắc chuyển vị từ 4 đồng hồ đo lún có độ chính xác 0,01 mm trong thí nghiệm nén tĩnh dọc trục hiện trường theo quy trình TCVN 269:2002.

Phương pháp chọn mẫu địa kỹ thuật được thực hiện thông qua lấy mẫu đất nguyên dạng bằng ống mẫu thành mỏng cho các tầng sét dẻo mềm và lấy mẫu xáo nhét bằng mũi thìa bổ đôi cho các tầng cát san lấp và cát pha sét. Nghiên cứu kết hợp hai nhóm phương pháp phân tích chính:

  • Phương pháp giải tích: Tính toán sức chịu tải cọc theo chỉ tiêu cường độ đất nền, chỉ tiêu cơ lý và kết quả thí nghiệm xuyên tiêu chuẩn SPT theo công thức Meyerhof và công thức của Viện Kiến trúc Nhật Bản.
  • Phương pháp mô phỏng số phần tử hữu hạn: Sử dụng phần mềm Plaxis 2D phiên bản V8.6 với bài toán đối xứng trục mô phỏng cọc bê tông mác B35 có mô đun đàn hồi $E = 3,45 \times 10^7\text{ kN/m}^2$. Lý do lựa chọn mô hình Plaxis 2D là khả năng tái hiện chân thực trường biến dạng phi tuyến, phân tách độc lập tỷ phần ma sát thành cọc và phản lực mũi cọc trong từng cấp gia tải từ 25% đến 200% tải trọng thiết kế, điều mà các thí nghiệm nén tĩnh thông thường không gắn cảm biến đo sâu khó có thể bóc tách chi tiết. Toàn bộ quá trình nghiên cứu và xử lý số liệu được triển khai liên tục từ tháng 01/2015 đến tháng 08/2015.

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

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

Quá trình đối chiếu giữa tính toán lý thuyết, mô phỏng số và dữ liệu nén tĩnh hiện trường tại công trình Kho chứa LPG Thị Vải đã đưa ra các phát hiện cơ bản sau:

  • Tỷ trọng đóng góp vượt trội của ma sát thành bên: Trong tổng sức chịu tải cực hạn của cọc dài 34 m, thành phần sức kháng ma sát bên chiếm tỷ lệ áp đảo từ 68% đến 73%, trong khi sức kháng mũi cọc chỉ đóng góp khoảng 27% đến 32%. Nguyên nhân do cọc đi qua tầng cát san lấp dày 7,3 m và tầng sét dẻo mềm dày 14,3 m tạo ra diện tích tiếp xúc xung quanh rất lớn.
  • Cơ chế huy động sức kháng theo chuyển vị: Sức ma sát thành bên đạt đến 90% giá trị cực hạn ngay khi độ lún đầu cọc ở mức 3,5 mm đến 4,8 mm. Ngược lại, sức kháng mũi chỉ mới huy động được khoảng 30% đến 40% khả năng chịu lực tại cùng cấp chuyển vị này và chỉ đạt trạng thái cực hạn khi độ lún vượt quá 8% đến 10% kích thước cạnh cọc, tương đương 28 mm đến 35 mm.
  • Độ phân tán giữa các công thức giải tích: Phương pháp SPT theo Meyerhof cho kết quả sức chịu tải giới hạn cao hơn khoảng 18% so với phương pháp của Viện Kiến trúc Nhật Bản. Khi so sánh với thí nghiệm nén tĩnh thực tế, công thức Viện Kiến trúc Nhật Bản có độ tương thích cao hơn với độ lệch chỉ khoảng 9,5%, trong khi công thức Meyerhof có xu hướng đánh giá cao sức kháng mũi ở các tầng cát chặt.
  • Vùng tập trung ứng suất và biến dạng dẻo: Mô phỏng Plaxis 2D cho thấy các điểm phá hoại dẻo theo tiêu chuẩn Mohr-Coulomb xuất hiện đầu tiên dọc theo bề mặt tiếp xúc cọc - đất ở độ sâu từ 2 m đến 8 m, sau đó lan dần xuống chân cọc khi gia tải lên mức 150% đến 200% tải trọng thiết kế.

Thảo luận kết quả

Các kết quả nghiên cứu phản ánh sự làm việc đồng thời nhưng không đồng pha giữa ma sát thân cọc và sức chống mũi cọc. Sự chênh lệch giữa các phương pháp tính toán giải tích xuất phát từ việc các công thức cổ điển xem đất nền là vật liệu đồng nhất và đẳng hướng, bỏ qua hiệu ứng suy giảm độ cứng khi biến dạng tăng dần.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị quan hệ tải trọng - độ lún dạng đường cong chữ S mượt mà, phản ánh quá trình chuyển tiếp từ giai đoạn đàn hồi tuyến tính sang giai đoạn biến dạng dẻo cục bộ. Đồng thời, biểu đồ phân bố ứng suất cắt dọc thân cọc theo độ sâu minh chứng rõ nét ảnh hưởng của từng tầng địa chất: ứng suất ma sát tăng vọt tại lớp cát pha sét cứng ở độ sâu từ 27,4 m đến 31,0 m với chỉ số SPT đạt 27 đến 34 búa, đóng góp hơn 35% tổng sức ma sát bên của cọc.

So với thí nghiệm nén tĩnh thực tế, kết quả mô phỏng Plaxis 2D có sai số chuyển vị đầu cọc dưới 8,5% ở các cấp tải trọng từ 50% đến 150%. Điều này khẳng định việc sử dụng phần mềm phần tử hữu hạn với các thông số địa kỹ thuật được hiệu chỉnh qua thí nghiệm SPT và CPTu hoàn toàn có khả năng thay thế hoặc hỗ trợ chuẩn xác cho công tác nén tĩnh thăm dò trước thi công.

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

Dựa trên các kết luận rút ra từ mô hình phân tích và kiểm chứng thực nghiệm, 4 nhóm giải pháp kỹ thuật cụ thể được đề xuất nhằm nâng cao hiệu quả thiết kế và thi công móng cọc:

  1. Ứng dụng mô phỏng số Plaxis 2D trong giai đoạn thiết kế cơ sở: Các đơn vị tư vấn thiết kế cần bắt buộc đưa mô phỏng phần tử hữu hạn vào giai đoạn lập dự án cho các công trình có tải trọng lớn trên nền đất yếu. Việc thiết lập mô hình đối xứng trục với hệ số tương tác tiếp xúc đất - cọc $R_{inter} = 0,65 - 0,75$ giúp dự báo chính xác đường cong quan hệ tải trọng - chuyển vị, hoàn thành trong khung thời gian 7 đến 10 ngày làm việc và giảm thiểu từ 10% đến 15% chi phí dự phòng phát sinh.
  2. Tối ưu hóa cao độ mũi cọc dựa trên phân tích ma sát thành bên: Đơn vị thiết kế cần tận dụng tối đa chiều dài cọc ngàm vào tầng cát chịu lực tốt từ 3 đến 5 lần đường kính cọc (tương đương 1,0 m đến 1,8 m đối với cọc 350x350 mm). Biện pháp này giúp huy động đồng thời trên 85% sức ma sát bên và 60% sức kháng mũi, kiểm soát độ lún công trình dưới ngưỡng cho phép 15 mm.
  3. Chuẩn hóa quy trình quan trắc nén tĩnh kết hợp cảm biến đo sâu: Chủ đầu tư và nhà thầu thí nghiệm cần lắp đặt bổ sung chuỗi cảm biến đo biến dạng cốt thép (Strain Gauge) ở các độ sâu 7 m, 21 m và 34 m trong các thí nghiệm nén tĩnh thăm dò. Mục tiêu là phân tách chính xác 100% tỷ phần ma sát và sức kháng mũi thực tế tại hiện trường, thực hiện trong vòng 48 giờ thử tải theo tiêu chuẩn TCVN 269:2002.
  4. Xây dựng cơ sở dữ liệu số hóa địa kỹ thuật khu vực ven biển: Các cơ quan quản lý xây dựng và viện nghiên cứu chuyên ngành cần chủ trì thiết lập hệ thống dữ liệu địa tầng số hóa cho khu vực Cái Mép - Thị Vải trong lộ trình 12 đến 24 tháng. Việc lưu trữ thống nhất các chỉ tiêu cơ lý, biểu đồ SPT và kết quả nén tĩnh của hơn 30 dự án lân cận sẽ rút ngắn 20% thời gian khảo sát và nâng cao độ an toàn cho các công trình hạ tầng cảng biển tương lai.

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

Công trình nghiên cứu mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế địa kỹ thuật và kết cấu nền móng: Nắm bắt phương pháp tính toán sức chịu tải theo tiêu chuẩn mới, hiểu rõ cơ chế huy động ma sát bên để lựa chọn chiều dài cọc tối ưu, giúp giảm từ 10% đến 20% khối lượng vật liệu cọc trong các dự án công nghiệp.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Công trình ngầm: Tiếp cận quy trình mô hình hóa bài toán tương tác đất - kết cấu trên phần mềm Plaxis 2D, phương pháp hiệu chỉnh thông số Mohr-Coulomb từ thí nghiệm hiện trường phục vụ các đề tài nghiên cứu chuyên sâu.
  • Đơn vị tư vấn thẩm tra và giám sát thi công: Sử dụng kết quả nghiên cứu làm cơ sở đối chuẩn, đánh giá độ tin cậy của các báo cáo khảo sát địa kỹ thuật và kiểm soát chặt chẽ các chỉ tiêu kỹ thuật trong thí nghiệm nén tĩnh cọc đạt tải trọng thử 200% đến 300% tải trọng thiết kế.
  • Ban quản lý dự án và chủ đầu tư hạ tầng cảng biển: Có thêm công cụ khoa học để thẩm định phương án móng sâu, kiểm soát chi phí thi công nền móng chiếm từ 25% đến 35% tổng mức đầu tư xây dựng công trình.

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

  1. Tại sao sức kháng ma sát bên lại huy động trước sức kháng mũi cọc trong quá trình gia tải? Theo nguyên lý truyền lực trong môi trường đất dẻo, lực ma sát thành bên được kích hoạt ngay khi có sự dịch chuyển tương đối rất nhỏ giữa bề mặt bê tông và các hạt đất xung quanh, đạt cực hạn ở độ lún 3 đến 5 mm. Trong khi đó, đất dưới mũi cọc chịu trạng thái ứng suất khối ba trục phức tạp, đòi hỏi độ lún lớn từ 25 mm đến 35 mm mới hình thành trọn vẹn nêm phá hoại dẻo.

  2. Mô hình Mohr-Coulomb trong Plaxis 2D có mô phỏng chính xác ứng xử của đất nền hay không? Mô hình Mohr-Coulomb phản ánh rất tốt trạng thái cân bằng giới hạn và phân bố ứng suất với độ lệch kết quả nén tĩnh dưới 10% nếu các thông số đầu vào như mô đun biến dạng $E$ và góc ma sát trong $\varphi$ được hiệu chỉnh chuẩn xác từ thí nghiệm SPT và CPTu. Tuy nhiên, đối với các bài toán dỡ tải hoặc biến dạng nhỏ, mô hình Hardening Soil sẽ mang lại độ chính xác cao hơn.

  3. Khi nào thí nghiệm nén tĩnh cọc hiện trường được phép dừng lại? Theo tiêu chuẩn TCVN 269:2002, thí nghiệm dừng khi tổng chuyển vị đầu cọc vượt quá 10% đường kính cọc (khoảng 35 mm đối với cọc 350x350 mm), độ lún tại một cấp tải tăng gấp 5 lần cấp trước đó, hoặc tốc độ lún không đạt trạng thái ổn định quy ước dưới 0,1 mm/h sau 24 giờ giữ tải liên tục.

  4. Phương pháp tính toán SPT theo Meyerhof khác gì so với công thức Viện Kiến trúc Nhật Bản? Công thức Meyerhof sử dụng trực tiếp chỉ số SPT trung bình với hệ số giới hạn sức kháng mũi $k_p = 400$, trong khi tiêu chuẩn Nhật Bản phân tách chi tiết thành phần đất rời và đất dính thông qua hệ số hiệu chỉnh độ mảnh $L/d$ và sức kháng cắt không thoát nước $c_u$, giúp giảm thiểu độ phân tán số liệu từ 10% đến 15% trong nền đất sét yếu.

  5. Cần lưu ý gì khi lựa chọn hệ số tương tác tiếp xúc đất - cọc trong phần mềm Plaxis? Hệ số tương tác tiếp xúc cọc - đất $R_{inter}$ quyết định trực tiếp đến giá trị lực ma sát thành bên. Đối với cọc bê tông cốt thép đúc sẵn đóng hoặc ép trong đất sét và cát chặt vừa, giá trị khuyến nghị nên dao động trong khoảng 0,65 đến 0,75 để mô phỏng chính xác sự trượt tương đối tại mặt tiếp xúc mà không làm méo mó lưới phần tử hữu hạn.

Kết luận

  • Luận văn đã phân tích toàn diện và làm rõ cơ chế làm việc của móng cọc bê tông cốt thép 350x350 mm dài 34 m ngàm qua nền đất yếu ven biển tại khu vực Cái Mép - Thị Vải.
  • Chứng minh thành phần ma sát thành bên đóng vai trò chủ đạo, chiếm từ 68% đến 73% tổng sức chịu tải cực hạn và được huy động tối đa ở cấp độ lún rất nhỏ từ 3 đến 5 mm.
  • Khẳng định mô hình phần tử hữu hạn Plaxis 2D là công cụ mô phỏng đáng tin cậy với độ sai lệch so với thí nghiệm nén tĩnh thực tế dưới 8,5%, hỗ trợ đắc lực cho công tác thiết kế sơ bộ.
  • Cung cấp bảng so sánh định lượng chi tiết giữa các phương pháp giải tích hiện hành, chỉ ra công thức Viện Kiến trúc Nhật Bản có độ tương thích cao nhất với thực tế địa chất miền Nam Việt Nam.
  • Đề xuất lộ trình 12 đến 24 tháng trong việc số hóa dữ liệu địa kỹ thuật và chuẩn hóa quy trình tính toán số cho các công trình ngầm và hạ tầng ven biển.

Các kỹ sư kết cấu, chuyên gia địa kỹ thuật và các nhà quản lý dự án hãy ứng dụng ngay phương pháp mô phỏng số kết hợp hiệu chỉnh số liệu nén tĩnh để nâng cao chất lượng thiết kế, đảm bảo an toàn tuyệt đối và tối ưu hóa ngân sách xây dựng nền móng công trình.