Tổng quan về luận án

Xây dựng công trình hạ tầng giao thông và đê đập trên nền đất sét mềm bão hòa nước có bề dày lớn là thách thức địa kỹ thuật hàng đầu tại khu vực Đồng bằng sông Cửu Long (ĐBSCL) và vùng phụ cận TP. Hồ Chí Minh. Luận án tiến sĩ kỹ thuật chuyên ngành Địa kỹ thuật xây dựng (Mã số: 62.11) của nghiên cứu sinh Lâm Ngọc Quí, thực hiện tại Trường Đại học Bách khoa – Đại học Quốc gia TP. Hồ Chí Minh dưới sự hướng dẫn khoa học của PGS.TS. Bùi Trường Sơn, mang tên: "Nghiên cứu phương pháp đánh giá độ lún cố kết nền đất yếu theo sơ đồ hai chiều có xét áp lực nước lỗ rỗng phụ thuộc ứng suất trung bình".

Research gap cốt lõi xuất phát từ sự bất cập của lý thuyết cố kết một chiều cổ điển (Terzaghi, 1925), vốn giả định tải trọng phân bố đều vô hạn, áp lực nước lỗ rỗng (ALNLR) thặng dư ban đầu ($u_0$) sinh ra tức thời và bằng đúng gia số ứng suất đứng ($\Delta \sigma_z$), bỏ qua chuyển vị ngang và biến dạng hình học trong điều kiện diện gia tải hẹp dạng tuyến. Luận án đặt ra ba câu hỏi nghiên cứu trọng tâm:

  1. RQ1: Quy luật phát sinh và tiêu tán ALNLR thặng dư thực tế trong đất sét mềm bão hòa nước theo thời gian và ứng suất tác dụng diễn ra như thế nào?
  2. RQ2: Mô hình giải tích cố kết hai chiều (2D plane strain) nào có khả năng tích hợp chính xác trường ALNLR phụ thuộc ứng suất trung bình ($\sigma_m$) và quy luật thấm phi tuyến?
  3. RQ3: Độ chính xác của mô hình phần tử đơn vị 2D khi áp dụng dự báo độ lún, lún lệch theo thời gian cho nền tự nhiên và nền gia cố bấc thấm/giếng cát kết hợp gia tải trước so với số liệu quan trắc thực tế đạt mức độ nào?

Hệ thống giả thuyết khoa học bao gồm: $H_1$ – ALNLR thặng dư cực đại ($u_{max}$) trong sét mềm không bằng $\Delta \sigma_z$ mà phụ thuộc tuyến tính vào ứng suất trung bình $\sigma_m = (\sigma_1 + \sigma_2 + \sigma_3)/3$; $H_2$ – Thời gian đạt $u_{max}$ có độ trễ do tính nhớt của cốt đất và bọt khí hòa tan (hiệu ứng ba pha), không triệt tiêu ngay khi bắt đầu thoát nước; $H_3$ – Tách độ lún thành biến dạng thể tích ($S_v$) và biến dạng trượt/hình dạng ($S_s$) trong không gian 2D cho phép triệt tiêu sai số dự báo lún lệch tại tâm và mép đắp. Phạm vi nghiên cứu bao gồm thực nghiệm chuyên sâu trên mẫu sét mềm TP.HCM (Quận 7) và ĐBSCL, đồng thời kiểm chứng tại 2 đại dự án: Quốc lộ 1A mở rộng đoạn Vĩnh Long – Cần Thơ và đường đầu cầu Cao Lãnh (Đồng Tháp). Đóng góp mang tính đột phá thể hiện qua việc lượng hóa biến dạng thể tích chiếm 80–95% tổng độ lún, kiểm soát chuyển vị ngang (chiếm 5–20%) và tối ưu hóa 10–15% khối lượng bù lún công trình đắp.

+-----------------------------------------------------------------------------------+
|                        KHUNG DẪN LUẬN ÁN TIẾN SĨ (LÂM NGỌC QUÍ)                   |
+-----------------------------------------------------------------------------------+
|  THỰC TRẠNG & RESEARCH GAP:                                                       |
|  - Terzaghi 1D: u0 = Δσz, tải rộng vô hạn, bỏ qua chuyển vị ngang (S_s).          |
|  - Thực tế đồng bằng: Chiều dày sét mềm lớn, tải đắp hẹp hình thang gây lún lệch. |
+-----------------------------------------------------------------------------------+
                                         │
                                         ▼
