Tổng quan về luận án

Công trình nghiên cứu tiến sĩ của tác giả Lương Nguyễn Hoàng Phương (2023) với đề tài "Nghiên cứu ảnh hưởng của nước đến đặc trưng cơ học của đất không bão hòa áp dụng cho nền đường đắp tại khu vực duyên hải miền Trung" thuộc chuyên ngành Kỹ thuật xây dựng công trình đặc biệt (Mã số: 95.206), Trường Đại học Giao thông Vận tải, dưới sự hướng dẫn khoa học của PGS. Trần Đình Nghiên và TS. Tống Anh Tuấn, đã tạo ra bước tiến quan trọng trong cơ học đất không bão hòa và kỹ thuật địa kỹ thuật hạ tầng giao thông thích ứng biến đổi khí hậu.

Khu vực duyên hải miền Trung Việt Nam là địa bàn chịu tác động khốc liệt của biến đổi khí hậu toàn cầu, nơi có lượng mưa trung bình năm lớn nhất cả nước (vượt mức 2800 mm/năm so với mức 1200–1800 mm/năm ở các vùng khác) và gánh chịu 8–10 trong tổng số 13–16 cơn bão đổ bộ hàng năm. Hiện tượng nắng nóng gay gắt kéo dài làm hạ thấp mực nước ngầm, gia tăng bề dày tầng đất không bão hòa; sau đó các đợt mưa lũ cực đoan với cường độ lớn gây ngập lụt, làm thay đổi cấu trúc lỗ rỗng và gây sạt lở nghiêm trọng các tuyến giao thông huyết mạch, điển hình như sự cố sạt lở taluy âm tại Km 20+315 trên tuyến cao tốc Đà Nẵng – Quảng Ngãi vào tháng 11/2017.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng phần lớn các phân tích ổn định nền đường đắp tại Việt Nam vẫn giả thiết môi trường đất ở trạng thái bão hòa hoàn toàn theo cơ học đất cổ điển của Terzaghi (1936), bỏ qua pha khí và thành phần lực hút dính ma trận $(u_a - u_w)$. Trong khi đó, các nghiên cứu quốc tế về đất không bão hòa chủ yếu tập trung vào mô hình hóa vĩ mô hoặc lý thuyết thấm một chiều, thiếu sự tích hợp đa quy mô giữa phân tích ổn định tổng thể mái dốc và cơ chế xói ngầm, dịch chuyển hạt mịn cục bộ trong cấu trúc ba pha hạt – khí – nước.

Nghiên cứu giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Quy luật biến thiên của đường cong đặc trưng đất – nước (SWCC) và hàm thấm của đất đắp duyên hải miền Trung chịu sự chi phối như thế nào bởi lực hút dính?
    H1: Hàm lượng nước và hệ số thấm giảm phi tuyến tính khi lực hút dính ma trận tăng vượt qua giá trị không khí đi vào (Air Entry Value - AEV).
  • RQ2: Cường độ chống cắt và các tham số góc ma sát trong hiệu quả thay đổi ra sao dưới các cấp áp suất hút dính khác nhau?
    H2: Góc ma sát biểu kiến $\phi^b$ liên hệ với lực hút dính mang tính phi tuyến và suy giảm đáng kể ở dải lực hút dính cao, trong khi góc ma sát hiệu quả $\phi'$ không đổi.
  • RQ3: Sự tương tác đa quy mô giữa dòng thấm không ổn định do mưa và cơ chế xói ngầm vi mô ảnh hưởng như thế nào đến hệ số an toàn mái dốc?
    H3: Mất ổn định nền đường bắt nguồn từ hiện tượng xói ngầm cục bộ ở quy mô lỗ rỗng (DEM-PFV) trước khi hình thành cung trượt tổng thể trong mô hình phần tử hữu hạn (FEM).

