Tổng quan về luận án

Công trình hố đào sâu tại các đô thị lớn Việt Nam đang đối mặt với nghịch lý kỹ thuật nghiêm trọng: kết quả dự báo chuyển vị tường chắn từ thiết kế thường lớn hơn đáng kể so với quan trắc thực tế. Tại dự án Nhiêu Lộc – Thị Nghè (TP. Hồ Chí Minh), mô hình Mohr-Coulomb dự báo chuyển vị lớn hơn quan trắc thực tế từ 15mm ở giai đoạn thi công thứ nhất; mô hình Hardening Soil cũng sai lệch từ 8mm trở lên ở cùng giai đoạn. Đây không phải sai số ngẫu nhiên mà là sai số hệ thống, bắt nguồn từ bộ thông số địa kỹ thuật đầu vào không phản ánh đúng trạng thái ứng suất thực tế của đất sau lưng tường chắn.

Research gap cốt lõi được xác định rõ ràng: vùng đất E (theo phân loại của Wong, 2009), tức vùng bên cạnh kết cấu chắn giữ nơi ứng suất thẳng đứng σ₁ giữ nguyên trong khi ứng suất ngang σ₃ giảm dần, là vùng "chưa được nghiên cứu đầy đủ" theo đánh giá của chính Wong (2009). Hầu hết phần mềm phổ biến – Plaxis, Geoslope, VEX, KCW – vẫn sử dụng tham số từ thí nghiệm ba trục thông thường (tăng σ₁, giữ nguyên σ₃), vốn không mô phỏng đúng đường ứng suất của vùng đất này.

Luận án tiến sĩ của Nguyễn Trường Huy (Viện Khoa học Công nghệ Xây dựng, 2015), chuyên ngành Kỹ thuật Xây dựng Công trình Ngầm, đặt ra và trả lời ba câu hỏi nghiên cứu cốt lõi: (1) Liệu thí nghiệm ba trục giảm ứng suất ngang có mô phỏng đúng trạng thái làm việc thực tế của đất bên cạnh hố đào? (2) Mô hình đất nền nào phản ánh chính xác quan hệ ứng suất–biến dạng theo đường ứng suất giảm σ₃? (3) Quy trình đồng bộ khảo sát–thí nghiệm–mô hình hóa nào phù hợp với hố đào sâu tại điều kiện đô thị Việt Nam?

Khung lý thuyết tích hợp ba trụ cột: lý thuyết trạng thái ứng suất trong đất theo Hans-Georg Kempfert và Berhane Geobreselassie (2006), mô hình đàn dẻo Lade (Lade và Jacobsen, 2002), và phương pháp phần tử hữu hạn với phần mềm LADEDEEP do tác giả tự lập trình. Phạm vi nghiên cứu bao trùm đất dính trong phạm vi ảnh hưởng đến 4 lần chiều sâu hố đào, điều kiện vách thẳng đứng, thi công nhanh tại khu đô thị mật độ cao, với dữ liệu thực địa từ ba công trình tại Hà Nội: Bệnh viện 108, Hoabinh Green City Minh Khai và đập thủy điện số 2.

Literature Review và Positioning

Nghiên cứu về hố đào sâu có lịch sử học thuật từ Terzaghi (1943) – người đầu tiên đề xuất tiêu chí phân loại hố đào nông/sâu dựa trên tỷ lệ chiều sâu/chiều rộng – đến Terzaghi & Peck (1967) và Peck et al. (1977) với ngưỡng định lượng 6m chiều sâu. Hai trường phái nghiên cứu quốc tế nổi bật gồm: nhóm tập trung phân tích ổn định và chuyển vị (Clough & Hansen, 1981; Mana, 1978; Mana & Clough, 1981; Ou & Shiau, 1998; Peck, 1969) và nhóm nghiên cứu mô hình đất nền và đường ứng suất (Bjerrum & Eide, 1956; Hashash, 1992; Bentler, 1998; Wong, 2009).

