Tổng quan về luận án

Bối cảnh phát triển cơ sở hạ tầng giao thông, khu công nghiệp và hệ thống đê điều tại Việt Nam, đặc biệt là vùng châu thổ sông Hồng và sông Mê Kông, đang đối mặt với thách thức nghiêm trọng từ các tầng trầm tích đất yếu có chiều dày lớn. Đất yếu tại các khu vực này mang đặc trưng cơ lý bất lợi: hệ số rỗng lớn ($e > 1.5$), mô đun biến dạng thấp ($E < 5000\text{ kN/m}^2$), sức kháng cắt không thoát nước rất nhỏ ($S_u \le 15\text{ kN/m}^2$), góc ma sát trong $\varphi_u \approx 0$, cùng tính nén lún và độ ẩm vượt quá giới hạn chảy ($B > 1$). Khi chịu tải trọng công trình lớn hoặc tải trọng động từ giao thông, nền đất thường xuất hiện biến dạng lún dư kéo dài và nguy cơ trượt trồi. Phương pháp truyền thống sử dụng cọc hạt rời (cọc cát, cọc đá không bọc) gặp giới hạn tới hạn: khi $S_u < 15\text{ kN/m}^2$, đất xung quanh không cung cấp đủ áp lực nở hông để giữ ổn định thân cọc, dẫn đến hiện tượng phình ngang và phá hủy cắt cục bộ ở phần đầu cọc.

Luận án tiến sĩ kỹ thuật địa chất của NCS Lê Quân với đề tài "Nghiên cứu ứng dụng cọc hạt rời bọc vải địa kỹ thuật trong cải tạo nền đất yếu" (2023) đặt nền móng tiên phong trong việc thiết lập hệ thống tính toán, mô phỏng số và kiểm chứng thực nghiệm hiện trường toàn diện cho công nghệ cọc hạt rời bọc vải địa kỹ thuật (Geosynthetic Encased Columns - GEC) tại Việt Nam.

Khoảng trống nghiên cứu (Research Gap) được xác định rõ ràng: Dù lý thuyết giải tích lăng trụ của Raithel và Kempfert (2000) cùng các nghiên cứu quốc tế (Ghionna & Jamiolkowski, 1981; Almeida et al., 2013; Castro & Sagaseta, 2011) đã chứng minh hiệu quả gia cường của áo bọc địa kỹ thuật, nhưng tại Việt Nam hoàn toàn chưa có công trình nào đánh giá có hệ thống sự tương tác phức hợp giữa cọc hạt rời bọc vải và đất bùn sét bão hòa nước thông qua mô hình số phi tuyến nâng cao kết hợp quan trắc thực địa.

Các câu hỏi nghiên cứu (Research Questions) và giả thuyết khoa học (Hypotheses) được xác lập:

  • RQ1: Cơ chế biến dạng ngang và tương tác truyền ứng suất giữa cọc hạt rời, vỏ bọc địa kỹ thuật và nền đất yếu diễn ra như thế nào dưới tác dụng của tải trọng ngoài?
  • RQ2: Mức độ chính xác và phạm vi ứng dụng của mô hình đàn dẻo Mohr-Coulomb so với mô hình đất yếu Soft Soil khi mô phỏng số bài toán cọc GEC trong phần mềm PLAXIS 2D?
  • RQ3: Mức độ sai lệch về dự báo lún và biến dạng ngang giữa phương pháp giải tích (Analytical Method - AM) và phương pháp phần tử hữu hạn (Finite Element Method - FEM) là bao nhiêu?
  • RQ4: Mức độ ảnh hưởng của các thông số hình học (đường kính $d_c$, khoảng cách $S$, chiều dày lớp đất yếu $L$) và vật liệu (độ cứng chịu kéo của vải $J$, mô đun biến dạng cọc $E_c$) đến hiệu quả giảm lún của hệ thống?

Hệ thống giả thuyết tương ứng:

  • H1: Lực kéo chu vi (Hoop tension) trong vải địa kỹ thuật cung cấp áp lực hông bổ sung $\sigma_{r,g}$, khống chế biến dạng phình ngang của cọc ở độ sâu $1\text{D} \div 3\text{D}$, giúp tăng khả năng chịu tải và giảm lún từ $30% \div 40%$ so với cọc hạt rời thông thường.
  • H2: Mô hình cấu hủy Soft Soil mô tả chính xác hơn mô hình Mohr-Coulomb đối với ứng xử nén lún phi tuyến của đất bùn sét yếu khi so sánh với quan trắc nén tĩnh hiện trường và phương pháp phân tích ngược Asaoka.
  • H3: Phương pháp giải tích theo Raithel & Kempfert (2000) cho kết quả dự báo biến dạng lún và biến dạng ngang lớn hơn có hệ thống so với phương pháp phần tử hữu hạn do các giả thiết bảo thủ về trạng thái ứng suất giới hạn chủ động.
  • H4: Tăng cường độ chịu kéo $J$ của vải địa kỹ thuật và giảm khoảng cách $S$ giữa các cọc tạo ra hiệu ứng tăng tỷ số tập trung ứng suất $n$, kiểm soát triệt để độ lún cố kết của nền gia cố.