Khung lý thuyết của luận án được xây dựng dựa trên lý thuyết biến trạng thái ứng suất độc lập của Fredlund và Morgenstern (1976), mô hình đường cong đặc trưng SWCC của Fredlund và Xing (1994), cùng lý thuyết dòng thấm không bão hòa Darcy mở rộng và cơ học môi trường rời rạc đa pha. Phạm vi thực nghiệm tập trung vào 2 mỏ vật liệu đất đắp đại diện gồm Mỏ Dương Cấm (sét pha cát) và Mỏ Cồn Lê, kết hợp phân tích mô phỏng kiểm chứng tại lý trình Km 24+980 tuyến cao tốc Đà Nẵng – Quảng Ngãi (huyện Đại Lộc, tỉnh Quảng Nam) dưới các kịch bản mưa kéo dài 24h–72h với cường độ $0,12\text{ m/h} - 0,144\text{ m/h}$.

Literature Review và Positioning

Lịch sử phát triển cơ học đất ghi nhận nền tảng từ nguyên lý ứng suất hiệu quả của Terzaghi (1936) với biểu thức $\sigma' = \sigma - u_w$, sau đó Biot (1941) mở rộng lý thuyết cố kết thấm cho môi trường đất hai pha bão hòa. Tuy nhiên, khi áp dụng cho đất không bão hòa tồn tại ba pha, mô hình này bộc lộ nhiều điểm hạn chế. Bishop (1960) đề xuất thông số ứng suất hiệu quả $\chi$ thông qua phương trình $\sigma' = (\sigma - u_a) + \chi(u_a - u_w)$, nhưng Jennings và Burland (1962) đã chứng minh thông số đơn $\chi$ không thể giải thích hiện tượng sụt lún cấu trúc (collapse) khi ngập nước dưới ứng suất chính không đổi.

Để giải quyết tranh luận này, Fredlund và Morgenstern (1976) đã tạo ra bước đột phá khi đề xuất tiếp cận bằng hai tenxơ biến trạng thái ứng suất độc lập gồm ứng suất pháp thực $(\sigma - u_a)$ và lực hút dính ma trận $(u_a - u_w)$. Dòng nghiên cứu về đường cong đặc trưng đất – nước (SWCC) được định hình bởi các mô hình bán thực nghiệm kinh điển của Van Genuchten (1980) và Fredlund & Xing (1994). Leong và Rahardjo (1997) khi so sánh các mô hình SWCC đã chứng minh phương trình của Fredlund & Xing mô tả chính xác nhất mối quan hệ giữa độ ẩm thể tích và lực hút dính nhờ bổ sung hàm hiệu chỉnh $C(\psi)$.

Về cường độ chống cắt, Fredlund, Morgenstern và Widger (1978) đã mở rộng mặt bao phá hoại Mohr–Coulomb thành không gian ba chiều $\tau = c' + (\sigma - u_a)\tan\phi' + (u_a - u_w)\tan\phi^b$. Mặc dù giả định ban đầu xem $\phi^b$ là hằng số, các nghiên cứu thực nghiệm sau đó của Fredlund và cs. (1987), Gan và cs. (1988), và Vanapalli và cs. (1996) đã khẳng định mối quan hệ giữa $\tau$ và $(u_a - u_w)$ là phi tuyến, trong đó $\phi^b$ suy giảm dần khi lực hút dính vượt qua giá trị AEV. Gần đây, Gofar và cs. (2014) cùng Schnellmann và cs. (2023) tiếp tục phát triển các công thức phi tuyến dự báo sức kháng cắt dựa trên độ bão hòa hữu hiệu $S_{re}$ và lực hút dính dư (RSV).

Tại Việt Nam, các nghiên cứu về đất không bão hòa đã được khởi xướng bởi Bạch Quốc Tiến (2009) về dòng thấm liên tục, Phạm Quang Hưng (2012) về ổn định mái dốc dưới mưa phân phối đều tại Hòa Bình, Nguyễn Thị Ngọc Hương (2013) về hiệu chỉnh tham số $m, n$ trong phương trình Fredlund & Xing cho đất đắp ba miền, Phạm Huy Dũng (2020) về ước lượng SWCC cho đất miền Bắc, và Nguyễn Thanh Quang (2021) về suy giảm chỉ tiêu cơ lý đất đập miền Trung theo thời gian.

