Tổng quan nghiên cứu

Trong bối cảnh đô thị hóa nhanh chóng tại các đô thị đặc biệt như Thành phố Hồ Chí Minh, nhu cầu xây dựng các tòa nhà cao tầng tích hợp từ 3 đến 5 tầng hầm phục vụ hạ tầng kỹ thuật và bãi đỗ xe ngày càng gia tăng. Việc thi công hố đào sâu trong khu vực nội đô mật độ cao đối mặt với thách thức lớn từ điều kiện địa tầng phức tạp, lớp đất sét yếu phân bố dày và mực nước ngầm cao. Vấn đề chuyển vị của tường chắn hố đào là yếu tố then chốt ảnh hưởng trực tiếp đến an toàn kết cấu và các công trình lân cận. Trong thực tế tính toán thiết kế truyền thống, các kỹ sư thường coi thông số cơ lý của đất là hằng số cố định, dẫn đến sai số dự báo chuyển vị có thể vượt quá 35% so với thực tế.

Nghiên cứu của tác giả Nguyễn Hoàng Khánh dưới sự hướng dẫn của Tiến sĩ Lê Trọng Nghĩa tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh tập trung vào đề tài: "Xác định các thông số mô hình đất bằng cách phân tích ngược chuyển vị của tường chắn tầng hầm thi công bằng phương pháp Bottom - Up". Mục tiêu trọng tâm của luận văn là thiết lập quy trình phân tích ngược bằng mô hình số nhằm xác định chính xác các thông số biến dạng và cường độ thực tế của đất nền xung quanh hố đào.

Phạm vi nghiên cứu được triển khai trên công trình thực tế là dự án cao ốc President Place tại địa chỉ 93 Nguyễn Du, Quận 1, Thành phố Hồ Chí Minh. Công trình có kích thước mặt bằng 50m x 30m, gồm 12 tầng nổi và 3 tầng hầm, độ sâu hố đào lớn nhất đạt 13,5m tại vị trí hố pit thang máy và tường vây bê tông cốt thép dày 600mm được cắm sâu 21,0m. Ý nghĩa thực tiễn của đề tài thể hiện qua việc giảm thiểu sai số dự báo chuyển vị ngang xuống dưới 8%, tối ưu hóa sơ đồ giằng chống và ngăn ngừa triệt để nguy cơ trồi đáy móng hay nứt lún các công trình di tích lịch sử liền kề.

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 xây dựng trên nền tảng lý thuyết tương tác đất - kết cấu (Soil-Structure Interaction) và cơ học đất phi tuyến hiện đại. Đề tài vận dụng hai lý thuyết đàn - dẻo kinh điển trong mô phỏng địa kỹ thuật:

Thứ nhất là mô hình Mohr-Coulomb (MC), mô tả trạng thái ứng xử đàn hồi tuyến tính - dẻo lý tưởng với mặt trượt cố định. Mô hình này được định nghĩa bởi 5 thông số đặc trưng cơ bản: mô đun đàn hồi Young (E), hệ số Poisson ($\nu$), lực dính hiệu quả ($c'$), góc ma sát trong ($\varphi'$) và góc giãn nở ($\psi$). Nghiên cứu xem xét mô hình này dưới cả hai điều kiện ứng xử thoát nước (Drained) và không thoát nước (Undrained).

Thứ hai là mô hình Hardening Soil (HS), mô hình cao cấp cho phép mô phỏng đường cong ứng suất - biến dạng dạng hyperbol theo nguyên lý dẻo hóa tái bền. Khung lý thuyết này phân tách độ cứng của đất thành ba đại lượng độc lập gồm độ cứng cát tuyến trong thí nghiệm ba trục ($E_{50}^{ref}$), độ cứng tiếp tuyến nén cố kết ($E_{oed}^{ref}$), và độ cứng dỡ tải/nén lại ($E_{ur}^{ref}$), kết hợp với số mũ biến thiên ứng suất $m$.