Khung lý thuyết (Theoretical Framework) tích hợp: Lý thuyết cân bằng lăng trụ đơn vị (Unit Cell Model), lý thuyết kéo màng mỏng đàn hồi tuyến tính (Hoop Stress Theory), mô hình cơ học đất đàn dẻo biến dạng dẻo cứng (Hardening/Soft Soil Plasticity Framework) và lý thuyết cố kết thấm Biot.

Phạm vi và quy mô nghiên cứu: Dữ liệu thực nghiệm, địa chất công trình và quan trắc lún thu thập từ dự án nhà máy VIFON II (tỉnh Long An), kết hợp khảo sát tham số trên 4 dải đường kính cọc ($d_c = 0.65; 0.8; 1.0; 1.2\text{ m}$), 4 khoảng cách cọc ($S = 1.3; 1.95; 2.6; 3.25\text{ m}$), 4 chiều dày đất yếu ($L = 3.5; 7.0; 14.0; 17.5\text{ m}$), 3 cấp mô đun cọc đá ($E_c = 48000; 80000; 120000\text{ kN/m}^2$) và 5 cấp cường độ chịu kéo của vải địa kỹ thuật ($J = 1000 \div 5000\text{ kN/m}$) dưới 6 cấp áp lực gia tải ($P = 41 \div 246\text{ kN/m}^2$).

Literature Review và Positioning

Lịch sử phát triển công nghệ cọc hạt rời bọc vải địa kỹ thuật (GEGC/GEC) bắt đầu từ các ý tưởng ban đầu của Van Impe và Silence, sau đó được Ghionna và Jamiolkowski (1981) công thức hóa thành lời giải giải tích dựa trên lý thuyết đàn hồi để xác định yêu cầu cường độ vải bao bọc ở trạng thái giới hạn. Tuy nhiên, mô hình của Ghionna & Jamiolkowski chưa xét đến tương tác biến dạng tương hỗ giữa vải và cọc. Bước đột phá lý thuyết diễn ra vào năm 2000 khi Raithel và Kempfert công bố phương pháp giải tích hoàn chỉnh dựa trên mô hình lăng trụ đối xứng trục (Unit Cell), thiết lập sự cân bằng biến dạng đứng giữa cọc và đất yếu xung quanh, đồng thời tích hợp ứng suất kéo chu vi của vải địa kỹ thuật $\sigma_{r,g} = J \cdot \Delta r_g / r_g^2$.

Tiến trình phát triển lý thuyết Cọc Hạt Rời Bọc Vải Địa Kỹ Thuật (1981 - 2023):

[1981] Ghionna & Jamiolkowski: Lời giải đàn hồi trạng thái giới hạn
       │
[2000] Raithel & Kempfert: Mô hình lăng trụ (Unit Cell Model) cân bằng biến dạng
       │
[2009-2011] Wu & Hong / Pulko / Castro & Sagaseta: Mô hình hóa phi tuyến & độ nhớt dẻo
       │
[2013-2017] Almeida / Hongyang Cheng: Ứng xử 3 trục, quan hệ ứng suất - độ giãn nở thể tích
       │
[2021] Yang Zhou et al.: Hiệu ứng vòm đất đắp và tính phi tuyến cọc hạt rời
       │
[2023] Luận án NCS Lê Quân: Tích hợp AM - FEM (Soft Soil) - Thí nghiệm bàn nén VIFON II

Các dòng nghiên cứu tiếp theo chứng kiến những tranh luận học thuật sâu sắc:

  1. Tranh luận về phân bố biến dạng ngang: Raithel và Kempfert (2000) giả định thể tích cọc không đổi và biến dạng ngang phân bố đồng đều dọc theo chiều dài cọc. Trái lại, các nghiên cứu thực nghiệm và mô phỏng của Lee et al. (2007), Khabbazian et al. (2009), Murugesan và Rajagopal (2010) chỉ ra rằng tải trọng đỉnh cọc gây ra biến dạng nén dọc trục đi kèm hiện tượng phình ngang bất đối xứng, tập trung cục bộ tại độ sâu từ $1\text{D} \div 3\text{D}$ tính từ đỉnh cọc, giảm nhanh về không ở các tầng sâu hơn.
  2. Tranh luận về ứng xử cơ học của cọc và đất: Một trường phái (Zhang & Zhao, 2015; Pulko et al., 2011) sử dụng mô hình đàn - dẻo lý tưởng Mohr-Coulomb, trong khi nhóm nghiên cứu của Hongyang Cheng et al. (2017) và Yang Zhou et al. (2021) nhấn mạnh rằng cọc hạt rời thể hiện rõ tính phi tuyến, sự thay đổi thể tích (dilatancy) và tỷ số ứng suất phụ thuộc biến dạng cắt. Khi bỏ qua ứng xử phi tuyến này, tỷ số tập trung ứng suất $n$ giữa cọc và đất sẽ bị dự báo vượt mức thực tế.