Tranh luận học thuật nổi bật tập trung vào hai quan điểm đối lập về độ chính xác của mô hình. Pearlmann et al. (2004) đã tổng kết nhiều công trình và kết luận: "Trong hầu hết các hiện trường đều cho thấy sự chuyển vị thực tế của tường chắn là nhỏ hơn dự báo theo tính toán" – điều này ngụ ý các mô hình hiện hành thiên về dự báo quá mức. Ngược lại, Ou (2006) và Wong (2009) nhấn mạnh bản chất phụ thuộc đường ứng suất của đất nền, cho rằng vấn đề không nằm ở mô hình mà ở sự không tương thích giữa đường ứng suất thí nghiệm và đường ứng suất thực tế.

Trong nước, các nghiên cứu của Đỗ Đình Đức (2002) và Chu Tuấn Hạ (2011) tập trung vào điều kiện địa kỹ thuật địa phương nhưng chưa giải quyết được vấn đề mô phỏng đường ứng suất giảm σ₃. Tổng kết hơn 160 sự cố hố đào tại Trung Quốc cho thấy nguyên nhân do thiết kế chiếm 46% và do thi công chiếm 41,5%, trong đó "thiếu hoặc không chính xác về thông tin nền đất" chiếm đa số. Bối cảnh Việt Nam có các sự cố điển hình: sập nhà Viện KHXH miền Nam do thi công Cao ốc Pacific (Nguyễn Thị Minh Khai, 2007) và nghiêng chung cư Cosaco do thi công Cao ốc Saigon Residences cùng năm – minh chứng trực quan cho khoảng trống nghiên cứu này.

So sánh quốc tế: trong khi nghiên cứu tại Trung Quốc (tổng kết của Ou, 2006) tập trung vào hệ chống đỡ và kiểm soát chuyển vị, và nghiên cứu của Hashash (1992) tại MIT nhấn mạnh mô hình phi tuyến cho đất sét Boston, luận án của Nguyễn Trường Huy định vị độc đáo ở nút thắt chưa được giải quyết: kết hợp phát triển thiết bị thí nghiệm mới với mô hình đất nền phù hợp cho đường ứng suất giảm σ₃ trong điều kiện đất dính đô thị Việt Nam.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Đóng góp lý thuyết quan trọng nhất là việc mở rộng và cải tiến mô hình Lade (Lade & Jacobsen, 2002) thành mô hình Lade cải tiến tích hợp bốn thành phần cấu trúc mới: (1) đặc trưng đàn hồi phi tuyến, (2) mặt phá hoại cải tiến trong không gian ứng suất chính, (3) hàm thế năng biến dạng dẻo điều chỉnh theo điều kiện giảm σ₃, và (4) mặt chảy dẻo phù hợp với cả đất cát giãn nở và đất sét không giãn nở thể tích.

Luận án chứng minh rằng mô hình Cam-Clay cải tiến (Roscoe & Burland, 1968) không mô tả đúng ứng xử dỡ tải khi giảm σ₃. Phân tích toán học cho thấy khi p'/p₀' < 0,11, đất nền theo mô hình Cam-Clay cải tiến trở về trạng thái đàn hồi với mô đun tiếp tuyến tăng lên – điều này mâu thuẫn với thực nghiệm vì mô đun tiếp tuyến phải luôn giảm khi độ lệch ứng suất tăng. Tương tự, mô hình Hyperbol cho thấy vi phân mô đun đàn hồi tiếp tuyến theo σ₃ có thể cho giá trị dương (đất trở nên cứng hơn khi dỡ tải ngang) trong điều kiện 2c·cosφ + 2σ₁·sinφ < 0 – một nghịch lý vật lý.

Về mặt khung phân tích, công trình tích hợp ba lý thuyết: (i) lý thuyết dẻo ứng suất tổng của Lade (kiểm soát mặt chảy và mặt phá hoại), (ii) lý thuyết đường ứng suất theo Stroh (1974) và Wong (2009) (phân vùng ứng xử đất), và (iii) phương pháp phần tử hữu hạn biến dạng phẳng (mô hình hóa quá trình thi công). Boundary conditions được xác định rõ: đất dính, hố đào vách thẳng đứng, chiều sâu > 6m, điều kiện không thoát nước tức thời, phạm vi ảnh hưởng ≤ 4H (H = chiều sâu hố đào).