Công trình của Lương Nguyễn Hoàng Phương (2023) định vị chính xác khoảng trống khoa học bằng cách: (1) thiết lập bộ thông số thực nghiệm đầy đủ về SWCC, hàm thấm và cường độ cắt ba trục CD không bão hòa đặc thù cho dải ven biển miền Trung; (2) so sánh đối chiếu với các mô hình quốc tế tiên tiến như Schnellmann và cs. (2023) và Sibille và cs. (2015); (3) tiên phong tích hợp phương pháp phần tử hữu hạn vĩ mô (FEM qua GeoStudio SEEP/W và SLOPE/W) với phương pháp phần tử rời rạc thể tích lỗ rỗng hữu hạn vi mô (DEM-PFV với $N = 50.000$ hạt) để làm sáng tỏ bản chất phá hoại mái dốc bắt nguồn từ xói ngầm cục bộ.

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

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

Nghiên cứu củng cố và phát triển lý thuyết biến trạng thái ứng suất độc lập của Fredlund và Morgenstern (1976), đồng thời bổ sung bằng chứng thực nghiệm khẳng định tính phi tuyến của góc ma sát lực hút dính $\phi^b$ trên các loại đất nguồn gốc duyên hải miền Trung. Luận án bác bỏ giả định góc $\phi^b$ không đổi trong tính toán truyền thống, chứng minh rằng ở dải lực hút dính thấp dưới giá trị AEV, $\phi^b \approx \phi'$, nhưng khi lực hút dính tăng lên $200\text{ kPa}$, $\phi^b$ suy giảm mạnh về mức $6,12^\circ$ trong khi $\phi'$ giữ nguyên $13,93^\circ$.

Mô hình lý thuyết được chuẩn hóa qua 4 mệnh đề (Propositions):

  1. Mệnh đề P1: Cường độ kháng cắt cực hạn của đất không bão hòa miền Trung tuân theo mặt bao phẳng Mohr–Coulomb mở rộng với độ dốc liên tục giảm theo trục $(u_a - u_w)$, tiến tới trạng thái tiệm cận ở dải hút dính dư.
  2. Mệnh đề P2: Hệ số thấm $k_w$ suy giảm theo hàm mũ đối với lực hút dính ma trận, kiểm soát trực tiếp quá trình dịch chuyển đường bão hòa trong thân nền đường khi xuất hiện gradient áp suất thủy lực.
  3. Mệnh đề P3: Hiện tượng xói ngầm vi mô tuân theo cơ chế mất ổn định cục bộ của các phần tử tứ diện Voronoi/Delaunay khi lực kéo thủy động của dòng thấm có áp vượt qua lực dính mao dẫn tại cổ họng lỗ rỗng (Pore Throats - PT).
  4. Mệnh đề P4: Sự mất ổn định tổng thể của nền đắp là hệ quả tích lũy cấu trúc từ sự gia tăng độ rỗng cục bộ do hạt mịn bị cuốn trôi, dẫn đến tái sắp xếp mạng lưới hạt thô và làm sụt giảm đột ngột hệ số an toàn $F_s$.

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

Khung phân tích của luận án tích hợp liền mạch 3 lý thuyết nền tảng: Lý thuyết cân bằng giới hạn dẻo (Limit Equilibrium Method), Lý thuyết dòng thấm không ổn định qua môi trường rỗng đa pha, và Cơ học môi trường hạt vi mô (Granular Micromechanics).

                                 KHUNG PHÂN TÍCH ĐA QUY MÔ
   [ QUY MÔ VĨ MÔ (MACRO-SCALE) ]                                 [ QUY MÔ VI MÔ (MICRO-SCALE) ]
   - Thí nghiệm áp suất trục (ASTM D6836)                         - Phân tích tương tác 3 pha hạt-khí-nước
   - Thí nghiệm ba trục CD cải tiến                               - Phương pháp DEM-PFV (N = 50.000 hạt)
   - Mô hình hóa FEM: GeoStudio (SEEP/W & SLOPE/W)                - Mạng lưới tứ diện Delaunay & biểu đồ Voronoi
   - Biến trạng thái: (σ - ua) và (ua - uw)                       - Lực mao dẫn & áp suất dòng thấm tại cổ họng rỗng (PT)
                             [ ĐÁNH GIÁ ỔN ĐỊNH TOÀN DIỆN ]
                             - Cơ chế xói ngầm sinh lỗ rỗng cục bộ
                             - Suy giảm cường độ chống cắt theo mưa
                             - Tiến triển hệ số an toàn Fs theo thời gian