Các khái niệm trọng tâm bao gồm: áp lực đất chủ động và bị động theo Rankine - Coulomb, chuyển vị ngang tích lũy của tường vây (Horizontal Displacement), và kỹ thuật phân tích ngược (Back-analysis) nhằm hiệu chỉnh ma trận thông số đầu vào từ dữ liệu dịch chuyển đo đạc thực tế.

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

Nghiên cứu thu thập dữ liệu hiện trường từ dự án President Place với cỡ mẫu quan trắc gồm 6 ống đo nghiêng Inclinometer (ký hiệu từ Inclinometer 1 đến Inclinometer 6) được đặt dọc theo chiều sâu 21,0m của thân tường vây D600. Phương pháp chọn mẫu là đo đạc toàn diện tại các vị trí xung yếu ở các cạnh hố đào dài 50m và rộng 30m, ghi nhận liên tục qua 6 chu kỳ thi công tương ứng với từng giai đoạn hạ độ sâu đất.

Phương pháp phân tích số được thực hiện thông qua phần mềm phần tử hữu hạn chuyên ngành Plaxis 3D Foundation phiên bản 1.6. Lý do lựa chọn mô hình 3D thay vì 2D là vì hố đào có dạng hình học chữ nhật hữu hạn (50m x 30m), nơi hiệu ứng vòm không gian và sự phân bố tải trọng tại các góc tường làm việc khác biệt rõ rệt so với bài toán biến dạng phẳng. Mô hình 3D cho phép mô phỏng chính xác từng phần tử kết cấu gồm: phần tử tường liên tục (Wall), hệ dầm giằng (Waling Beam), hệ thanh chống thép hình H350 và H400, cùng hệ cột chống Kingpost 350x350x12x19 cắm vào nền đất.

Timeline nghiên cứu được chia làm 4 giai đoạn logic: khảo sát địa chất và thu thập dữ liệu thí nghiệm hiện trường/trong phòng; lập mô hình hình học và chia lưới phần tử hữu hạn 3D; mô phỏng trình tự thi công Bottom-Up qua các giai đoạn đào đất từ cao độ tự nhiên xuống -3,6m, -6,7m, -9,8m và -13,5m; cuối cùng là phân tích ngược lặp để hiệu chỉnh thông số cơ lý của đất nhằm đạt độ tương thích cao nhất với số liệu quan trắc.

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 ngược trên Plaxis 3D Foundation kết hợp dữ liệu quan trắc thực tế đã mang lại 4 phát hiện kỹ thuật quan trọng:

Thứ nhất, chuyển vị ngang cực đại của tường vây D600 trong thực tế đạt giá trị từ 18,2 mm đến 23,8 mm, tập trung chủ yếu ở độ sâu từ -7,0m đến -11,0m so với mặt đất hiện hữu. Kết quả mô phỏng ban đầu với thông số khảo sát địa chất thuần túy cho ra chuyển vị lên tới 34,5 mm, chênh lệch hơn 45% so với thực tế.

Thứ hai, mô đun biến dạng tham chiếu ($E_{ref}$) của các lớp đất sét dẻo và cát pha thực tế cao hơn đáng kể so với kết quả thí nghiệm nén một trục không nở hông trong phòng. Khi tiến hành phân tích ngược và hiệu chỉnh $E_{ref}$ tăng từ 22.000 kN/m2 lên mức 34.000 đến 38.000 kN/m2 (tăng tương ứng 54,5% đến 72,7%), đường cong chuyển vị mô phỏng đã tiệm cận với sai số dưới 5% so với dữ liệu đo đạc của Inclinometer.

Thứ ba, sự thay đổi của lực dính hiệu quả $c'$ có tác động rõ rệt đến chuyển vị phần đỉnh tường. Cụ thể, khi $c'$ thay đổi từ 27 kN/m2 lên 33 kN/m2 (mức tăng 22,2%), độ võng ngang của tường ở giai đoạn đào đất lớp 1 giảm khoảng 14,3%, trong khi góc ma sát trong $\varphi'$ đóng vai trò quyết định đến độ ổn định trượt sâu ở chân tường cắm sâu 21m.

