Tổng quan nghiên cứu

Khu vực Đồng bằng sông Cửu Long và Thành phố Hồ Chí Minh có hơn 70% diện tích tự nhiên là nền đất yếu trầm tích với chiều dày tầng bùn sét dao động phổ biến từ 15 m đến hơn 30 m. Trong bối cảnh hạ tầng giao thông giai đoạn 2020 đến tầm nhìn 2030 đang được đẩy mạnh quy hoạch, việc xây dựng hàng loạt tuyến cao tốc và cầu vượt sông đặt ra thách thức kỹ thuật rất lớn. Vấn đề cốt lõi phát sinh tại khu vực chuyển tiếp giữa mố cầu và đường dẫn: chiều cao khối đất đắp sau mố thường đạt từ 4 m đến 8 m tạo ra tải trọng thẳng đứng rất lớn, làm nền đất yếu bên dưới dịch chuyển ngang và tác động áp lực bị động lên hệ móng cọc mố cầu. Hàng loạt sự cố nghiêm trọng từng xảy ra trong thực tế như mố cầu Kỳ Hà IV bị chuyển vị ngang từ 5,6 cm đến 9,0 cm, cầu An Nghĩa dịch chuyển đỉnh mố tới 6,0 m trên lớp sét dày 17 m, hay vụ trượt mố cầu Kênh Ngang xô dạt 4,5 m gây thiệt hại hơn 3 tỷ đồng.

Mục tiêu cụ thể của công trình này là phân tích chuyên sâu ảnh hưởng của khối đất đắp đến trạng thái ứng suất và biến dạng của móng cọc khoan nhồi mố cầu thông qua mô phỏng số không gian ba chiều. Phạm vi nghiên cứu tập trung vào tương tác phi tuyến giữa cọc và nền đất yếu đặc thù tại miền Nam trong giai đoạn nghiên cứu từ tháng 01/2016 đến tháng 12/2016. Ý nghĩa khoa học và thực tiễn của đề tài thể hiện qua việc lượng hóa chính xác hiệu ứng nhóm cọc với khoảng cách từ 2,0D đến 4,0D, giúp giảm thiểu sai số dự báo chuyển vị xuống dưới 5%, nâng cao tuổi thọ công trình trên 50 năm và tiết kiệm khoảng 15% chi phí xử lý nền móng.

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 tích hợp các nền tảng lý thuyết cơ học đất và kết cấu công trình tiên tiến nhằm phản ánh sát nhất ứng xử thực tế của hệ tương tác:

  1. Mô hình đàn dẻo Mohr - Coulomb: Đây là mô hình vật liệu chủ đạo dùng để mô phỏng đất nền với 5 thông số cơ bản bao gồm mô đun đàn hồi Young E (dao động từ 1.500 kPa đến 25.000 kPa xác định từ thí nghiệm nén ba trục và chỉ số SPT), hệ số Poisson v (từ 0,25 đến 0,45), lực dính đơn vị c (từ 12 kPa đến 45 kPa), góc ma sát trong phi (từ 5 độ đến 28 độ) và góc giãn nở psi. Mô hình xác định chính xác mặt phá hoại dẻo và thế năng biến dạng dẻo của đất sét trong điều kiện chịu tải không thoát nước.
  2. Lý thuyết dầm đàn hồi trên nền đàn hồi Winkler và phương pháp đường cong p - y: Khảo sát phản lực đất nền phi tuyến p tương ứng với chuyển vị ngang y của thân cọc theo các hàm phi tuyến bậc ba, khắc phục nhược điểm của giả thuyết phản lực tuyến tính cổ điển.
  3. Mô hình đàn hồi tuyến tính (Linear Elastic): Áp dụng cho kết cấu đài mố và cọc khoan nhồi bê tông cốt thép với mô đun đàn hồi vật liệu đạt 30.000 MPa, khối lượng riêng 25 kN/m³ và hệ số Poisson 0,20.
  4. Lý thuyết cơ học tiếp xúc bề mặt: Sử dụng phần tử tiếp xúc phi tuyến nhằm mô phỏng hiện tượng trượt, truyền ứng suất cắt tiếp xúc và khả năng tách hở giữa bề mặt thân cọc bê tông và khối đất xung quanh.

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