Cách tiếp cận phân tích vi mô sử dụng mô hình thể tích tỉ lệ rỗng hữu hạn (Pore-scale Finite Volume - PFV) ghép nối với phương pháp phần tử rời rạc (DEM) trên nền tảng phân rã tam giác Delaunay và ô mạng Voronoi. Khung phân tích xác định rõ các điều kiện biên: bề mặt thoáng chịu phân phối mưa theo thời gian thực, biên thấm có áp thay đổi theo mực nước ngầm thượng/hạ lưu, và các biên cứng giả định tiếp xúc hạt.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái nhận thức luận 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 định lượng thực nghiệm đa cấp (Multi-level quantitative experimental design). Thiết kế kết hợp giữa thí nghiệm cơ học đất tiêu chuẩn nâng cao trong phòng thí nghiệm và mô phỏng số hai cấp độ: quy mô công trình (Macro-scale FEM) và quy mô hạt (Micro-scale DEM-PFV).

Mẫu đất nghiên cứu được thu thập từ hai mỏ đất đắp tiêu biểu tại duyên hải miền Trung phục vụ thi công cao tốc Đà Nẵng – Quảng Ngãi:

  • Mỏ Dương Cấm: Loại đất sét pha cát, chỉ số dẻo $I_p$ trung bình, dung trọng khô tiêu chuẩn đầm nén lớn nhất $\gamma_{d,\max}$, độ ẩm tối ưu $w_{opt}$.
  • Mỏ Cồn Lê: Đất đắp có cấp phối hạt chọn lọc, hàm lượng hạt thô chiếm tỷ lệ cao hơn.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu tuân thủ nghiêm ngặt các tiêu chuẩn quốc tế và thiết bị hiện đại:

  • Xác định đường cong SWCC: Sử dụng thiết bị buồng chiết suất áp suất cao với đĩa gốm sứ tiếp nhận áp suất khí cao (High-Air-Entry Ceramic Disk) theo tiêu chuẩn ASTM D2325-68 và ASTM D6836-02. Ứng dụng kỹ thuật tịnh tiến trục (Axis Translation Technique), duy trì áp suất nước lỗ rỗng $u_w = 0$ và tác dụng các cấp áp suất khí dương $u_a$ tăng dần, đo lường chính xác lượng nước thoát qua đĩa gốm tại mỗi trạng thái cân bằng.
  • Thí nghiệm cắt trực tiếp (DST) cải tiến: Xác định nhanh mặt bao phá hoại Mohr–Coulomb mở rộng ở các cấp lực hút dính khác nhau.
  • Thí nghiệm nén ba trục cố kết thoát nước + khí (CD) cải tiến: Thực hiện trên hệ thống buồng nén ba trục đặc biệt tại Phòng thí nghiệm Địa kỹ thuật (Đại học Thủy lợi). Quy trình gồm 4 giai đoạn nghiêm ngặt: (1) Bão hòa mẫu và đĩa gốm; (2) Cố kết đẳng hướng dưới áp suất buồng; (3) Tạo và cân bằng lực hút dính ma trận $(u_a - u_w)$ mục tiêu (0, 100, 200 kPa); (4) Cắt mẫu với tốc độ biến dạng chậm nhằm đảm bảo tiêu tán hoàn toàn áp suất nước và khí lỗ rỗng dư.
                    QUY TRÌNH THÍ NGHIỆM BA TRỤC CD KHÔNG BÃO HÒA