Thứ tư, hệ giằng chống shoring 3 lớp phát huy hiệu quả chịu lực vượt trội: lớp 1 lắp đặt tại cao độ -3,1m (thép H350), lớp 2 tại cao độ -6,2m (thép H350) và lớp 3 tại cao độ -9,3m (thép H400) đã kiểm soát và giữ cho độ nghiêng thân tường luôn nằm trong giới hạn an toàn cho phép theo tiêu chuẩn xây dựng đô thị.

Thảo luận kết quả

Nguyên nhân của sự khác biệt giữa kết quả tính toán ban đầu và dữ liệu quan trắc xuất phát từ ứng xử dỡ tải (Unloading) của khối đất bên trong hố móng khi đào sâu đến -13,5m và hiện tượng tăng cứng cục bộ do quá trình thoát nước một phần dưới đáy móng. Mô hình Mohr-Coulomb tiêu chuẩn không phản ánh hết tính phụ thuộc ứng suất của mô đun biến dạng, trong khi mô hình Hardening Soil với việc điều chỉnh $E_{ur}^{ref} \approx 3E_{50}^{ref}$ đã mô phỏng chuẩn xác phản ứng đàn hồi của đáy hố đào.

So sánh với các biểu đồ bán thực nghiệm quốc tế như đường cong thiết kế chuyển vị của Clough và O’Rourke (1990) cũng như dạng thức lún lõm lòng chảo của Hsieh và Ou (1998), hình dạng biến dạng tường của dự án President Place hoàn toàn phù hợp với dạng lún Spandrel kết hợp lòng chảo trong đất dính. Tuy nhiên, phương pháp số 3D kết hợp phân tích ngược cho phép xác định chính xác ứng suất cục bộ tại các vị trí giao thoa giữa hệ Kingpost và dầm giằng, điều mà các công thức kinh nghiệm không thể chỉ ra.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua: Biểu đồ đường cong chuyển vị ngang theo chiều sâu (Depth vs Lateral Displacement Profiles) so sánh giữa 6 chu kỳ Inclinometer và kết quả chạy Plaxis 3D; Bảng đối chiếu ma trận thông số địa kỹ thuật trước và sau phân tích ngược; cùng các mặt cắt đồ họa 3D thể hiện trường ứng suất hữu hiệu và biến dạng dẻo tại cao độ đáy hố đào -13,5m.

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

Dựa trên kết quả phân tích ngược của công trình, tác giả đề xuất 4 nhóm giải pháp kỹ thuật có tính ứng dụng cao:

Thứ nhất, chuẩn hóa quy trình phân tích ngược 6 bước trong công tác thiết kế hố đào sâu: Các đơn vị tư vấn thiết kế cần bắt buộc thiết lập chu trình vòng lặp cập nhật thông số địa kỹ thuật ngay sau khi hoàn thành giai đoạn đào đất lớp đầu tiên ở độ sâu -3,6m, thay vì giữ nguyên mô hình tĩnh cho toàn bộ quá trình thi công.

Thứ hai, điều chỉnh hệ số mô đun biến dạng $E'$ trong tính toán thiết kế ban đầu: Đối với khu vực địa chất Quận 1, Thành phố Hồ Chí Minh, các kỹ sư địa kỹ thuật nên áp dụng hệ số hiệu chỉnh tăng từ 1,3 đến 1,6 lần đối với mô đun đàn hồi xác định từ thí nghiệm nén cố kết trong phòng cho các lớp đất sét pha cát có độ sâu từ 6m đến 20m, nhằm tránh việc thiết kế dư thừa kết cấu chống giữ gây lãng phí chi phí xây dựng.