Phương pháp nghiên cứu kết hợp giữa phân tích lý thuyết giải tích và mô phỏng số phần tử hữu hạn phi tuyến 3D trên phần mềm ANSYS MAPDL phiên bản 17.0:

  • Cỡ mẫu và phương pháp chọn mẫu mô hình: Nghiên cứu thiết lập hệ thống 7 sơ đồ hình học đại diện gồm mô hình 1 cọc chịu tải ngang (1P_LL), 1 cọc chịu tải đứng (1P), và 5 mô hình cụm cọc đôi chịu tải đứng với khoảng cách trục cọc biến thiên có hệ thống gồm 2,0D; 2,5D; 3,0D; 3,5D và 4,0D (với D là đường kính cọc khoan nhồi). Mỗi sơ đồ được phân tích qua 4 kịch bản vật liệu tiếp xúc khác nhau, tạo nên cỡ mẫu 28 trường hợp mô phỏng số chi tiết. Phương pháp chọn mẫu phân tầng theo cự ly cọc giúp bao quát toàn bộ dải làm việc thực tế của móng mố cầu theo tiêu chuẩn ngành.
  • Nguồn dữ liệu đầu vào: Trích xuất từ số liệu khảo sát địa chất công trình các dự án cầu thực tế tại khu vực Thành phố Hồ Chí Minh và Nam Bộ với tầng sét yếu dày 12 m đến 18 m, xử lý qua các thí nghiệm ba trục không thoát nước UU và cố kết thoát nước CD.
  • Lý do lựa chọn phương pháp phân tích: Sử dụng phần tử khối SOLID187 bậc hai (10 nút) cho cọc và đất kết hợp cặp phần tử tiếp xúc TARGE170 và CONTA174. Thuật toán giải phi tuyến Full Newton - Raphson với kiểm soát dung sai hội tụ khắt khe 0,001 (0,1%) được lựa chọn vì đảm bảo khả năng hội tụ cao khi đất chuyển sang trạng thái dẻo cục bộ.
  • Timeline nghiên cứu: Toàn bộ quá trình lập mô hình, kiểm chuẩn độc lập với tiêu chuẩn thiết kế và trích xuất dữ liệu hoàn thành trong chu kỳ 11 tháng (từ tháng 01/2016 đến tháng 12/2016).

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ố 28 kịch bản trên phần mềm ANSYS đã đem lại các phát hiện định lượng mang tính bước ngoặt:

  1. Chuyển vị ngang đỉnh cọc có sự phân hóa mạnh giữa các hàng cọc: Tại cự ly cọc chuẩn 4,0D, chuyển vị ngang đỉnh cọc trước đạt 14,2 mm trong khi cọc sau đạt 18,6 mm (cao hơn khoảng 31,0%). Khi thu hẹp cự ly cọc về mức 2,0D, chuyển vị cọc sau tăng vọt lên 24,8 mm, chênh lệch hơn 74,6% so với khi cọc làm việc đơn lẻ dưới cùng cấp tải.
  2. Vùng tập trung ứng suất nén cực đại: Ứng suất nén lớn nhất trên thân cọc xuất hiện tập trung tại cao độ từ âm 4,0 m đến âm 6,0 m (ngay tại ranh giới tiếp giáp giữa đáy khối đất đắp và đỉnh lớp bùn sét yếu). Ở khoảng cách 2,0D, ứng suất nén cọc sau đạt đỉnh 18,4 MPa, cao hơn 42,6% so với ứng suất nén cọc trước (chỉ đạt 12,9 MPa).
  3. Hiện tượng chắn tải và hiệu ứng vòm đất: Hàng cọc phía sau đóng vai trò như một bức tường rào cản, gánh chịu tới 62% đến 68% tổng áp lực ngang do khối đất đắp truyền vào móng. Tỷ số chênh lệch ứng suất đất giữa hai lớp đất tăng mạnh từ 1,18 lần ở cự ly 4,0D lên 1,85 lần ở cự ly 2,0D.
  4. Ảnh hưởng của phần tử tiếp xúc bề mặt: Mô hình gắn kết cứng giả định (bỏ qua phần tử tiếp xúc) đánh giá thừa sức chịu tải ngang từ 22% đến 35%, đồng thời đánh giá thấp chuyển vị ngang thực tế khoảng 28% so với mô hình có xét tiếp xúc trượt Mohr - Coulomb.

Thảo luận kết quả

Nguyên nhân cốt lõi của sự chênh lệch nội lực là do tải trọng khối đất đắp phía sau mố ép lớp sét yếu dịch chuyển dẻo theo phương ngang hướng về phía lòng sông. Khi khoảng cách cọc nhỏ hơn 3,0D, các nêm đất dẻo xung quanh từng cọc giao thoa mạnh mẽ, tạo nên hiệu ứng vòm chắn đất khiến cọc sau chịu tải uốn và lực cắt vượt trội. Kết quả này tương đồng với các nghiên cứu thực nghiệm mô hình của Zhao Ming-hua và khẳng định rõ khi cọc đặt trong nhóm chịu tải bị động, sự phân phối tải trọng hoàn toàn không đồng đều như giả định trong các phương pháp tính toán dầm trên nền đàn hồi đơn giản.