Khung phân tích độc đáo

Khung phân tích mới xây dựng chuỗi liên tục: Trạng thái ứng suất thực tế (vùng E theo Wong, 2009) → Thí nghiệm ba trục giảm σ₃ (RTC path) → Tham số mô hình Lade cải tiến → Phần mềm LADEDEEP (PTHH) → Dự báo chuyển vị ngang. Tính mới nằm ở sự đồng bộ giữa đường ứng suất thí nghiệm và đường ứng suất trong thực tế – điều mà các khung phân tích hiện hành (dùng CTC path cho vùng E) đã bỏ qua.

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

Thiết kế nghiên cứu

Triết lý nghiên cứu định vị theo chủ nghĩa hiện thực phê phán (critical realism): thừa nhận tồn tại khách quan của cơ chế vật lý (trạng thái ứng suất–biến dạng) trong khi nhận thức về nó luôn bị giới hạn bởi công cụ đo lường. Thiết kế nghiên cứu là hỗn hợp tuần tự (sequential mixed methods): nghiên cứu thực nghiệm phòng thí nghiệm → mô hình hóa số học → kiểm chứng thực địa đa công trình. Đây là thiết kế đa tầng với ba cấp độ: mẫu đất đơn lẻ (tầng vi mô), hố đào đơn công trình (tầng trung mô), và tổng hợp đa công trình (tầng vĩ mô).

Quy trình nghiên cứu rigorous

Thiết bị thí nghiệm: Máy ba trục giảm ứng suất ngang được thiết kế và chế tạo đặc biệt tại Trường Đại học Kiến trúc Hà Nội, với bộ phận giảm áp buồng có kiểm soát tốc độ và hệ đo biến dạng thể tích tích hợp. Thiết bị được so sánh với thiết bị GDS LVDT (GDS Instruments) và hệ thống cảm biến tiệm cận hiện đại nhưng với chi phí thực tế phù hợp. Phụ lục A của luận án mô tả chi tiết sáu thành phần: buồng ba trục, thiết bị duy trì áp lực nước, thiết bị tạo lực dọc trục, thiết bị giảm áp buồng, thiết bị đo biến dạng/áp suất, và thiết bị ghi chép số liệu.

Quy trình thí nghiệm theo sơ đồ RTC (Reduced Triaxial Compression): (1) Áp đặt trạng thái ứng suất thủy tĩnh ban đầu; (2) Áp dụng ứng suất lệch ban đầu để mô phỏng điều kiện Ko thực tế (phù hợp với đường TX_116 từ nghiên cứu tham chiếu); (3) Giữ σ₁ không đổi; (4) Giảm dần σ₃ đến khi phá hoại với tốc độ kiểm soát. Đây là sự khác biệt quan trọng so với thí nghiệm CTC thông thường (tăng σ₁, giữ nguyên σ₃).

Tam giác hóa (triangulation): dữ liệu (số liệu thí nghiệm phòng + quan trắc hiện trường + báo cáo khảo sát địa chất), phương pháp (thí nghiệm nén ba trục + PTHH + phân tích hồi quy), và lý thuyết (mô hình Lade + Cam-Clay + Hyperbol để so sánh).

Phần mềm: LADEDEEP (tự lập trình, mã nguồn tại Phụ lục C), PLAXIS (so sánh), ngôn ngữ lập trình cho PTHH với phần tử tứ giác đẳng tham số, phần tử thanh dàn, phần tử thanh dầm, phần tử tiếp xúc và phương pháp giải phi tuyến lặp Newton-Raphson.

Data và phân tích

Ba công trình kiểm chứng với đặc điểm khác nhau: (1) Bệnh viện 108 – dữ liệu chuyển vị đỉnh tường theo nhiều giai đoạn đào, so sánh giữa Lade cải tiến và Mohr-Coulomb; (2) Hoabinh Green City Minh Khai – hố đào nhiều giai đoạn, so sánh ở giai đoạn đào thứ 4 với ba mô hình; (3) Đập thủy điện số 2 – điều kiện địa chất khác biệt, kiểm chứng tính tổng quát của phương pháp. Phân tích ma trận đàn dẻo [Dp] và ma trận độ cứng tổng thể [K]s được thực hiện qua giải hệ phương trình phi tuyến theo từng giai đoạn thi công, mô phỏng trình tự đào thực tế.