+-----------------------------------------------------------------------------------+
|  ĐỘT PHÁ THỰC NGHIỆM:                                                             |
|  - Hộp nén cố kết cải tiến (H = 2, 4, 6 cm) có sensor đo ALNLR đáy mẫu.           |
|  - Phát hiện u_max ≈ σ_m (nhỏ hơn Δσ_z); trễ 50-70 phút mới đạt đỉnh.             |
+-----------------------------------------------------------------------------------+
                                         │
                                         ▼
+-----------------------------------------------------------------------------------+
|  MÔ HÌNH HÓA & KIỂM CHỨNG:                                                        |
|  - Thiết lập phương trình cố kết 2D bài toán phẳng: u = f(σ_m, t, x, z).           |
|  - Unit cell biến dạng phẳng cho giếng cát/PVD kết hợp tải trước (11 giai đoạn).  |
|  - Kiểm chứng: QL1A Vĩnh Long và Cầu Cao Lãnh (Km 8+100) khớp quan trắc R² > 0.95.|
+-----------------------------------------------------------------------------------+

Literature Review và Positioning

Lịch sử phát triển cơ học đất ghi nhận nền tảng từ lý thuyết cố kết một chiều của Terzaghi (1925), phát triển tiếp nối bởi Biot (1941) với lý thuyết cố kết ba chiều đàn hồi nhiều pha kết hợp tương tác giữa khung cốt đất và chất lưu lỗ rỗng. Barron (1948) và Hansbo (1979, 1981) giải quyết bài toán cố kết thấm xuyên tâm đối xứng trục quanh bấc thấm (PVD) và giếng cát có xét đến hiệu ứng xáo động (smear zone, $k_s$, $r_s$) và sức cản giếng (well resistance, $q_w$). Tuy nhiên, tranh luận học thuật kéo dài giữa hai trường phái:

  • Trường phái 1D cổ điển và đàn nhớt: Biến dạng lún thuần túy do ép co thể tích dọc trục z, ALNLR ban đầu tiêu tán ngay lập tức (Barden, 1965; Mesri, 1973; Yin & Graham, 1996; Tsytovich, 1976).
  • Trường phái Biến dạng phẳng 2D và Chuyển đổi thấm: Nhấn mạnh chuyển vị ngang đóng vai trò sống còn trong nền đất đắp hẹp. Các nghiên cứu chuyển đổi hệ số thấm từ đối xứng trục ($k_{ax}$) sang bài toán phẳng ($k_{pl}$, $k_{hp}$) của Hird (1992, 1995), Bergado & Long (1994), Chai (1995), Indraratna & Redana (1997), Lin et al. (2000) đã mở đường cho mô hình hóa 2D nhưng phần lớn vẫn mặc định $u_0 = \Delta \sigma_z$ cục bộ.
Tiêu chí So sánh Terzaghi (1925) / Barron (1948) Indraratna & Redana (1997) / Lin et al. (2000) Mô hình Luận án (Lâm Ngọc Quí, 2022)
Sơ đồ không gian 1D thẳng đứng / Đối xứng trục 3D 2D biến dạng phẳng quy đổi tương đương 2D biến dạng phẳng theo phân tố đơn vị
Áp lực nước lỗ rỗng ban đầu ($u_0$) $u_0 = \Delta \sigma_z$ (Tức thời, đẳng hướng) $u_0 = \Delta \sigma_z$ hoặc theo hàm ứng suất đứng $u_{max} \approx \sigma_m = (\sigma_x + \sigma_y + \sigma_z)/3$
Động học ALNLR thặng dư Đạt cực đại tại $t=0$, tiêu tán đơn điệu Tiêu tán đơn điệu theo thời gian thấm Gia tăng từ từ, đạt đỉnh sau 50–70 phút rồi mới tiêu tán
Thành phần độ lún $S = S_v$ (Thuần túy biến dạng thể tích) $S = S_v + S_s$ (Mô phỏng phần tử hữu hạn) Phân tích tách biệt: $S_v$ (80–95%) và $S_s$ (5–20%)
Quy luật thấm & nén Tuyến tính ($k, a_v = const$) Thường giả định tuyến tính từng khoảng Phi tuyến: $e - \ln(p')$, $k(\sigma')$ hàm số mũ