Để phục vụ công tác thiết kế trực quan, toàn bộ dữ liệu phân tích có thể được mô tả hiệu quả qua hệ thống biểu đồ phân bố mô men uốn và chuyển vị ngang theo chiều sâu thân cọc từ 0 m đến âm 24 m. Bảng so sánh ứng suất nén lớn nhất giữa các kịch bản 1P, 2P_2D đến 2P_4D thể hiện rõ xu hướng phi tuyến: khoảng cách cọc càng nhỏ thì độ dốc đường cong gia tăng ứng suất càng cao. Biểu đồ đường cong quan hệ áp lực - biến dạng đất tại độ sâu âm 6,0 m giúp các kỹ sư nhận biết chính xác thời điểm đất nền đạt tới giới hạn chảy dẻo cực hạn để đưa ra các biện pháp gia cường kịp thời.

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

Dựa trên các phát hiện khoa học từ mô hình phân tích số, luận văn đề xuất 4 nhóm giải pháp kỹ thuật mang tính thực thi cao:

  • Ban hành quy trình tính toán móng mố 3D phi tuyến bắt buộc: Yêu cầu các đơn vị tư vấn thiết kế chuyển đổi việc tính toán móng mố cầu trên nền đất yếu từ mô hình 2D đàn hồi tuyến tính sang mô hình không gian 3D áp dụng tiêu chuẩn Mohr - Coulomb có phần tử tiếp xúc, nhằm kiểm soát sai số chuyển vị dưới 5% (thời gian áp dụng: ngay từ quý 1/2027 cho các dự án nhóm A và nhóm B).
  • Quy định khoảng cách tối ưu giữa các cọc khoan nhồi trong mố: Các kỹ sư thiết kế cần bố trí khoảng cách giữa cọc trước và cọc sau tối thiểu từ 3,5D đến 4,0D theo phương ngang, giúp phân bổ đồng đều trên 85% áp lực đất bị động, hạn chế tình trạng ứng suất nén vượt ngưỡng cho phép 15 MPa ở hàng cọc sau (thời gian thực hiện: tích hợp trực tiếp vào hồ sơ thiết kế kỹ thuật thi công).
  • Triển khai giải pháp xử lý triệt để nền đất đắp đầu mố: Chủ đầu tư và nhà thầu thi công cần bắt buộc áp dụng các giải pháp giảm tải như sử dụng vật liệu đắp nhẹ EPS, bấc thấm kết hợp hút chân không hoặc cọc xi măng đất để nâng sức chống cắt không thoát nước của đất nền lên trên 35 kPa trước khi đắp đất cao trên 4 m (thời gian xử lý: kiểm soát quá trình cố kết trong 90 đến 120 ngày trước khi gác dầm).
  • Lắp đặt hệ thống quan trắc chuyển vị ngang và ứng suất tự động: Đơn vị tư vấn giám sát cần bố trí các ống đo độ nghiêng thân cọc (inclinometer) và hộp đo áp lực đất tại các độ sâu trọng yếu từ âm 2 m đến âm 15 m, thiết lập ngưỡng cảnh báo sớm khi chuyển vị ngang vượt quá 20 mm (thời gian theo dõi: liên tục 24 tháng trong suốt quá trình thi công và vận hành khai thác).

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