So sánh với hai dự án và nghiên cứu quốc tế điển hình:

  • Dự án mở rộng nhà máy Airbus tại Hamburg (Đức, 2002): Công trình cải tạo đê bao $2500\text{ m}$ trên tầng bùn sét dày $4 \div 14\text{ m}$ ($S_u = 0.4 \div 10\text{ kN/m}^2$) bằng $60,000$ cọc GEC ($d = 0.8\text{ m}$, $J = 1700 \div 2800\text{ kN/m}$) đã giúp "tiết kiệm 35,000 tấn thép do không cần tường cừ, giảm 1.1 tỷ $\text{m}^3$ cát san lấp do không phải tạo mái dốc và rút ngắn thời gian thi công từ 3 năm xuống 8 tháng" [trích dẫn nguyên văn dữ liệu từ Luận án].
  • Nghiên cứu nén ba trục kích thước lớn của Jianfeng Xue et al. (2018): Thí nghiệm trên 16 cọc đá bọc vải địa kỹ thuật ($J = 33, 43, 52, 65\text{ kN/m}$) dưới các áp lực buồng $50, 100, 150, 200\text{ kPa}$ đã chứng minh lớp bọc vải đóng vai trò như một lực dính biểu kiến hoặc áp lực buồng bổ sung, làm thay đổi căn bản đường cong quan hệ ứng suất - biến dạng từ dạng hyperbolic sang tuyến tính tăng bền.

Tại Việt Nam, các nghiên cứu trước đây của Lê Bá Vinh và Lê Bá Khánh (2017), Lương Phương Hợp (2017), Đỗ Mạnh Hùng (2015), Lương Văn Hoàng (2021) chủ yếu tập trung vào cọc cát, cọc đá đầm rung không bọc và chỉ áp dụng mô hình Mohr-Coulomb đơn giản. Luận án của Lê Quân định vị chính xác khoảng trống học thuật này bằng cách thiết lập nghiên cứu có hệ thống đầu tiên tại Việt Nam, kết nối phương pháp giải tích Raithel & Kempfert với mô phỏng số PLAXIS v2018 sử dụng mô hình Soft Soil, được kiểm chứng trực tiếp bằng chuỗi dữ liệu nén tĩnh cọc đơn và quan trắc lún bàn nén tại công trình nhà máy VIFON II.

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

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

Luận án đóng góp vào cơ học đất và kỹ thuật nền móng các luận điểm nền tảng:

  • Mở rộng lý thuyết giam giữ hông (Confinement Theory): Làm sáng tỏ bản chất của áo bọc địa kỹ thuật trong việc chuyển đổi ứng suất nén dọc trục thành ứng suất kéo vòng (Hoop stress), tạo ra ứng suất nén hông bổ sung $\sigma_{r,g}$. Ứng suất này bù đắp trực tiếp sự thiếu hụt áp lực ngang của đất yếu bão hòa nước có $S_u < 15\text{ kN/m}^2$.
  • Mô hình hóa tương tác phi tuyến cọc - vải - đất: Khắc phục hạn chế của lý thuyết đàn hồi cổ điển thông qua việc chứng minh rằng cơ chế truyền lực thực tế nằm ở trạng thái trung gian giữa cân bằng ứng suất ($\Delta \sigma_c = \Delta \sigma_s$) và cân bằng biến dạng ($S_{cl} = S_{sl}$).
  • Hệ thống mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề P1: Dưới tác dụng của tải trọng thẳng đứng $\Delta \sigma_z$, biến dạng phình hông của cọc hạt rời kích hoạt sức căng $T_g$ trong vải địa kỹ thuật, làm tăng góc ma sát trong hiệu dụng và lực dính biểu kiến của khối hạt rời.
    • Mệnh đề P2: Vùng biến dạng trượt dẻo cực đại không phân bố đều mà định vị tập trung trong khoảng $1.0\text{D} \div 3.0\text{D}$ từ đỉnh cọc, trùng khớp với quy luật mặt trượt cắt Rowe & Barden (1964).
    • Mệnh đề P3: Tỷ số tập trung ứng suất $n = \Delta \sigma_c / \Delta \sigma_s$ tỷ lệ thuận với mô đun độ cứng của vải bọc $J$ và độ cứng cọc $E_c$, nhưng tỷ lệ nghịch với diện tích lăng trụ thay thế $A_e$.
    • Mệnh đề P4: Mô hình đàn dẻo phi tuyến Soft Soil phản ánh chính xác quy luật suy giảm thể tích và gia tăng áp lực nước lỗ rỗng so với mô hình Mohr-Coulomb trong bài toán tải trọng nền đắp.