So sánh với hai công trình quốc tế kinh điển:

  1. Nghiên cứu hiện trường Lilla Mellösa của Viện Địa kỹ thuật Thụy Điển (Hansbo, 1960; Larsson, 2002): Quan trắc đê đắp trên sét mềm bão hòa nước suốt 57 năm (1945–2002) chỉ ra rằng sau 1 năm gia tải, ALNLR thặng dư phân bố dưới độ sâu 5m xấp xỉ tải trọng tác dụng, và có tới 0.2m độ lún sinh ra thuần túy do chuyển dịch ngang của đất nền ra ngoài chân taluy.
  2. Nghiên cứu xử lý nền cảng Thị Vải (Pham Van Long et al., 2013): Phân tích ngược 1D cho thấy tỷ số xáo trộn $d_s/d_m$ và biến thiên hệ số thấm $k_h/k_v$ làm sai lệch hệ số cố kết $C_h$ và hệ số từ biến thứ cấp $C_\alpha$.

Định vị của luận án nằm tại điểm giao thoa giữa cơ học môi trường nhiều pha (poromechanics) và địa kỹ thuật công trình đắp thực hành: Khắc phục triệt để giả thiết $u_0 = \Delta \sigma_z$ bằng việc tích hợp trường ứng suất trung bình $\sigma_m$ vào hệ phương trình vi phân biến dạng phẳng 2D, giải quyết hoàn chỉnh bài toán lún lệch và tiến độ tiêu tán áp lực nước thực tế.

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

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

Luận án mở rộng và định hình lại hệ tiên đề của lý thuyết cố kết Terzaghi – Biot:

  1. Bác bỏ tiên đề ALNLR tức thời bằng tải trọng đứng: Luận án chứng minh trong môi trường sét mềm bão hòa, phân bố ALNLR tức thời không thể tách rời trạng thái ứng suất đa trục. Gia số áp lực nước lỗ rỗng thặng dư liên hệ trực tiếp với bất biến ứng suất bậc nhất thông qua ứng suất trung bình: $$\sigma_m = \frac{\sigma_1 + \sigma_2 + \sigma_3}{3} = \frac{\sigma_z + 2\sigma_x}{3}$$ Vì $\sigma_x < \sigma_z$, giá trị $\sigma_m$ luôn nhỏ hơn đáng kể so với gia số ứng suất đứng $\Delta \sigma_z$. Do đó, việc dùng $\Delta \sigma_z$ tính toán trong bài toán 1D dẫn đến ước lượng vượt quá (overestimation) ALNLR ban đầu và đánh giá sai quy luật tiêu tán.
  2. Phát hiện hiệu ứng trễ pha của cực đại ALNLR: Thực nghiệm xác lập định đề $u(t)$ không đạt đỉnh tại $t=0$ mà hình thành đường cong có cực trị, leo thang từ từ và đạt đỉnh sau 50–70 phút do sự nén ép của bọt khí hòa tan (hệ số hòa tan Henry $H$) và từ biến trượt của cốt đất nhớt.
  3. Mô hình cấu trúc phân rã độ lún 2D: $$S(t) = S_s(t) + S_v(t) = \int_0^h \frac{\sigma_z - \sigma_m}{2G} dz + \int_0^h \frac{\sigma_m - u(t)}{K} dz$$ Trong đó $G = \frac{E_0}{2(1+\nu)}$ là mô đun trượt, $K = \frac{E_0}{3(1-2\nu)}$ là mô đun nén thể tích, $E_0$ là mô đun tổng biến dạng, $\nu$ là hệ số Poisson.
        CƠ CHẾ TRUYỀN TẢI TRỌNG VÀ BIẾN DẠNG THEO SƠ ĐỒ 2D
 
         Tải trọng ngoài hình thang q(x) (Nền đường đắp)
  ─────────────────────────────────────────────────────────────
                             │
                             ▼
  ┌───────────────────────────────────────────────────────────┐
  │         TRƯỜNG ỨNG SUẤT 2D TRONG NỀN ĐẤT YẾU              │
  │     Ứng suất đứng σ_z(x,z) & Ứng suất ngang σ_x(x,z)       │
  └──────────────────────────┬────────────────────────────────┘
                             │
            ┌────────────────┴────────────────┐
            ▼                                 ▼
  ┌───────────────────┐             ┌───────────────────┐
  │ Ứng suất lệch     │             │ Ứng suất trung    │
  │ q = σ_z - σ_x     │             │ bình σ_m          │
  └─────────┬─────────┘             └─────────┬─────────┘
            │                                 │
            ▼                                 ▼
  ┌───────────────────┐             ┌───────────────────┐
  │ Biến dạng hình    │             │ Biến dạng thể     │
  │ dạng / Cắt (S_s)  │             │ tích (S_v)        │
  │ (Chiếm 5 - 20%)   │             │ (Chiếm 80 - 95%)  │
  └─────────┬─────────┘             └─────────┬─────────┘
            │ Chuyển vị ngang                 │ Tiêu tán ALNLR
            │ trượt trồi mép                  │ u(t) -> 0
            └────────────────┬────────────────┘
                             ▼
  ┌───────────────────────────────────────────────────────────┐
  │  TỔNG ĐỘ LÚN VÀ LÚN LỆCH THEO THỜI GIAN: S(t) = S_s + S_v │
  └───────────────────────────────────────────────────────────┘

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