[ Giai đoạn 1: Bão hòa ]       [ Giai đoạn 2: Cố kết ]      [ Giai đoạn 3: Cân bằng Suction ]
Bão hòa mẫu và đĩa gốm áp    Cố kết đẳng hướng mẫu dưới    Tác dụng áp suất khí ua và nước uw
suất cao (High AEV disk)      áp suất buồng hiệu dụng        để thiết lập (ua - uw) = 0, 100, 200 kPa
                         [ Giai đoạn 4: Gia tải cắt chậm ]
                         Thoát nước và khí tự do (CD)
                         Đo biến dạng dọc trục và ứng suất lệch đỉnh

Quy trình kiểm chứng độ tin cậy được thực hiện qua tam giác hóa phương pháp (Triangulation): đối chiếu số liệu thí nghiệm ba trục CD thực tế với công thức dự báo thực nghiệm của Fredlund & Vanapalli (1996), so sánh nghiệm giải tích thấm của Terzaghi với kết quả mô phỏng DEM-PFV.

Data và phân tích

Dữ liệu thực nghiệm được xử lý và mô phỏng trên các nền tảng tính toán chuyên sâu:

  • Phân tích hàm SWCC và độ thấm: Sử dụng thuật toán khớp phi tuyến theo phương trình Fredlund & Xing (1994) để xác định 4 tham số $a, n, m, C(\psi)$ với giá trị lực hút dính giới hạn $10^6\text{ kPa}$, từ đó tích phân xác định hàm hệ số thấm phi tuyến $k_w(u_a - u_w)$.
  • Mô phỏng phần tử hữu hạn (FEM): Phần mềm GeoStudio kết hợp hai module SEEP/W (phân tích trường áp suất nước lỗ rỗng không ổn định theo thời gian) và SLOPE/W (phân tích ổn định mái dốc theo phương pháp cân bằng giới hạn Morgenstern-Price) cho mặt cắt nền đắp cao tại Km 24+980 cao tốc Đà Nẵng – Quảng Ngãi.
  • Mô phỏng phần tử rời rạc kết hợp thể tích rỗng hữu hạn (DEM-PFV): Sử dụng hệ thống $N = 50.000$ hạt cầu, thiết lập các liên kết tiếp xúc phi tuyến theo định luật II Newton, giải bài toán dòng thấm có áp cục bộ và cố kết hạt mịn.

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

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