Phát hiện đột phá và implications

Những phát hiện then chốt

Phát hiện 1: Sự khác biệt định lượng giữa hai đường ứng suất. Thí nghiệm ba trục giảm σ₃ cho quan hệ ứng suất–biến dạng khác biệt rõ ràng so với thí nghiệm ba trục thông thường tăng σ₁. Vòng tròn Mohr từ hai loại thí nghiệm trên cùng mẫu đất cho đường bao phá hoại không trùng nhau, khẳng định tính phụ thuộc đường ứng suất của các thông số cường độ.

Phát hiện 2: Mô hình Lade cải tiến vượt trội. Tại Bệnh viện 108 và Hoabinh Green City, kết quả mô hình Lade cải tiến với tham số từ thí nghiệm RTC cho chuyển vị đỉnh tường sát hơn quan trắc thực tế so với mô hình Mohr-Coulomb và Hardening Soil. Mô hình Mohr-Coulomb dự báo chuyển vị lớn hơn quan trắc thực tế (tương tự xu hướng đã ghi nhận tại dự án Nhiêu Lộc – Thị Nghè: sai lệch từ 15mm ở giai đoạn 1).

Phát hiện 3: Trạng thái ứng suất lệch ban đầu là yếu tố quyết định. Thí nghiệm TX_116 (có ứng suất lệch ban đầu trước khi giảm σ₃) cho kết quả phù hợp hơn thực tế so với TX_123 (không có ứng suất lệch ban đầu), xác nhận rằng điều kiện Ko trong đất thực tế phải được mô phỏng trong thí nghiệm.

Phát hiện 4: Hàm thế năng dẻo phân biệt theo loại đất. Tham số hàm thế năng dẻo (ψ₁, ψ₂) xác định được khác nhau có ý nghĩa thống kê giữa đất cát giãn nở thể tích và đất sét không giãn nở thể tích, cung cấp cơ sở cho việc phân loại ứng xử dẻo theo loại đất.

Phát hiện 5: Mô hình Cam-Clay cải tiến thất bại có hệ thống ở điều kiện giảm σ₃. Phân tích toán học (phương trình 1.6 trong luận án) chứng minh rằng khi p'/p₀' < 0,11 theo mô hình Cam-Clay cải tiến, mô đun tiếp tuyến tăng lên khi giảm σ₃ – nghịch lý với thực nghiệm. Đây là phát hiện có giá trị cảnh báo cho cộng đồng kỹ thuật.

Implications đa chiều

Lý thuyết: Kết quả mở rộng lý thuyết đàn dẻo của Lade sang điều kiện dỡ tải ngang, đồng thời cung cấp bằng chứng thực nghiệm về giới hạn áp dụng của mô hình Cam-Clay cải tiến và Hyperbol trong bài toán hố đào. Phương pháp: Quy trình thí nghiệm RTC có thể được áp dụng với thiết bị ba trục thông thường qua cải tiến bộ phận giảm áp buồng – điều này có nghĩa chi phí triển khai thấp. Thực tiễn: Tư vấn thiết kế có thể áp dụng ngay bộ tham số Lade cải tiến từ thí nghiệm RTC vào phần mềm LADEDEEP hoặc tích hợp vào PLAXIS qua user-defined model. Chính sách: Kết quả hỗ trợ việc cập nhật TCVN 9363-2012 và TCXDVN 194:2006 với yêu cầu bổ sung về loại thí nghiệm ba trục phù hợp với bài toán hố đào sâu.

Limitations và Future Research