Hệ phương trình vi phân cố kết biến dạng phẳng liên tục đa pha được thiết lập từ nguyên lý bảo toàn khối lượng và định luật Darcy: $$\frac{k_x}{\gamma_w} \frac{\partial^2 u}{\partial x^2} + \frac{k_z}{\gamma_w} \frac{\partial^2 u}{\partial z^2} = m_v \frac{\partial u}{\partial t} - m_v \frac{\partial \sigma_m}{\partial t} + \alpha_c \frac{\partial^2 u}{\partial z \partial t}$$ Trong đó:

  • $m_v$: Hệ số nén thể tích tương đối của cốt đất ($m_v = \frac{a_0}{1+e_0}$).
  • $\alpha_c$: Hệ số từ biến nhớt cốt đất.
  • Điều kiện biên thấm và tải trọng phi đối xứng: $u(x, 0, t) = 0$ tại mặt thoát nước; $\frac{\partial u}{\partial x} = 0$ tại trục đối xứng tâm; $\sigma_m(x,z)$ xác định theo lời giải giải tích đàn hồi cho tải hình thang có độ dốc mái taluy ($1:m$).

Boundary Conditions (Điều kiện biên giới hạn): Mô hình áp dụng tối ưu cho nền sét mềm đồng nhất hoặc phân lớp có bề dày $H \ge 10$ m, hệ số quá cố kết $OCR \le 1.2$, bề rộng tải đắp $B \le 40$ m, chịu tải trọng tĩnh phân bố hình thang.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành trên nền tảng triết học Thực chứng Phản biện (Critical Realism) kết hợp Cơ học môi trường liên tục (Continuum Mechanics). Thiết kế nghiên cứu đa cấp độ (multi-level design) bao gồm:

  1. Cấp độ vi mô - thực nghiệm trong phòng (Laboratory Testing): Chế tạo buồng nén cải tiến kiểm soát ALNLR đáy mẫu.
  2. Cấp độ trung mô - phân tích giải tích & mô hình hóa (Analytical & Unit Cell Modeling): Xây dựng thuật toán phân chia lớp phân tố và giải số hệ phương trình cố kết 2D.
  3. Cấp độ vĩ mô - kiểm chứng hiện trường (Field Validation): Phân tích đối chiếu trên các công trình thực tế quan trắc dài hạn.
                           THIẾT KẾ NGHIÊN CỨU ĐA CẤP ĐỘ
 
  ┌───────────────────────────────────────────────────────────────────────────┐
  │ LEVEL 1: THỰC NGHIỆM TRONG PHÒNG (Micro-scale)                           │
  │ - Chế tạo buồng nén Oedometer tích hợp sensor ALNLR đáy.                  │
  │ - Mẫu đất: M1 -> M7 sét mềm TP.HCM & ĐBSCL (Đường kính D=60mm, H=2,4,6cm).│
  │ - Cấp tải: 50 -> 100 -> 200 -> 400 -> 800 kPa. Thoát nước 1D & 2D.       │
  └─────────────────────────────────────┬─────────────────────────────────────┘
                                        │
                                        ▼
  ┌───────────────────────────────────────────────────────────────────────────┐
  │ LEVEL 2: PHÁT TRIỂN MÔ HÌNH TOÁN HỌC 2D (Meso-scale)                      │
  │ - Phương trình cố kết 2D liên kết: u(x,z,t) gắn với ứng suất trung bình.  │
  │ - Mô hình phần tử đơn vị (Unit Cell) phẳng cho giếng cát/bấc thấm.        │
  │ - Tích hợp phi tuyến e - ln(p') và quy luật thấm biến thiên k(σ').        │
  └─────────────────────────────────────┬─────────────────────────────────────┘
                                        │
                                        ▼
  ┌───────────────────────────────────────────────────────────────────────────┐
  │ LEVEL 3: KIỂM CHỨNG THỰC ĐỊA QUAN TRẮC (Macro-scale)                      │
  │ - Case 1: QL1A Vĩnh Long - Cần Thơ (Nền tự nhiên, đắp cao 1.5 - 2.5m).     │
  │ - Case 2: Đường dẫn Cầu Cao Lãnh (Gia cố giếng cát d=0.4m, gia tải 11 đợt)│
  │ - Đo đạc: Bàn đo lún sâu, lún bề mặt, piezometer khí nén & điện tử.      │
  └───────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

  • Chế tạo thiết bị thực nghiệm: Thiết kế hộp nén cố kết chuyên dụng có gắn đầu đo áp lực rỗng điện tử tại đáy dao vòng, bổ sung vòng đệm cao su kín (O-ring) chống rò rỉ áp lực tại khớp nối. Dao vòng chế tạo bằng inox chống ăn mòn với 3 chiều cao khác nhau ($H = 2$ cm, $4$ cm, $6$ cm, đường kính $D = 60$ mm) nhằm thay đổi chiều dài đường thấm $H_{dr}$.
  • Giao thức thu thập dữ liệu: Tự động hóa hoàn toàn bằng bộ xử lý tín hiệu đa kênh kết nối máy tính, ghi nhận đồng thời biến dạng lún ($\Delta H$) với độ chính xác $0.001$ mm và biến thiên ALNLR ($\Delta u$) với độ nhạy $0.1$ kPa ở các bước thời gian logarit: $0.1, 0.25, 0.5, 1, 2, 4, 8, 15, 30, 60, 120, \dots, 1440$ phút (mở rộng tới 72 giờ cho nén một chiều).
  • Quy trình kiểm chuẩn và Tam giác đạc (Triangulation):
    • Data triangulation: So sánh dữ liệu lún phòng thí nghiệm với thí nghiệm hiện trường xuyên tĩnh đo áp lực nước rỗng (CPTu) tại dự án Tân Sơn Nhất – Bình Lợi.
    • Method triangulation: Kiểm định hệ số cố kết ngang $C_h$ và đứng $C_v$ qua 3 phương pháp độc lập: Casagrande ($\log t$), Taylor ($\sqrt{t}$), và Asaoka quan trắc tiêu tán.