Nghiên cứu mang lại 5 phát hiện cốt lõi mang tính định lượng cao:

  1. Sự phụ thuộc phi tuyến của cường độ chống cắt vào lực hút dính:
    Đối với mẫu đầm nén mỏ Dương Cấm, ở trạng thái bão hòa hoàn toàn $(u_a - u_w = 0\text{ kPa})$, các thông số hiệu quả đạt lực dính $c' = 14,4\text{ kPa}$ và góc ma sát trong $\phi' = 13,93^\circ$. Khi lực hút dính ma trận tăng lên các cấp $100\text{ kPa}$ và $200\text{ kPa}$, lực hút dính thực tế đo được tại thời điểm phá hoại (ứng suất lệch đỉnh) đạt giá trị trung bình lần lượt là $41,17\text{ kPa}$ và $66,88\text{ kPa}$. Góc ma sát biểu kiến $\phi^b$ suy giảm phi tuyến mạnh mẽ từ $\phi^b = 13,93^\circ$ (tại $0\text{ kPa}$) xuống $\phi^b = 6,12^\circ$ (tại $200\text{ kPa}$), trong khi góc ma sát trong hiệu quả $\phi'$ bảo toàn ở mức $13,93^\circ$.

  2. Dạng đường cong SWCC đặc thù của đất duyên hải miền Trung:
    Đường cong SWCC của đất mỏ Dương Cấm và Cồn Lê phản ánh rõ rệt 3 vùng trạng thái: vùng hiệu ứng biên có giá trị không khí đi vào (AEV) tương đối thấp do tỷ lệ hạt cát pha sét đặc trưng, vùng chuyển tiếp có độ dốc cao cho thấy tốc độ khử bão hòa diễn ra nhanh chóng, và vùng tàng dư đạt tới giá trị lực hút dính dư (RSV) trong khoảng $1500 - 3000\text{ kPa}$.

  3. Mối tương quan giữa mô hình phân phối mưa và sự suy giảm hệ số ổn định $F_s$:
    Kết quả mô phỏng ghép nối SEEP/W và SLOPE/W cho thấy mô hình mưa phân phối chuẩn (mưa tập trung cường độ đỉnh ở giữa chu kỳ) gây ra mức độ suy giảm hệ số an toàn mái dốc $F_s$ nghiêm trọng hơn nhiều so với mô hình mưa phân phối đều với cùng tổng lượng mưa ($0,12\text{ m/h}$ và $0,144\text{ m/h}$ trong 72h). Đỉnh điểm mất ổn định xuất hiện khi mặt thấm bão hòa dâng cao từ chân taluy kết hợp nước mưa thấm từ đỉnh và vai đường làm triệt tiêu hoàn toàn áp suất âm $(u_a - u_w)$ tại vùng trượt nguy hiểm.

  4. Làm sáng tỏ cơ chế xói ngầm vi mô bằng mô hình DEM-PFV:
    Mô phỏng vi mô với 50.000 hạt đã chứng minh rằng khi gradient thủy lực cục bộ vượt ngưỡng tới hạn, dòng thấm có áp bứt các hạt mịn ra khỏi khung xương chịu lực tại các cổ họng rỗng (PT). Hiện tượng này làm gia tăng đột biến độ rỗng cục bộ, tạo dòng chảy tập trung dạng đường ống (piping) trước khi mái dốc biểu hiện biến dạng trượt vĩ mô.

  5. Sự phù hợp giữa mô hình nén 3 trục CD và dự báo lý thuyết:
    Mối quan hệ giữa ứng suất cắt $\tau$ và lực hút dính $(u_a - u_w)$ thu được từ thí nghiệm ba trục cải tiến CD tương thích cao với phương trình thực nghiệm của Fredlund & Vanapalli (1996), khẳng định độ chính xác của hệ thống thiết bị chế tạo tại Việt Nam.

                  BIẾN THIÊN CÁC THÔNG SỐ ĐỊA KỸ THUẬT THEO LỰC HÚT DÍNH

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở khoa học để hiệu chỉnh các mô hình dự báo sức kháng cắt phi tuyến trong điều kiện khí hậu nhiệt đới gió mùa ven biển.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn phối hợp mô phỏng đa quy mô DEM-PFV và FEM trong phân tích tương tác đất – nước, có thể áp dụng cho các bài toán đập đất, sạt lở đồi núi và đê biển.
  • Về mặt thực tiễn công trình: Cho phép kỹ sư địa kỹ thuật tối ưu hóa kích thước hình học mái dốc taluy nền đắp, xác định chính xác vị trí cần bố trí hệ thống rãnh ngầm thu nước và tầng lọc ngược chống xói ngầm.
  • Về mặt chính sách: Đề xuất cơ quan quản lý giao thông (Bộ GTVT) bổ sung quy chuẩn tính toán ổn định công trình giao thông có xét đến điều kiện không bão hòa và biến đổi khí hậu tại dải ven biển miền Trung.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận 4 giới hạn:

  1. Giả thiết hình học trong mô hình rời rạc DEM hiện tại xem các hạt đất là hình cầu lý tưởng, chưa mô phỏng được độ góc cạnh và hình dạng phức tạp của hạt khoáng sét tự nhiên.
  2. Năng lực phần cứng tính toán giới hạn số lượng hạt mô phỏng DEM-PFV ở mức $N = 50.000$ hạt, chưa thể mở rộng mô hình vi mô này trực tiếp ra toàn bộ thể tích nền đắp hàng chục mét khối.
  3. Thí nghiệm ba trục CD chưa khảo sát hiện tượng trễ (hysteresis) giữa quá trình làm ướt (wetting) và làm khô (drying) của đường cong SWCC dưới các chu kỳ bốc hơi – ngập lụt lặp lại nhiều lần.
  4. Chưa tích hợp đầy đủ tác động của áp lực nước biển dâng kết hợp triều cường vào điều kiện biên thủy lực chân mái dốc phía hạ lưu.