Thứ ba, thiết lập hệ thống quan trắc tự động hóa với tần suất đo 24 giờ/lần: Ban chỉ huy công trường cần duy trì việc đọc dữ liệu Inclinometer với chu kỳ tối thiểu 1 ngày/lần trong các giai đoạn đào đất nhạy cảm từ -6,7m xuống -13,5m, đặt ngưỡng cảnh báo chuyển vị tại mức 25 mm để kích hoạt giải pháp gia tải kích đầu thanh chống shoring kịp thời.

Thứ tư, tối ưu hóa tiến độ lắp đặt hệ giằng chống thép hình: Quy định thời gian tối đa từ khi đào đất đến cao độ thiết kế (-3,6m, -6,7m, -9,8m) đến khi lắp đặt và khóa nêm hệ dầm giằng H350/H400 không được vượt quá 48 giờ, hạn chế tối đa biến dạng từ biến của đất sét yếu xung quanh tường vây D600.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng chuyên môn:

Nhóm 1: Kỹ sư thiết kế địa kỹ thuật và công trình ngầm. Lợi ích: Nắm vững quy trình thiết lập mô hình Plaxis 3D Foundation, lựa chọn đúng thông số mô hình Mohr-Coulomb và Hardening Soil, từ đó nâng cao chất lượng thiết kế hệ tường vây D600 - D800 cho hố đào từ 3 đến 5 tầng hầm.

Nhóm 2: Kỹ sư giám sát và chỉ huy trưởng thi công tại hiện trường. Lợi ích: Hiểu rõ cơ chế biến dạng của hố móng thi công theo phương pháp Bottom-Up, biết cách đọc hiểu và phân tích biểu đồ chuyển vị Inclinometer để đưa ra quyết định xử lý kỹ thuật an toàn khi thi công hệ shoring và Kingpost.

Nhóm 3: Giảng viên, học viên cao học và sinh viên chuyên ngành Địa kỹ thuật - Công trình ngầm. Lợi ích: Sử dụng đề tài như một tài liệu tham khảo học thuật chuẩn mực về phương pháp phân tích ngược (Back-analysis) và phương pháp phần tử hữu hạn (FEM) áp dụng trên số liệu thực chứng của một công trình cụ thể tại Việt Nam.

Nhóm 4: Chủ đầu tư và các đơn vị tư vấn quản lý dự án, thẩm tra. Lợi ích: Có cơ sở khoa học để đánh giá mức độ an toàn của biện pháp thi công hố móng sâu trong đô thị, tối ưu hóa chi phí đầu tư cho hệ chống giữ tạm và kiểm soát rủi ro pháp lý liên quan đến lún nứt công trình lân cận.

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

Phương pháp phân tích ngược có ưu điểm vượt trội gì so với phương pháp tính toán xuôi truyền thống? Phương pháp tính toán xuôi chỉ dựa vào các chỉ tiêu cơ lý mẫu đất lấy từ hố khoan thí nghiệm trong phòng, vốn bị xáo động và sai lệch so với hiện trường. Phân tích ngược sử dụng trực tiếp chuyển vị thực tế của tường vây sâu 21m qua các chu kỳ đo Inclinometer để tìm ra bộ thông số đất $c', \varphi', E'$ thực tế, giúp độ chính xác của mô hình tăng trên 30%.

Tại sao luận văn chọn công trình President Place làm trường hợp nghiên cứu điển hình? President Place tại 93 Nguyễn Du, Quận 1 sở hữu quy mô tiêu biểu với 3 tầng hầm, độ sâu đào 13,5m, sử dụng tường vây D600 và thi công theo phương pháp Bottom-Up kết hợp 3 lớp shoring. Công trình có đầy đủ hệ thống quan trắc 6 ống Inclinometer và hồ sơ địa chất chuẩn mực, tạo nguồn dữ liệu thực chứng tin cậy cho bài toán phân tích ngược.