Data và phân tích

Mẫu thí nghiệm đại diện cho trầm tích Holocene sét mềm châu thổ Mekong:

  • Đặc tính cơ lý nền sét Quận 7 và Vĩnh Long: Độ ẩm tự nhiên $w = 70 \div 130%$; Khối lượng thể tích $\gamma_w = 1.3 \div 1.8$ g/cm³; Sức chống cắt không thoát nước cánh cắt hiện trường (VST) $S_u = 8 \div 15$ kPa; Chỉ số nén $C_c = 0.6 \div 1.2$; Hệ số rỗng ban đầu $e_0 = 1.8 \div 3.2$; $OCR = 1.0 \div 1.15$ (cố kết bình thường đến tiền cố kết nhẹ); Tỷ số $C_h/C_v$ đo từ CPTu dao động trong khoảng $1.6 \div 1.8$.
  • Kỹ thuật tính toán: Lập trình thuật toán phân tích tổng lớp phân tố trong bài toán phẳng 2D. Nền đất được chia thành mạng lưới vi phân các lớp mỏng $\Delta z = 0.5 \div 1.0$ m và dải $\Delta x = 1.0$ m. Trường ứng suất $\sigma_z(x,z), \sigma_x(x,z), \tau_{xz}(x,z)$ được tính toán giải tích, sau đó xác định $\sigma_m(x,z)$. Quá trình gia tải nhiều giai đoạn (staged construction) được mô phỏng bằng cách cộng dồn gia số ứng suất và giải bài toán tiêu tán ALNLR theo từng bước thời gian $\Delta t$.

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

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

  1. Quy luật $u_{max} \approx \sigma_m$: Thí nghiệm trên các mẫu sét M1 đến M7 ở các cấp áp lực $50 \to 100$ kPa, $100 \to 200$ kPa và $400 \to 800$ kPa khẳng định ALNLR thặng dư lớn nhất không bao giờ đạt tới giá trị gia số ứng suất đứng $\Delta p$. Giá trị $u_{max}$ đo được chỉ bằng $55% \div 75%$ giá trị $\Delta p$, tương đương hoàn toàn với ứng suất trung bình $\sigma_m = \frac{\Delta \sigma_z + 2\Delta \sigma_x}{3}$ (với hệ số áp lực hông $K_0 \approx 0.5 \div 0.6$).
  2. Hiện tượng trễ pha đạt đỉnh ALNLR (50–70 phút): Đối với mọi chiều cao mẫu $H = 2$ cm, $4$ cm, $6$ cm, ALNLR thặng dư không tăng vọt tức thời mà tăng dần đều, đạt giá trị đỉnh $u_{max}$ sau thời gian $t = 50 \div 70$ phút sau khi chất tải. Điều này chứng minh quá trình truyền tải từ chất lưu sang cốt hạt bị chi phối bởi độ nhớt cấu trúc và sự hòa tan của bọt khí, hoàn toàn độc lập với chiều dài đường thấm $H_{dr}$.
  3. Tỷ lệ phân chia biến dạng thể tích – biến dạng trượt: Ứng dụng mô hình 2D tính toán cho nền đường đắp Quốc lộ 1A (Vĩnh Long) cho thấy biến dạng thể tích do cố kết ($S_v$) chiếm $80% \div 95%$ tổng độ lún; biến dạng trượt do chuyển vị ngang ($S_s$) chiếm $5% \div 20%$. Chuyển vị ngang tập trung chủ yếu ở độ sâu từ $0 \div 6$ m dưới chân taluy, gây lún bổ sung từ $15 \div 30$ cm.
  4. Khả năng dự báo lún lệch chính xác: Phương pháp đề nghị triệt tiêu hiện tượng ước lượng độ lún quá mức ở tâm tải trọng ($S_{center}$) và dưới mức ở mép taluy ($S_{edge}$) của phương pháp 1D. Độ lún lệch $\Delta S = S_{center} - S_{edge}$ theo thời gian phản ánh chính xác quy luật biến dạng hình học thực tế.
  5. Sự tương thích vượt trội tại công trình gia cố giếng cát (Cầu Cao Lãnh): Mô hình phần tử đơn vị phẳng 2D mô phỏng chính xác quá trình thi công đắp 11 giai đoạn (từ cao độ $+2.0$ m đến $+7.5$ m). Độ lún tính toán tại Km 8+100 đạt $182$ cm sau 450 ngày, sai lệch dưới $4.5%$ so với độ lún quan trắc bàn đo ($188$ cm), vượt trội hoàn toàn so với mô hình 1D (sai số $>22%$).
                    SO SÁNH TIẾN TRÌNH TIÊU TÁN ALNLR VÀ ĐỘ LÚN
 
      ALNLR (u)                              Độ lún (S)
        ▲                                      ▲
   Δσ_z ┤- - - - (Terzaghi 1D: u0 = Δσz)       │              (Mô hình 2D đề nghị)
        │                                      │           ┌──────────────────────
    σ_m ┼────┐ (Đỉnh trễ 50-70 phút)           │         ┌─┘
        │    │╲                                │       ┌─┘    (Terzaghi 1D:
        │   ┌┘  ╲                              │     ┌─┘       Lún nhanh & tắt sớm)
        │  ┌┘    ╲  (Mô hình 2D đề tài)        │    ┌┘- - - - - - - - - - - - - - -
        │ ┌┘      ╲                            │   ┌┘
        │┌┘        ╲                           │ ┌─┘
      0 └┴──────────┴──────────────► Thời gian 0 └┴─────────────────────────► Thời gian
        0 50-70m      t_cố_kết                 0        t_ổn_định

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung vào cơ học đất không bão hòa và bão hòa luận điểm liên kết giữa trường ứng suất tensor bậc hai với phương trình truyền chất lưu lỗ rỗng, tạo bước chuyển dịch từ bài toán 1D đơn giản hóa sang bài toán phẳng 2D có cơ sở vật lý vững chắc.
  • Về mặt phương pháp luận: Cung cấp quy trình thí nghiệm nén cố kết có đo ALNLR chuẩn hóa, cho phép các phòng thí nghiệm LAS-XD xác định đồng thời hệ số thấm phi tuyến $k(e)$, hệ số cố kết $C_v, C_h$, và hệ số áp lực rỗng $B$.
  • Về ứng dụng thực tiễn và chính sách:
    • Cho phép tư vấn thiết kế tính toán chính xác khối lượng đất đắp bù lún (vốn thường xuyên gây tranh chấp thanh quyết toán giữa Chủ đầu tư và Nhà thầu do lún vượt dự toán $15 \div 30%$).
    • Tối ưu hóa thời gian gia tải chờ lún cho các dự án cao tốc Bắc – Nam phía Đông, đoạn qua ĐBSCL (Cần Thơ – Cà Mau, Châu Đốc – Cần Thơ – Sóc Trăng).