Chương trình nghiên cứu 5–10 năm tới mở ra các hướng đi mới:

  • Phát triển phần tử DEM đa hình thái (non-spherical particles) để mô phỏng chính xác tương tác tiếp xúc góc cạnh.
  • Xây dựng mô hình ghép nối song phương liên tục – rời rạc (Continuous-Discontinuous Coupled FEM-DEM) chạy trên hệ thống siêu máy tính hiệu năng cao (GPU/Parallel Computing).
  • Nghiên cứu cơ chế trễ SWCC dưới tác động của biến đổi nhiệt độ và chu kỳ khô – ướt cực đoan mô phỏng biến đổi khí hậu dài hạn.
  • Mở rộng ngân hàng dữ liệu thực nghiệm SWCC và cường độ cắt không bão hòa cho toàn bộ các mỏ vật liệu đất đắp dọc tuyến cao tốc Bắc – Nam qua 14 tỉnh duyên hải miền Trung.

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

  • Ảnh hưởng học thuật: Luận án tạo tiền đề cho các công bố quốc tế uy tín trong danh mục ISI/Scopus về địa kỹ thuật môi trường, ước tính mang lại hàng trăm trích dẫn học thuật nhờ giải quyết khoảng trống thực nghiệm đất không bão hòa miền Trung và mô hình hóa DEM-PFV.
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp cho các tổng công ty tư vấn thiết kế giao thông (TEDI, TEDIS), giúp tối ưu hóa khối lượng đào đắp từ 5–10% trong khi vẫn đảm bảo tuyệt đối an toàn công trình trước thiên tai.
  • Tác động kinh tế – xã hội: Giảm thiểu nguy cơ đứt gãy mạch máu giao thông do sạt lở nền đường sau mưa bão, tiết kiệm hàng trăm tỷ đồng chi phí sửa chữa khắc phục sự cố hàng năm tại các tỉnh miền Trung.
  • Ý nghĩa quốc tế: Đóng góp dữ liệu thực nghiệm quý giá từ khu vực nhiệt đới gió mùa châu Á vào ngân hàng dữ liệu cơ học đất không bão hòa toàn cầu của Hội Địa kỹ thuật Quốc tế (ISSMGE TC106).

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận quy trình vận hành chi tiết buồng thí nghiệm ba trục cải tiến CD cho đất không bão hòa và thuật toán mô phỏng DEM-PFV mã nguồn mở.
  • Giảng viên & Chuyên gia Địa kỹ thuật cao cấp: Có thêm tài liệu tham khảo chuyên sâu và bộ tham số thực tế phục vụ công tác giảng dạy sau đại học và cố vấn công trình trọng điểm.
  • Kỹ sư R&D & Tư vấn thiết kế hạ tầng: Sở hữu phương pháp luận kết hợp SEEP/W và SLOPE/W với đường cong SWCC thực đo, nâng cao độ chính xác so với phương pháp tính toán bão hòa truyền thống.
  • Nhà hoạch định chính sách & Ban quản lý dự án: Cơ sở khoa học vững chắc để phê duyệt định mức kỹ thuật, tiêu chuẩn thiết kế nền đường cao tốc thích ứng kịch bản biến đổi khí hậu và nước biển dâng đến năm 2050.

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 làm sáng tỏ quy luật suy giảm phi tuyến của góc ma sát lực hút dính $\phi^b$ khi lực hút dính tăng cao trên nền đất sét pha cát duyên hải miền Trung, mở rộng trực tiếp lý thuyết cường độ chống cắt ba chiều của Fredlund, Morgenstern và Widger (1978). Luận án chứng minh $\phi^b$ không phải là hằng số mà suy giảm từ $13,93^\circ$ về $6,12^\circ$ ở dải hút dính $200\text{ kPa}$, xác nhận tính đúng đắn của các mô hình dự báo phi tuyến hiện đại.

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