Khung Phân Tích Cơ Chế Truyền Lực & Tương Tác Ứng Suất:

   [ Ngoại lực thẳng đứng Δσ_z ]
                 │
      ┌──────────┴──────────┐
      ▼                     ▼
[ Cọc Hạt Rời: Δσ_c ]    [ Đất Yếu Xung Quanh: Δσ_s ]
      │                     │
   (Nén dọc trục)        (Chịu ứng suất ngang σ_r,s)
      ▼                     │
[ Biến dạng ngang Δr_c ]    │
      │                     │
      ▼                     │
[ Sức căng Hoop T_g ]       │
      │                     │
      ▼                     │
[ Áp lực hông bổ sung σ_r,g ├───────────────────────┘
      │                     │
      └──────────┬──────────┘
                 ▼
[ Tái lập Cân bằng Ứng suất Ngang: Δσ_r = σ_r,c - σ_r,s - σ_r,g = 0 ]
                 ▼
[ Kiểm soát Biến dạng Lún S' và Phình ngang Tối đa ]

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

Khung phân tích của luận án tích hợp đồng thời 3 trường phái lý thuyết:

  1. Lý thuyết mô hình lăng trụ (Unit Cell Approach) theo Raithel & Kempfert (2000): Thiết lập phương trình cân bằng ngoại lực: $$\Delta \sigma_s A_s + \Delta \sigma_c A_c = \Delta \sigma_z A_e$$ với tỷ số diện tích thay thế $a_s = A_c / A_e = C (d_c / s)^2$.
  2. Lý thuyết đàn dẻo Soft Soil (Cam-Clay cải biên): Mô tả đất yếu qua chỉ số nén biến đổi theo áp lực tiền cố kết, phân tách biến dạng đàn hồi và dẻo.
  3. Lý thuyết sức căng màng mỏng (Hoop Tensile Mechanics): $$\sigma_{r,g} = \frac{T_g}{r_g} = J \frac{\Delta r_c - (r_g - r_c)}{r_g^2}$$

Điều kiện biên xác định (Boundary Conditions): Bài toán được giải trong không gian đối xứng trục 2D; điều kiện biên cơ học triệt tiêu biến dạng ngang tại biên ngoài lăng trụ ($\varepsilon_r = 0$), biên đáy cố kết cứng hoàn toàn ($u_x = u_y = 0$); điều kiện biên thủy lực cho phép thoát nước tự do theo phương đứng và hướng tâm vào lõi cọc cát/đá.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết học thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), sử dụng phương pháp hỗn hợp (Mixed-Methods) kiểm chứng đa giác quan giữa mô hình toán học giải tích, mô phỏng số học phần tử hữu hạn và quan trắc thực địa.

Thiết kế nghiên cứu đa cấp độ (Multi-Level Design):

  • Cấp độ phần tử (Element Level): Mô phỏng hành vi ứng xử đàn dẻo của vật liệu hạt rời và bùn sét yếu dưới các đường dẫn ứng suất khác nhau.
  • Cấp độ đơn cọc / Lăng trụ đơn vị (Unit Cell Level): Phân tích tương tác cọc - vải - đất trong khối lăng trụ đối xứng trục.
  • Cấp độ toàn thể (Global Foundation Level): Đánh giá độ lún tổng thể của nền móng công trình nhà máy VIFON II dưới tác dụng của tải trọng sàn và hoạt tải sử dụng.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu tuân thủ các tiêu chuẩn kỹ thuật nghiêm ngặt:

  • Khảo sát địa chất hiện trường: Dữ liệu địa chất thu thập từ công trình nhà máy VIFON II (Long An) bao gồm hố khoan địa chất sâu, thí nghiệm xuyên tĩnh CPT đo áp lực nước lỗ rỗng ($u$), xác định tầng bùn sét yếu bão hòa nước dày từ $3.5\text{ m}$ đến $17.5\text{ m}$.
  • Thí nghiệm nén tĩnh cọc hiện trường (Static Plate Load Test): Tiến hành trên 3 cọc thử nghiệm điển hình (S1-1-317, S5-2-266, S8-1-202) với quy trình gia tải từng cấp, đo chuyển vị lún đầu cọc bằng đồng hồ so có độ chính xác $0.01\text{ mm}$ theo chu kỳ gia tải và giảm tải tiêu chuẩn.
  • Quan trắc lún dài hạn: Đo đạc độ lún bề mặt nền đất tại 4 phân khu quan trắc (S7, S8, S9, S10) trong suốt quá trình đắp gia tải và chờ cố kết.
  • Tam giác hóa phương pháp (Methodological Triangulation): Kết quả độ lún thực tế được xử lý bằng phương pháp quan trắc ngược Asaoka, đối sánh trực tiếp với kết quả chạy mô hình số PLAXIS và công thức giải tích AM.
Sơ đồ Tam Giác Hóa Dữ Liệu & Phương Pháp Nghiên Cứu:

                [ THỰC NGHIỆM HIỆN TRƯỜNG ]
                 (Công trình VIFON II, Long An)
                 ┌──────────────┴──────────────┐
                 ▼                             ▼
        [ Thí nghiệm Bàn nén ]       [ Quan trắc lún nền ]
        (Cọc S1, S5, S8)             (Phân khu S7, S8, S9, S10)
                 │                             │
                 │                             ▼
                 │                   [ Phân tích ngược Asaoka ]
                 │                             │
                 └──────────────┬──────────────┘
                                ▼
                   [ KIỂM CHỨNG & ĐỐI SÁNH ]
                                ▲
                 ┌──────────────┴──────────────┐
                 ▼                             ▼
        [ MÔ PHỎNG SỐ FEM ]          [ PHƯƠNG PHÁP GIẢI TÍCH ]
        (PLAXIS 2D v2018)            (Mô hình Raithel & Kempfert)
        - Mohr-Coulomb               - Cân bằng lăng trụ Unit Cell
        - Soft Soil Model            - Ứng suất kéo Hoop

Data và phân tích

Mô phỏng số được thực hiện trên phần mềm bản quyền PLAXIS 2D v2018. Lưới phần tử hữu hạn sử dụng phần tử tam giác 15 nút bậc cao (15-node triangular elements) với 12 điểm lấy tích phân Gauss Gauss, đảm bảo độ chính xác vượt trội trong việc tính toán gradient ứng suất và biến dạng trượt tại vùng tiếp xúc giữa cọc và đất.

Bảng thông số cơ lý đầu vào chuẩn hóa từ báo cáo địa chất VIFON II:

Thông số cơ lý Ký hiệu (Đơn vị) Đất bùn sét yếu Cọc hạt rời (Đá 1x2) Lớp vải địa kỹ thuật
Trọng lượng thể tích tự nhiên $\gamma\text{ (kN/m}^3\text{)}$ $15.5 \div 16.5$ $19.0 \div 20.0$ -
Trọng lượng thể tích bão hòa $\gamma_{sat}\text{ (kN/m}^3\text{)}$ $16.0 \div 17.0$ $20.0 \div 22.0$ -
Mô đun đàn hồi / biến dạng $E\text{ (kN/m}^2\text{)}$ $1500 \div 3500$ $48000 \div 120000$ -
Hệ số Poisson $\nu\text{ (-)}$ $0.35 \div 0.40$ $0.25 \div 0.30$ $0.10$
Lực dính hiệu dụng $c'\text{ (kN/m}^2\text{)}$ $4.0 \div 8.0$ $0.1$ -
Góc ma sát trong $\varphi'\text{ (độ)}$ $5.5^\circ \div 9.0^\circ$ $38.0^\circ \div 45.0^\circ$ -
Chỉ số nén dẻo (Soft Soil) $\lambda^*$ hoặc $C_c$ $0.12 \div 0.18$ - -
Chỉ số nở (Soft Soil) $\kappa^*$ hoặc $C_s$ $0.02 \div 0.04$ - -
Độ cứng chịu kéo của vải $J\text{ (kN/m)}$ - - $1000 \div 5000$

Kiểm tra độ vững chắc (Robustness Checks): Tiến hành phân tích độ nhạy với các dạng chia lưới từ thô (Coarse), trung bình (Medium) đến rất mịn (Very Fine), xác nhận kết quả chuyển vị hội tụ hoàn toàn với lưới mịn quanh thân cọc.

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

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

  1. Hiệu ứng kiểm soát độ lún vượt trội của vải địa kỹ thuật: Kết quả tính toán giải tích và mô phỏng số chỉ ra rằng, việc bọc vải địa kỹ thuật giúp giảm độ lún của nền gia cố từ $30% \div 40%$ so với việc sử dụng cọc hạt rời thông thường không bọc. Khi cường độ chịu kéo $J$ của vải tăng từ $1000\text{ kN/m}$ lên $5000\text{ kN/m}$, độ lún của nền giảm liên tục và đạt trạng thái ổn định tiệm cận.
  2. Quy luật phân bố biến dạng phình ngang: Khác với giả thiết biến dạng đều của các mô hình giải tích sơ khai, kết quả PLAXIS 2D chứng minh biến dạng ngang lớn nhất của cọc GEC luôn tập trung cục bộ tại độ sâu từ $1.0\text{D} \div 3.0\text{D}$ tính từ đỉnh cọc (tương đương độ sâu $1.0 \div 2.5\text{ m}$). Dưới độ sâu $4.0\text{D}$, chuyển vị ngang của cọc tiệm cận về $0$ do áp lực đất xung quanh gia tăng theo độ sâu.
  3. Ưu thế tuyệt đối của mô hình Soft Soil (SSM) so với Mohr-Coulomb (MCM): Khi đối sánh với dữ liệu thí nghiệm nén tĩnh cọc hiện trường tại 3 cọc S1-1-317, S5-2-266 và S8-1-202, mô hình Soft Soil cho kết quả đường cong tải trọng - độ lún bám sát thực nghiệm với độ lệch trung bình dưới $8.5%$, trong khi mô hình Mohr-Coulomb đánh giá thấp độ lún (độ cứng bị phóng đại) với sai số vượt quá $25% \div 35%$.
Vị trí cọc thử tải Độ lún thực nghiệm (mm) Độ lún PLAXIS - Soft Soil (mm) Sai số SSM (%) Độ lún PLAXIS - Mohr-Coulomb (mm) Sai số MCM (%)
S1-1-317 $18.42$ $17.15$ $-6.89%$ $12.30$ $-33.22%$
S5-2-266 $21.15$ $19.60$ $-7.33%$ $14.10$ $-33.33%$
S8-1-202 $16.80$ $15.50$ $-7.74%$ $11.85$ $-29.46%$
  1. Độ lệch có hệ thống giữa phương pháp Giải tích (AM) và Số (FEM): Độ lún và biến dạng ngang tính toán theo phương pháp giải tích Raithel & Kempfert luôn lớn hơn kết quả mô phỏng phần tử hữu hạn từ $12% \div 22%$. Nguyên nhân do phương pháp AM giả thiết cọc làm việc ở trạng thái cân bằng giới hạn chủ động hoàn toàn trên toàn bộ chiều dài, trong khi mô hình số phản ánh chính xác sự phát triển từng bước của các điểm dẻo (plastic points) cục bộ.
  2. Độ nhạy của các tham số hình học và vật liệu:
    • Khoảng cách cọc $S$: Giảm khoảng cách cọc từ $3.25\text{ m}$ ($a_s = 0.05$) xuống $1.3\text{ m}$ ($a_s = 0.31$) làm giảm độ lún tổng thể đến $58%$.
    • Đường kính cọc $d_c$: Tăng $d_c$ từ $0.65\text{ m}$ lên $1.2\text{ m}$ làm tăng tỷ số diện tích thay thế, chuyển đổi cơ chế truyền lực chủ yếu qua thân cọc xuống tầng chịu lực.
    • Mô đun vật liệu cọc $E_c$: Tăng $E_c$ từ $48,000\text{ kN/m}^2$ lên $120,000\text{ kN/m}^2$ giúp giảm lún $18%$, tuy nhiên mức độ đóng góp của $E_c$ giảm dần khi cường độ vải $J$ đạt giá trị cao.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Luận án bổ sung bằng chứng thực nghiệm vững chắc xác nhận tính đúng đắn của lý thuyết vòng kéo Hoop và mô hình tương tác lăng trụ trong điều kiện đất bùn sét bão hòa nước tại Đông Nam Bộ; chứng minh sự cần thiết phải sử dụng mô hình cấu hủy phi tuyến Soft Soil thay thế cho Mohr-Coulomb trong tính toán nền đất yếu.
  • Về phương pháp luận: Thiết lập quy trình chuẩn hóa kết hợp giữa mô hình hóa 2D đối xứng trục trong PLAXIS và phân tích ngược Asaoka từ dữ liệu quan trắc, cung cấp công cụ kiểm định chéo tin cậy cho các kỹ sư địa kỹ thuật.
  • Về thực tiễn kỹ thuật công trình: Cung cấp giải pháp kỹ thuật khả thi cho các công trình đắp cao trên đất cực yếu ($S_u < 15\text{ kN/m}^2$), cho phép tận dụng vật liệu rời địa phương (cát hạt thô, đá dăm nghiền) thay thế các giải pháp cọc bê tông cốt thép đắt tiền, giảm chi phí xây dựng từ $20% \div 30%$.
  • Về chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học định lượng để các cơ quan quản lý chuyên ngành (Bộ Xây dựng, Bộ Giao thông Vận tải) xây dựng và ban hành tiêu chuẩn thiết kế, thi công và nghiệm thu cọc hạt rời bọc vải địa kỹ thuật tại Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ rõ các giới hạn học thuật:

  1. Giả thiết điều kiện thoát nước hoàn toàn: Mô hình giải tích và số trong luận án giả định điều kiện thoát nước (Drained analysis) để triệt tiêu áp lực nước lỗ rỗng thặng dư, chưa mô phỏng chi tiết quá trình tiêu tán áp lực nước lỗ rỗng theo thời gian thực (Coupled consolidation analysis) dưới tải trọng đắp theo từng giai đoạn thi công.
  2. Ứng xử đàn hồi tuyến tính của vải địa kỹ thuật: Vải địa kỹ thuật được giả thiết có ứng xử đàn hồi tuyến tính không suy giảm theo thời gian, chưa xét đến hiện tượng từ biến (creep), chùng ứng suất và tổn thương cơ học cục bộ trong quá trình thi công nén rung.
  3. Bỏ qua hiệu ứng xáo trộn khi thi công (Smear Effect): Quá trình đóng ống vách rung ép tạo cọc gây ra vùng xáo trộn xung quanh thân cọc làm giảm tính thấm và sức kháng cắt của đất yếu lân cận, yếu tố này chưa được đưa vào mô hình số.
  4. Giới hạn bài toán đối xứng trục phẳng: Nghiên cứu tập trung vào mô hình lăng trụ đơn cọc 2D đối xứng trục, chưa phản ánh đầy đủ hiệu ứng tương tác 3D không gian khi cọc chịu tải trọng ngang phức tạp hoặc tải trọng uốn tại góc móng và mái dốc taluy.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Hướng 1: Phát triển mô hình cố kết thấm liên kết (Biot Coupled Consolidation) mô phỏng chính xác sự gia tăng và tiêu tán áp lực nước lỗ rỗng trong quá trình thi công cọc GEC theo thời gian thực.
  • Hướng 2: Nghiên cứu ứng xử từ biến dài hạn (Long-term Creep Behavior) của vải địa kỹ thuật dệt cường độ cao dưới tác dụng của ứng suất kéo chu vi liên tục trong môi trường kiềm/axit của đất phèn miền Tây Nam Bộ.
  • Hướng 3: Mở rộng mô phỏng 3D phần tử hữu hạn đánh giá khả năng chống trượt sâu và ổn định tổng thể của mái dốc đường đầu cầu đắp cao gia cố bằng hệ nhóm cọc GEC.
  • Hướng 4: Nghiên cứu thực nghiệm quy mô lớn tại hiện trường sử dụng các loại vật liệu tái chế (xỉ thép, bê tông nghiền) làm vật liệu hạt rời cho cọc GEC nhằm thúc đẩy kinh tế tuần hoàn.

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

  • Tác động học thuật: Công trình tạo tiền đề học thuật cho các nghiên cứu chuyên sâu về công nghệ cọc bọc địa kỹ thuật tại các trường đại học khối kỹ thuật; dự kiến trở thành tài liệu tham khảo then chốt với khả năng trích dẫn cao trong các nghiên cứu địa kỹ thuật công trình ven biển và đồng bằng.
  • Chuyển đổi công nghiệp xây dựng: Cung cấp giải pháp thay thế kỹ thuật hiệu quả cho các tổng thầu xây dựng hạ tầng, khu công nghiệp, bến bãi logistics tại Đồng bằng sông Cửu Long (Long An, Tiền Giang, Cần Thơ) và ven biển Hải Phòng, Quảng Ninh.
  • Lợi ích xã hội và môi trường: Giảm thiểu nhu cầu sử dụng cọc bê tông cốt thép (hạn chế phát thải khí nhà kính từ sản xuất xi măng và thép), tiết kiệm hàng triệu mét khối cát đắp nhờ tối ưu hóa mái dốc công trình, giảm ô nhiễm tiếng ồn và rung chấn so với giải pháp đóng cọc cừ truyền thống.

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

  • Nghiên cứu sinh và Giảng viên Địa kỹ thuật: Tiếp cận khung lý thuyết hoàn chỉnh, phương pháp giải tích Raithel & Kempfert cải tiến và phương pháp thiết lập mô hình Soft Soil trong PLAXIS 2D.
  • Kỹ sư tư vấn thiết kế và R&D Doanh nghiệp: Sở hữu biểu đồ thiết kế sẵn sàng áp dụng, tra cứu tương quan giữa độ lún, đường kính cọc $d_c$, khoảng cách $S$ và độ cứng vải $J$ cho các tầng đất bùn sét tương đương VIFON II.
  • Nhà thầu thi công hạ tầng: Nắm vững các thông số kiểm soát chất lượng độ chặt cọc hạt rời ($> 75%$), lựa chọn thiết bị thi công nén rung hoặc khoan xoắn phù hợp với điều kiện địa chất đô thị.
  • Cơ quan quản lý nhà nước: Có cơ sở khoa học định lượng để thẩm định an toàn công trình, phê duyệt dự án và xây dựng hệ thống quy chuẩn kỹ thuật quốc gia.

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ì?

Đóng góp độc đáo nhất là việc chứng minh và định lượng cơ chế tương tác trung gian giữa cân bằng ứng suất và cân bằng biến dạng thông qua việc mở rộng lý thuyết cân bằng lăng trụ Raithel & Kempfert (2000). Luận án làm rõ vai trò của ứng suất kéo chu vi màng mỏng $\sigma_{r,g}$ trong việc ngăn chặn phá hoại trượt cắt dẻo của cọc đá trong nền đất cực yếu ($S_u < 15\text{ kN/m}^2$), xác lập mối quan hệ phi tuyến giữa độ cứng chịu kéo $J$ của vải và sự phân bố ứng suất hông bổ sung.

2. Điểm mới về phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu trong nước (Lê Bá Vinh, 2017; Lương Văn Hoàng, 2021) vốn chỉ áp dụng mô hình Mohr-Coulomb tuyến tính cho cọc không bọc, luận án đã:

  • Tích hợp thành công mô hình cấu hủy Soft Soil mô tả chính xác ứng xử nén lún dẻo phi tuyến của bùn sét.
  • Thiết lập quy trình đối sánh tam giác (Triangulation Protocol): Đối chiếu kết quả giữa Phương pháp Giải tích (AM), Phương pháp Phần tử hữu hạn (FEM 15 nút) và Phương pháp Quan trắc hiện trường (Bàn nén tĩnh & Phân tích ngược Asaoka).

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là phương pháp giải tích kinh điển (AM) luôn dự báo độ lún và biến dạng ngang lớn hơn phương pháp số FEM từ $12% \div 22%$. Điều này chứng minh rằng việc giả thiết trạng thái ứng suất chủ động phát triển đồng thời trên toàn bộ chiều dài cọc trong tính toán giải tích là quá bảo thủ, trong khi thực tế cọc chỉ biến dạng cục bộ ở đoạn đầu cọc ($1\text{D} \div 3\text{D}$).

4. Quy trình lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp đầy đủ và chi tiết toàn bộ hồ sơ thiết kế, kích thước hình học lăng trụ, bảng chỉ tiêu cơ lý đất nền và cọc, quy trình gán phần tử Geogrid đàn hồi trong PLAXIS 2D, các bước lặp giải phương trình giải tích một ẩn $\Delta \sigma_s$, cho phép các nhà nghiên cứu độc lập tái lập mô hình với độ chính xác $100%$.

5. Định hướng chương trình nghiên cứu 10 năm tới?

Xây dựng khung mô hình số 3D cố kết thấm hoàn chỉnh xét đến tải trọng xe chạy lặp động (Cyclic traffic loading), nghiên cứu thực nghiệm độ bền suy giảm do mỏi và môi trường hóa chất của vải địa kỹ thuật dệt tại các vùng ngập mặn Đồng bằng sông Cửu Long.

Kết luận

  1. Khẳng định tính ưu việt vượt trội của công nghệ cọc hạt rời bọc vải địa kỹ thuật trong việc gia cố nền đất cực yếu ($S_u < 15\text{ kN/m}^2$), giúp giảm lún tổng thể từ $30% \div 40%$ và nâng cao sức chịu tải của nền đất.
  2. Xác lập tính đúng đắn và sự phù hợp vượt bậc của mô hình cấu hủy Soft Soil so với mô hình Mohr-Coulomb trong phân tích số bài toán cải tạo đất yếu bằng cọc GEC, giảm thiểu sai số dự báo lún xuống dưới $8.5%$.
  3. Khám phá quy luật phân bố biến dạng phình ngang cục bộ tập trung tại độ sâu $1\text{D} \div 3\text{D}$ từ đỉnh cọc, cung cấp cơ sở tối ưu hóa chiều dài bọc vải địa kỹ thuật trong thực hành kỹ thuật.
  4. Làm sáng tỏ sự sai lệch có hệ thống ($12% \div 22%$) giữa phương pháp giải tích Raithel & Kempfert và phương pháp phần tử hữu hạn, định hướng hệ số an toàn thực tế khi tính toán thiết kế.
  5. Mở ra 3 hướng nghiên cứu mới: Mô hình hóa cố kết thấm thời gian thực theo không gian 3D, phân tích động học cọc GEC dưới tải trọng giao thông lặp, và ứng dụng vật liệu cốt liệu tái chế thân thiện môi trường trong gia cố nền móng tại Việt Nam.