Limitations và Future Research

  1. Giới hạn về đặc tính tải trọng: Mô hình hiện tại mới xét đến tác dụng của tải trọng tĩnh phân bố hình thang (tải trọng bản thân nền đắp), chưa tích hợp thành phần tải trọng động trùng phục (cyclic traffic loading) do các đoàn xe tải nặng lưu thông trên mặt đường.
  2. Giới hạn về giả định môi trường: Giả định vật liệu đất sét đẳng hướng từng lớp về mặt đàn hồi ($E, \nu$), chưa xét đến tính bất đẳng hướng cấu trúc tự nhiên phức tạp (anisotropy) và hiệu ứng phá hoại dẻo phi tuyến bậc cao khi nền đạt trạng thái giới hạn dẻo (plastic yield state).
  3. Định hướng nghiên cứu 5–10 năm tới:
    • Phát triển mô hình cố kết 2D/3D liên kết động – nhiệt – thủy – cơ (THMC Coupled Modeling) cho nền đất yếu ven biển chịu biến đổi khí hậu và triều cường.
    • Ứng dụng mạng nơ-ron nhân tạo (PINNs - Physics-Informed Neural Networks) để giải ngược trường hệ số thấm $k_h(x,z)$ từ chuỗi số liệu quan trắc piezometer thời gian thực.
    • Mở rộng kiểm chứng cho các công nghệ xử lý đất mới như cọc đất trộn xi măng (CDM), cọc cát đầm chặt (SCP) kết hợp vải địa kỹ thuật cường độ cao.

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

  • Học thuật & Trích dẫn: Đặt nền móng lý thuyết mới cho các đề tài nghiên cứu sinh và học viên cao học chuyên ngành Địa kỹ thuật – Công trình ngầm. Tiềm năng trích dẫn cao trong các tạp chí ISI/Scopus thuộc nhóm Q1/Q2 (như Geotechnique, Computers and Geotechnics, Canadian Geotechnical Journal).
  • Chuyển đổi ngành xây dựng giao thông: Thay đổi căn bản quy trình tính toán thiết kế nền đắp trên đất yếu trong các tiêu chuẩn ngành (hướng tới cập nhật thay thế hoặc bổ sung cho TCVN 9355:2012 và TCVN 4200:2012).
  • Hiệu quả kinh tế – xã hội: Giảm thiểu rủi ro trượt lở mái taluy đường đầu cầu, giảm chi phí duy tu sửa chữa hư hỏng lượn sóng mặt đường sau khi đưa vào khai thác, tiết kiệm hàng trăm tỷ đồng cho ngân sách nhà nước từ việc chuẩn hóa khối lượng đắp bù lún.

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

  • Nghiên cứu sinh & Giảng viên Địa kỹ thuật: Tiếp cận mô hình giải tích 2D hoàn chỉnh, hệ phương trình vi phân nhiều pha và phương pháp chế tạo buồng nén cải tiến để mở rộng nghiên cứu học thuật.
  • Kỹ sư Thiết kế & Tư vấn Giám sát: Nắm bắt công cụ tính toán phân lớp phân tố chính xác để dự báo lún tâm, lún mép và tốc độ tiêu tán ALNLR cho các công trình đắp thực tế.
  • Chủ đầu tư & Ban Quản lý Dự án Giao thông: Có căn cứ pháp lý và kỹ thuật minh bạch để nghiệm thu khối lượng đất đắp bù lún và chuẩn hóa tiến độ dỡ tải trước khi thi công móng mặt đường.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc chứng minh và toán học hóa quy luật: Áp lực nước lỗ rỗng thặng dư cực đại sinh ra trong nền sét mềm bão hòa nước dưới công trình đắp phụ thuộc vào ứng suất trung bình $\sigma_m = (\sigma_z + 2\sigma_x)/3$ thay vì gia số ứng suất đứng $\Delta \sigma_z$. Luận án đã mở rộng lý thuyết cố kết đàn hồi nhiều pha của Biot (1941) và hiệu chỉnh lý thuyết cố kết một chiều Terzaghi (1925), giải quyết được mâu thuẫn tồn tại hàng thập kỷ giữa lý thuyết và kết quả đo đạc piezometer thực địa.