So với nghiên cứu của Phạm Quang Hưng (2012) và Nguyễn Thị Ngọc Hương (2013) vốn chỉ dùng mô hình phần tử hữu hạn vĩ mô hoặc thí nghiệm cơ học truyền thống, luận án tiên phong kết hợp hai quy mô: mô hình hóa vĩ mô FEM (GeoStudio SEEP/W & SLOPE/W) với mô phỏng vi mô phần tử rời rạc thể tích rỗng hữu hạn DEM-PFV ($N = 50.000$ hạt). Phương pháp này cho phép theo dõi trực tiếp dòng thấm cục bộ và cơ chế cuốn trôi hạt mịn gây xói ngầm mà FEM thuần túy không thể mô tả được.

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

Phát hiện bất ngờ nhất là sự bất đối xứng nguy hiểm của mô hình phân phối mưa: cùng một lượng mưa tích lũy và cường độ đỉnh, mô hình mưa phân phối chuẩn (tập trung mưa lớn ở giữa đợt) gây suy giảm hệ số an toàn mái dốc $F_s$ nhanh và sâu hơn nhiều so với mô hình mưa phân phối đều. Hơn nữa, quá trình phá hoại tổng thể luôn được kích hoạt trước bởi hiện tượng gia tăng độ rỗng cục bộ do xói ngầm vi mô ở chân taluy.

4. Luận án có cung cấp quy trình tái lập thực nghiệm chi tiết không?

Có. Luận án mô tả chi tiết quy trình chế bị mẫu đầm nén, kỹ thuật tịnh tiến trục theo ASTM D6836-02, các bước bão hòa đĩa gốm áp suất cao và quy trình 4 bước vận hành buồng nén ba trục cải tiến CD tại PTN Đại học Thủy lợi, đảm bảo đầy đủ điều kiện để các nhà nghiên cứu độc lập tái lập thực nghiệm.

5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình hướng tới: (1) Ứng dụng trí tuệ nhân tạo (Machine Learning) kết hợp mạng thần kinh vật lý (PINNs) để dự báo tức thời SWCC từ các chỉ tiêu cơ lý đơn giản; (2) Phát triển phần mềm tính toán ghép nối DEM-FEM song song trên GPU cho toàn bộ khối đắp công trình; (3) Mở rộng tiêu chuẩn thiết kế nền đường thích ứng biến đổi khí hậu trên quy mô toàn quốc.

Kết luận

  1. Luận án đã xác định thành công bộ thông số cơ học đất không bão hòa đặc trưng (SWCC, hàm thấm phi tuyến, $c' = 14,4\text{ kPa}$, $\phi' = 13,93^\circ$, $\phi^b$ phi tuyến) cho đất đắp mỏ Dương Cấm và Cồn Lê tại duyên hải miền Trung.
  2. Xây dựng và kiểm chứng hoàn thiện quy trình thí nghiệm nén ba trục CD cải tiến và kỹ thuật tịnh tiến trục đo lực hút dính ma trận chính xác cao.
  3. Thiết lập mô hình mô phỏng vi mô DEM-PFV với $50.000$ hạt, giải mã cơ chế xói ngầm và biến dạng thể tích cục bộ trong cấu trúc ba pha hạt – khí – nước.
  4. Đánh giá toàn diện ổn định nền đường đắp cao tốc Đà Nẵng – Quảng Ngãi (Km 24+980) dưới các kịch bản mưa phân phối chuẩn và mưa đều cực đoan kéo dài 72h.
  5. Chứng minh vai trò quyết định của lực hút dính ma trận trong việc duy trì ổn định mái dốc và cảnh báo nguy cơ sạt lở khi đường bão hòa dâng cao làm triệt tiêu áp suất nước lỗ rỗng âm.
  6. Mở ra 3 hướng nghiên cứu mới về cơ học môi trường rời rạc đa quy mô, địa kỹ thuật thích ứng biến đổi khí hậu toàn cầu và tiêu chuẩn hóa thiết kế công trình giao thông ven biển.