Bốn giới hạn được thừa nhận minh bạch: (1) Phạm vi loại đất hẹp – nghiên cứu tập trung đất dính có thể lấy mẫu ba trục; đất cát và đất hỗn hợp chưa được kiểm chứng đầy đủ. (2) Điều kiện thoát nước – phân tích chủ yếu cho điều kiện thi công nhanh (không thoát nước tức thời); ảnh hưởng cố kết dài hạn chưa được mô hình hóa. (3) Số lượng công trình kiểm chứng – ba công trình, cả ba tại Hà Nội, hạn chế khả năng tổng quát hóa cho điều kiện địa chất TP. Hồ Chí Minh và Đà Nẵng. (4) Mô hình 2D biến dạng phẳng – không mô phỏng được hiệu ứng góc hố đào và điều kiện 3D.

Agenda nghiên cứu tương lai gồm năm hướng: (i) Mở rộng thí nghiệm RTC cho đất rời (cát Hà Nội, cát Đà Nẵng) với dải σ₃ rộng hơn; (ii) Phát triển phiên bản 3D của phần mềm LADEDEEP; (iii) Tích hợp mô hình cố kết để phân tích dài hạn sau thi công; (iv) Kiểm chứng tại các hố đào sâu tuyến metro Nhổn–Ga Hà Nội và Bến Thành–Suối Tiên; (v) Nghiên cứu hệ số Ko thực tế tại các vùng địa chất đặc thù và tác động đến tham số mô hình Lade cải tiến.

Tác động và ảnh hưởng

Về học thuật, công trình lần đầu tiên ở Việt Nam cung cấp quy trình đồng bộ khép kín: thí nghiệm RTC → tham số Lade cải tiến → PTHH → dự báo chuyển vị, có thể được tái lập và mở rộng bởi các nhóm nghiên cứu khác. Mã nguồn phần mềm LADEDEEP (Phụ lục C) được công bố đầy đủ, tạo nền tảng cho các phát triển tiếp theo.