2. Điểm mới về phương pháp luận nghiên cứu khi so sánh với các công bố quốc tế trước đây?

So với các nghiên cứu quy đổi bấc thấm của Hird (1995) hay Indraratna & Redana (1997) vốn chủ yếu biến đổi hệ số thấm tương đương ($k_{hp}$) trên nền bài toán đàn hồi tuyến tính cố định, luận án đã:

  • Thiết kế thành công buồng nén cố kết có sensor đo áp lực rỗng đáy mẫu ở 3 cấp chiều cao ($2, 4, 6$ cm), bắt được thời gian trễ pha đạt đỉnh 50–70 phút.
  • Xây dựng giải thuật 2D biến dạng phẳng tích hợp phân lớp phân tố, cho phép tính toán lún dựa trên tách rời biến dạng cắt ($S_s$) và biến dạng thể tích ($S_v$) theo tiến trình đắp từng giai đoạn thực tế.

3. Phát hiện thực nghiệm nào bất ngờ nhất và bằng chứng dữ liệu đi kèm?

Phát hiện bất ngờ nhất là thời gian đạt ALNLR cực đại ($t \approx 50 \div 70$ phút) hoàn toàn không đổi khi thay đổi chiều cao mẫu đất từ $2$ cm lên $4$ cm và $6$ cm. Theo lý thuyết thấm thuần túy, thời gian truyền áp lực phải tỷ lệ thuận với bình phương chiều dài đường thấm ($H_{dr}^2$). Bằng chứng thực nghiệm này chứng minh hiện tượng tích tụ áp lực nước lỗ rỗng ban đầu bị khống chế bởi động học co nén của bọt khí hòa tan và tính nhớt dẻo của khung cốt hạt sét, chứ không đơn thuần là quá trình thủy lực.