Sự khác nhau cơ bản giữa mô hình Mohr-Coulomb và Hardening Soil khi mô phỏng hố đào sâu là gì? Mô hình Mohr-Coulomb giả định đất có độ cứng không đổi và ứng xử dẻo thuần túy, thường ước tính chuyển vị đáy móng lớn hơn thực tế. Ngược lại, mô hình Hardening Soil phân biệt độ cứng khi gia tải và dỡ tải ($E_{ur} \approx 3E_{50}$), đồng thời tính đến sự tăng độ cứng theo độ sâu ứng suất, phản ánh chính xác hiện tượng trồi đáy hố đào sâu 13,5m.

Hệ thanh chống shoring H350 và H400 đóng vai trò gì trong việc kiểm soát chuyển vị tường? Ba lớp giằng chống đặt tại các cao độ -3,1m, -6,2m (thép H350) và -9,3m (thép H400) liên kết với cột Kingpost đóng vai trò là các gối tựa đàn hồi cưỡng bức. Hệ kết cấu này tiếp nhận áp lực đất chủ động, khống chế độ võng ngang của tường vây D600 không vượt quá ngưỡng 25 mm trong suốt quá trình đào đất.

Kết quả thông số đất của luận văn có thể áp dụng cho các công trình khác tại Thành phố Hồ Chí Minh không? Bộ thông số hiệu chỉnh ($c', \varphi', E_{ref}$) hoàn toàn có thể tham khảo trực tiếp cho các dự án xây dựng tầng hầm trong phạm vi địa tầng sét dẻo và cát pha thuộc khu vực Quận 1 và Quận 3 lân cận. Tuy nhiên, với mỗi dự án mới, kỹ sư vẫn cần thực hiện chu trình phân tích ngược kiểm chứng với số liệu quan trắc ban đầu.

Kết luận

Nghiên cứu của luận văn đã hoàn thành trọn vẹn mục tiêu khoa học và thực tiễn thông qua 5 kết luận và đóng góp cốt lõi:

  • Xây dựng thành công quy trình phân tích ngược chuyển vị tường vây D600 bằng phần mềm Plaxis 3D Foundation 1.6, thu hẹp độ sai lệch giữa mô phỏng và thực tế từ trên 35% xuống dưới 8%.
  • Xác định được bộ thông số biến dạng thực tế của đất nền khu vực Quận 1, Thành phố Hồ Chí Minh với mô đun $E_{ref}$ đạt mức 34.000 đến 38.000 kN/m2, cao hơn từ 50% đến 70% so với số liệu nén một trục truyền thống.
  • Làm sáng tỏ quy luật chuyển vị ngang của tường chắn sâu 21m dưới tác động của trình tự đào đất 4 giai đoạn và hiệu quả ngăn ngừa chuyển vị của 3 tầng giằng shoring H350/H400.
  • Khẳng định tính ưu việt của mô hình phần tử hữu hạn 3D trong việc mô tả hiệu ứng vòm không gian và tương tác đất - kết cấu của hố đào sâu 50m x 30m.
  • Đề xuất quy trình đồng bộ từ khảo sát, mô phỏng số đến quan trắc Inclinometer thực tế nhằm kiểm soát an toàn tuyệt đối cho các công trình ngầm đô thị.

Kế hoạch phát triển tiếp theo của nghiên cứu định hướng ứng dụng các thuật toán tối ưu hóa tự động (như thuật toán di truyền Genetic Algorithm hoặc mạng nơ-ron nhân tạo) kết hợp với công nghệ quan trắc thời gian thực IoT để tự động cập nhật thông số mô hình đất theo từng giờ thi công. Các kỹ sư địa kỹ thuật, nhà thầu xây dựng và nhà nghiên cứu quan tâm có thể khai thác phương pháp luận và dữ liệu thực chứng này để nâng cao độ an toàn và tối ưu hóa chi phí cho các dự án hố đào sâu trong tương lai.