Về ngành xây dựng, kết quả trực tiếp ảnh hưởng đến các dự án hố đào sâu trong hai tuyến đường sắt đô thị đang thi công (Nhổn–Ga Hà Nội và Bến Thành–Suối Tiên), các tổ hợp thương mại–văn phòng có nhiều tầng hầm (mô hình Lotte Center 5 tầng ngầm, D'.Palais de Louis 4 tầng hầm) và các trạm bơm ngầm trong dự án chống ngập đô thị.

Về chính sách, nghiên cứu cung cấp cơ sở khoa học cho việc bổ sung yêu cầu thí nghiệm đường ứng suất vào tiêu chuẩn khảo sát địa kỹ thuật quốc gia, giảm tỷ lệ sự cố hố đào hiện đang ở mức đáng lo ngại (nguyên nhân thiết kế 46%, thi công 41,5% theo thống kê Trung Quốc áp dụng tham chiếu cho Việt Nam).

Đối tượng hưởng lợi

Nghiên cứu sinh và học giả nhận được một gap rõ ràng – vùng ứng suất E theo phân loại Wong (2009) – với phương pháp giải quyết có thể tiếp tục mở rộng sang điều kiện đất khác, khu vực khác, và bài toán 3D. Kỹ sư địa kỹ thuật cấp cao được cung cấp quy trình thí nghiệm RTC thực thi được trên thiết bị hiện có qua cải tiến đơn giản. Đơn vị tư vấn thiết kế có thể tích hợp phần mềm LADEDEEP hoặc tham số Lade cải tiến vào quy trình kiểm tra thiết kế hố đào nhằm giảm hệ số dự phòng dư thừa và tối ưu chi phí kết cấu chắn giữ. Cơ quan quản lý nhà nước (Bộ Xây dựng, Bộ GTVT) có nền tảng để cập nhật TCVN về khảo sát địa kỹ thuật hố đào sâu, đặc biệt cho các dự án metro đô thị.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất? Cải tiến và mở rộng mô hình Lade (Lade & Jacobsen, 2002) với bốn thành phần mới – đặc biệt là hàm thế năng biến dạng dẻo điều chỉnh theo điều kiện giảm σ₃ – tạo thành mô hình Lade cải tiến với bộ tham số xác định trực tiếp từ thí nghiệm RTC, thay vì dùng tham số từ thí nghiệm CTC như tất cả mô hình hiện hành.

2. Innovation về phương pháp? So với Hashash (1992) dùng mô hình MIT-E3 cần thí nghiệm phức tạp, và so với Ou (2006) dùng Hardening Soil trong PLAXIS với thông số CTC thông thường, công trình này phát triển thiết bị RTC có thể chế tạo tại chỗ từ máy ba trục hiện có, với quy trình thí nghiệm được tiêu chuẩn hóa và phần mềm LADEDEEP mã nguồn mở.

3. Phát hiện bất ngờ nhất? Mô hình Cam-Clay cải tiến – vốn được coi là mô hình tiên tiến cho đất sét – thất bại có hệ thống trong điều kiện giảm σ₃ khi p'/p₀' < 0,11, dự báo đất trở nên cứng hơn khi dỡ tải ngang. Phân tích toán học (phương trình 1.6) cho thấy đây là giới hạn cấu trúc của mặt chảy ellipsoid, không thể khắc phục bằng hiệu chỉnh tham số.

4. Giao thức tái lập? Có. Phụ lục A mô tả đầy đủ cấu tạo thiết bị; Phụ lục B trình bày cơ sở lý thuyết PTHH; Phụ lục C cung cấp mã nguồn LADEDEEP – đủ để các phòng thí nghiệm khác tái lập nghiên cứu.

5. Agenda nghiên cứu 10 năm? Có thể định hình từ năm hướng đã đề xuất: mở rộng sang đất rời → phát triển 3D → tích hợp cố kết → kiểm chứng metro → nghiên cứu Ko vùng → tích hợp vào tiêu chuẩn quốc gia, tạo lộ trình từ nghiên cứu phòng thí nghiệm đến chuẩn hóa thực hành kỹ thuật.

Kết luận

Nghiên cứu này tạo ra sáu đóng góp cụ thể có thể đo lường:

(1) Thiết lập và kiểm chứng thí nghiệm ba trục giảm ứng suất ngang (RTC) như một phương pháp đặc trưng hóa đất nền phù hợp với bài toán hố đào sâu, giải quyết gap nghiên cứu mà Wong (2009) đã chỉ ra nhưng chưa có lời giải.

(2) Phát triển mô hình Lade cải tiến với bộ tham số xác định từ thí nghiệm RTC, chứng minh vượt trội so với Mohr-Coulomb và Hardening Soil trong ba công trình kiểm chứng thực tế tại Hà Nội.

(3) Chứng minh bằng toán học và thực nghiệm giới hạn áp dụng của mô hình Cam-Clay cải tiến và Hyperbol trong điều kiện dỡ tải ngang – đóng góp cảnh báo có giá trị cho cộng đồng kỹ thuật quốc tế.

(4) Lập trình phần mềm LADEDEEP theo phương pháp PTHH với mã nguồn công bố đầy đủ, hỗ trợ cộng đồng tái lập và mở rộng nghiên cứu.

(5) Đề xuất quy trình đồng bộ ba bước: khảo sát địa kỹ thuật → thí nghiệm RTC → mô hình hóa PTHH với Lade cải tiến, tạo nền tảng cho cập nhật tiêu chuẩn TCVN về hố đào sâu.

(6) Mở ra ba luồng nghiên cứu mới: (a) ứng xử đất theo đường ứng suất dỡ tải ngang cho các loại đất khác tại Việt Nam; (b) phát triển thiết bị ba trục giảm σ₃ thế hệ tiếp theo với đo biến dạng ngang trực tiếp; (c) tích hợp mô hình Lade cải tiến vào các nền tảng phần mềm thương mại như PLAXIS qua giao diện user-defined soil model.

Vị thế quốc tế của công trình nằm ở điểm tiếp nối độc đáo: trong khi các nghiên cứu MIT (Hashash, 1992), NUS Singapore (Ou, 2006) và HKUST (Wong, 2009) phân tích vấn đề ở cấp độ lý thuyết, công trình này cung cấp giải pháp thực thi với chi phí hợp lý cho điều kiện phòng thí nghiệm đang phát triển – một đóng góp đặc biệt có ý nghĩa cho các quốc gia đang đô thị hóa nhanh tại Đông Nam Á.