4. Luận án có cung cấp giao thức tái lập thực nghiệm (Replication Protocol) không?

Có. Bản thuyết minh luận án trình bày chi tiết bản vẽ cấu tạo hộp nén cố kết cải tiến (vòng đệm cao su kín O-ring, đầu kết nối transducer áp lực đáy), thông số kỹ thuật bộ thu nhận dữ liệu, các bước gia tải đẳng hướng và bất đẳng hướng, cùng toàn bộ thuật toán lập trình giải tích phân lớp phân tố cho bài toán phẳng 2D.

5. Lộ trình nghiên cứu phát triển 10 năm tiếp theo được vạch ra ra sao?

Lộ trình 10 năm bao gồm: Tích hợp hiệu ứng tải trọng động giao thông vào mô hình cố kết 2D; số hóa thuật toán thành module chuyên dụng tích hợp vào các phần mềm mã nguồn mở; mở rộng ứng dụng cho bài toán cố kết kết hợp cọc đất xi măng sâu và các giải pháp gia cố nền bền vững thích ứng với sụt lún quy mô vùng tại ĐBSCL.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán xác định độ lún và lún lệch theo thời gian cho nền đất yếu có bề dày lớn dưới công trình đắp dạng tuyến theo sơ đồ hai chiều biến dạng phẳng.
  2. Xác lập bằng chứng thực nghiệm vững chắc khẳng định ALNLR thặng dư cực đại $u_{max}$ xấp xỉ ứng suất trung bình $\sigma_m$, có độ trễ đạt đỉnh 50–70 phút do tương tác đa pha cốt đất – nước – bọt khí.
  3. Đề xuất thành công phương pháp phân chia độ lún thành 2 thành phần độc lập: Biến dạng thể tích $S_v$ (chiếm 80–95%) và biến dạng trượt do chuyển vị ngang $S_s$ (chiếm 5–20%), triệt tiêu sai số của sơ đồ 1D cổ điển.
  4. Xây dựng mô hình phần tử đơn vị phẳng cho nền gia cố giếng cát/bấc thấm kết hợp gia tải trước nhiều giai đoạn, đạt độ khớp quan trắc hiện trường tại Quốc lộ 1A và Cầu Cao Lãnh với sai số dưới $5%$.
  5. Mở ra 3 nhánh nghiên cứu học thuật mới: Cố kết tương tác động lực học giao thông, địa kỹ thuật tính toán với vật liệu bất đẳng hướng, và mô hình hóa lan truyền chất lưu đa pha trong đất yếu bão hòa.
  6. Kết quả nghiên cứu tạo nền tảng khoa học trực tiếp phục vụ thiết kế, kiểm soát ổn định, và quản lý chi phí đắp bù lún cho các đại dự án hạ tầng giao thông trọng điểm tại khu vực phía Nam.