Nội dung luận văn mang lại giá trị học thuật và ứng dụng thực tiễn to lớn cho 4 nhóm đối tượng chuyên môn:

  1. Kỹ sư thiết kế cầu đường và địa kỹ thuật: Nắm vững kỹ thuật lập mô hình phần tử hữu hạn 3D trên ANSYS MAPDL, làm chủ phương pháp khai báo phần tử tiếp xúc TARGE170/CONTA174 và mô hình vật liệu Mohr - Coulomb để tối ưu hóa đường kính cọc từ 1,0 m đến 1,5 m, tiết kiệm 10% đến 15% khối lượng cốt thép.
  2. Nhà thầu thi công công trình giao thông: Hiểu sâu về cơ chế trượt nứt đất đắp đầu cầu cao từ 4 m đến 8 m để lập biện pháp thi công đắp đất phân lớp hợp lý dày không quá 0,5 m mỗi đợt, ngăn ngừa rủi ro gãy nứt cọc khoan nhồi hoặc xô lệch mố cầu.
  3. Cán bộ quản lý dự án và tư vấn giám sát: Sử dụng tài liệu làm cơ sở kỹ thuật để thẩm tra phương án thiết kế móng mố cầu trên đất yếu, giám sát chặt chẽ tiêu chuẩn chuyển vị ngang dưới 25 mm và đánh giá an toàn kết cấu vùng giáp ranh cầu - đường.
  4. Giảng viên, học viên cao học và nghiên cứu sinh ngành Công trình giao thông: Sử dụng như một tài liệu tham khảo chuẩn mực về phương pháp phân tích phi tuyến tương tác kết cấu - đất, kế thừa mã lệnh APDL và phương pháp xử lý số liệu với hơn 30 tài liệu tham khảo chuyên ngành uy tín.

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

  1. Tại sao khối đất đắp sau mố lại gây nguy cơ mất an toàn lớn cho cọc khoan nhồi? Khối đất đắp đầu cầu cao từ 4 m đến 8 m tạo ra tải trọng thẳng đứng rất lớn làm lớp bùn sét yếu bên dưới bị ép trồi và dịch chuyển ngang về phía sông. Lực đẩy bị động này tác động trực tiếp dọc thân cọc gây uốn và cắt nghiêm trọng, từng làm mố cầu An Nghĩa dịch chuyển 6,0 m và gây sụp đổ mố cầu Kênh Ngang.

  2. Mô hình đất Mohr - Coulomb vượt trội hơn mô hình đàn hồi tuyến tính ở điểm nào? Mô hình Mohr - Coulomb mô phỏng chính xác trạng thái biến dạng dẻo của đất qua 5 thông số gồm mô đun E, hệ số Poisson v, lực dính c, góc ma sát phi và góc nở psi. Mô hình đàn hồi tuyến tính không xét đến giới hạn chảy dẻo nên thường đánh giá sai lệch sức chịu tải từ 20% đến 35%.

  3. Bố trí khoảng cách cọc mố cầu bao nhiêu là an toàn và kinh tế nhất? Khoảng cách cọc tối ưu theo phương ngang là từ 3,5D đến 4,0D. Khi khoảng cách cọc nhỏ hơn 3,0D, hiệu ứng vòm đất làm cọc sau phải chịu hơn 65% tổng áp lực ngang, gây tập trung ứng suất nén vượt 18 MPa và tăng nguy cơ nứt vỡ cọc bê tông.

  4. Vai trò của cặp phần tử tiếp xúc TARGE170 và CONTA174 trong mô phỏng là gì? Cặp phần tử này cho phép mô phỏng chính xác hiện tượng trượt tiếp xúc và tách hở giữa bề mặt cọc bê tông và khối đất nền. Nếu giả định dính kết cứng tuyệt đối, chuyển vị ngang tính toán sẽ bị đánh giá thấp hơn thực tế khoảng 28%, dẫn đến thiết kế thiếu an toàn.

  5. Giải pháp thực tế nào hiệu quả nhất để hạn chế chuyển vị ngang mố cầu? Giải pháp tối ưu là kết hợp gia cố nền đất yếu bằng cọc xi măng đất nâng sức chống cắt lên trên 35 kPa, duy trì khoảng cách cọc trên 3,5D, thi công đắp đất gia tải từng lớp mỏng dưới 0,5 m và quan trắc chuyển vị liên tục không để vượt ngưỡng 20 mm.

Kết luận

  • Công trình đã thiết lập thành công mô hình số 3D phi tuyến trên ANSYS MAPDL để đánh giá toàn diện ảnh hưởng của khối đất đắp đến móng cọc mố cầu trên nền đất yếu.
  • Phát hiện quy luật phân bố tải trọng bất đối xứng: cọc sau chịu ứng suất nén và chuyển vị ngang lớn hơn cọc trước tới 42,6% do hiện tượng chắn và hiệu ứng vòm đất khi cự ly cọc dưới 3,0D.
  • Khẳng định cự ly bố trí cọc từ 3,5D đến 4,0D là giải pháp hình học tối ưu giúp cân bằng nội lực giữa các hàng cọc và giảm thiểu biến dạng công trình.
  • Chứng minh tính tất yếu của việc áp dụng mô hình Mohr - Coulomb kết hợp phần tử tiếp xúc bề mặt nhằm loại bỏ sai số tính toán lên đến 35% của các phương pháp truyền thống.
  • Đóng góp luận cứ khoa học quan trọng cho việc hoàn thiện tiêu chuẩn thiết kế và xử lý sự cố móng mố cầu trên tầng đất yếu dày trên 15 m tại khu vực phía Nam.

Hướng nghiên cứu tiếp theo sẽ mở rộng đánh giá cho các dạng móng cọc xiên và sơ đồ cọc bố trí bất đối xứng trong khung thời gian 12 tháng tới. Các chuyên gia, kỹ sư và nhà quản lý hãy tham khảo ngay luận văn này để ứng dụng các giải pháp mô phỏng số tiên tiến vào dự án thực tế, đảm bảo an toàn tuyệt đối và nâng cao hiệu quả kinh tế cho công trình cầu đường!