Tổng quan về luận án

Bảo vệ và phát triển bền vững dải ven biển trước tác động cực đoan của biến đổi khí hậu và hiện tượng nước biển dâng là nhiệm vụ địa kỹ thuật chiến lược tại Việt Nam. Hệ thống đê biển Việt Nam có quy mô đồ sộ với tổng chiều dài xấp xỉ 2.700 km (bao gồm khoảng 1.400 km đê trực diện biển và 1.300 km đê cửa sông ven biển), kéo dài từ Quảng Ninh đến Kiên Giang. Do điều kiện kinh tế - kỹ thuật và đặc thù hình thành tự nhiên, các tuyến đê biển nước ta phần lớn được thiết kế dưới dạng "công trình bán vĩnh cửu" (chu kỳ quai đê lấn biển 15–20 năm với diện tích 1.000–2.000 ha/lần tại các vùng bồi tụ như Hải Phòng, Thái Bình, Nam Định, Ninh Bình; hoặc chịu tốc độ xói lở nghiêm trọng lên tới 17 m/năm như tại Hải Hậu - Nam Định). Thách thức cốt lõi là thân đê và nền đê chủ yếu cấu tạo từ nền đất yếu, vật liệu tại chỗ rời rạc, hạt mịn (bùn sét, cát pha, á cát, á sét có cấp phối hạt $0,005 \div 0,5\text{ mm}$), góc ma sát trong và lực dính rất bé, dễ mất ổn định khi chịu sóng tràn với lưu tốc lên đến $4\text{ m/s}$ hoặc khi xuất hiện áp lực nước lỗ rỗng tăng vọt.

Khoảng trống học thuật (Research Gap) trung tâm được xác định: Trong lý thuyết phân thỏi (Method of Slices) tính toán ổn định mái dốc công trình đất, bài toán xác định hệ số an toàn $F_s$ là bài toán siêu tĩnh bậc cao. Với $n$ thỏi đất trượt, số lượng ẩn số lực và điểm đặt lực cần tìm lên tới $(6n - 2)$ đại lượng, trong khi hệ phương trình cân bằng tĩnh học chỉ thiết lập được $4n$ phương trình, dẫn đến sự thiếu hụt $(2n - 2)$ phương trình độc lập. Để giải quyết, các phương pháp kinh điển (Fellenius, Bishop đơn giản, Janbu, Spencer, Morgenstern - Price, General Limit Equilibrium - GLE) đều phải đưa ra các giả thiết đơn giản hóa hoặc gán trước quy luật hàm góc nghiêng lực tương tác $\lambda f(x)$, hoặc giả định đường tương tác lực. Khi đưa cốt vải địa kỹ thuật (VĐKT) vào gia cố thân đê, các mô hình tính toán hiện hành thường chỉ coi lực neo của cốt như một ngoại lực tập trung cục bộ hoặc đưa vào dưới dạng mômen chống trượt thuần túy ($M_{V/0} = V_i \cdot R$), mà hoàn toàn bỏ qua điều kiện tương thích lực tương tác giữa các thỏi đất lân cận có xét đến sự phân bố lại ứng suất nội tại do lực căng cốt vải truyền vào đất.

Luận án tiến sĩ kỹ thuật chuyên ngành Xây dựng Công trình Thủy (Mã số: 62 58 40 01) của NCS. Phan Tiến An dưới sự hướng dẫn khoa học của GS.TS. Vũ Đình Hùng tại Viện Khoa học Thủy lợi Việt Nam (2011) đã giải quyết trọn vẹn khoảng trống này với ba câu hỏi nghiên cứu (RQs) và ba giả thuyết khoa học ($H_s$) cụ thể:

  • RQ1: Làm thế nào để tĩnh định hóa bài toán phân tích ổn định mái dốc phân thỏi mà không cần áp đặt các giả thiết nhân tạo, phi vật lý về góc nghiêng hay vị trí điểm đặt lực tương tác?
    • H1: Việc tích hợp lý thuyết phân tích hệ thống kết hợp điều kiện tương thích lực đẩy - lực chống trượt giữa các thỏi đất liền kề (dựa trên định lý Gvozdev và nguyên lý cực trị Coulomb) cho phép khép kín hệ phương trình tĩnh học một cách tường minh và duy nhất.
  • RQ2: Cơ chế tương tác truyền ứng suất và huy động lực kéo neo của lớp cốt VĐKT bên trong khối trượt làm thay đổi trạng thái cân bằng giới hạn của từng thỏi đất như thế nào?
    • H2: Lực kéo căng cốt vải $V_i$ không chỉ tạo mômen giữ tổng thể mà trực tiếp tái phân bố trường số gia lực pháp tuyến $\Delta E_i$ và số gia lực tiếp tuyến $\Delta X_i$ giữa các thỏi, từ đó nâng cao cục bộ và toàn cục hệ số kháng trượt $F_s$.
  • RQ3: Việc áp dụng công nghệ đất có cốt VĐKT kết hợp phương pháp tính tương thích có khả thi và kiểm chứng được trên mô hình vật lý quy mô lớn và công trình thực tế hay không?
    • H3: Mô hình số được lập trình chuyên dụng sẽ hội tụ ổn định, khớp nối chính xác với kiểm chứng phần mềm thương mại GeoSlope (Slope/W) và phản ánh sát thực chuyển vị ngang đo đạc tại mô hình vật lý 1:1 cũng như đoạn đê biển thử nghiệm KCT2 Bình Minh 3 (Kim Sơn, Ninh Bình).

Khung lý thuyết nền tảng tích hợp chặt chẽ giữa: (1) Cơ học đất dẻo Mohr - Coulomb; (2) Lý thuyết cân bằng giới hạn phân thỏi; (3) Lý thuyết phân tích hệ thống tối ưu hóa trong địa kỹ thuật của GS. Phan Trường Phiệt và TS. Phan Trường Giang; (4) Lý thuyết đất có cốt của Henri Vidal (1963) và Schlosser (1972). Quy mô kiểm chứng trải rộng từ mô hình giải tích, mô hình vật lý tỷ lệ 1:1 đo đạc chuyển vị ngang ($S_n$), thí nghiệm hộp cắt phẳng kéo tuột cốt vải xác định hệ số tương tác ma sát - dính ($f_{ds}, f_{po}$), cho đến thi công thực nghiệm toàn diện trên tuyến đê biển cấp bách chống bão lũ tại bờ biển Ninh Bình.

Literature Review và Positioning

Lịch sử phát triển của các phương pháp phân tích ổn định mái dốc theo trường phái cân bằng giới hạn (Limit Equilibrium Methods - LEM) đã định hình nền tảng địa kỹ thuật công trình thế giới qua hơn một thế kỷ:

  • Fellenius (1936): Đơn giản hóa triệt để bằng cách bỏ qua toàn bộ lực tương tác giữa các thỏi ($E_i = X_i = 0$), chỉ xét cân bằng của từng thỏi riêng lẻ. Phương pháp này chỉ mang giá trị lịch sử vì đánh giá hệ số $F_s$ quá thiên về an toàn sai lệch (dưới ước tính sức chịu tải thực tế từ 30% đến 50%).
  • Bishop đơn giản (1955): Thừa nhận có lực đẩy ngang $E_i$ nhưng triệt tiêu thành phần lực cắt thẳng đứng giữa các thỏi ($X_i = 0$). Dù thỏa mãn điều kiện cân bằng mômen tổng thể và cho kết quả chấp nhận được trong mái dốc đồng nhất, phương pháp Bishop bỏ qua tương tác cắt giữa các thỏi, dẫn đến sai số lớn khi mái đê có phân lớp phức tạp hoặc gia cường cốt chịu kéo nằm ngang.
  • Spencer (1967) và Morgenstern - Price (1965): Khôi phục cả hai thành phần $E_i$ và $X_i$ thông qua giả thiết góc nghiêng lực tương tác không đổi ($X/E = \tan\theta = \text{const}$) ở Spencer, hoặc hàm góc biến thiên tùy biến $X/E = \lambda f(x)$ ở Morgenstern - Price.
  • Janbu tổng quát (1973): Thiết lập hệ 5 phương trình chứa 6 ẩn số ($E_i, X_i, h_{pi}, N_i, T_i, F_s$) bằng cách định vị đường tác dụng lực tương tác (đường đẩy $h_{pi}$). Tuy nhiên, theo phân tích của Fredlund (1977, 1981), phương pháp Janbu tổng quát dù chặt chẽ về lý thuyết nhưng thuật toán số cực kỳ khó hội tụ và rất nhạy cảm với việc chọn hình học đường lực.
  • Phương pháp cân bằng giới hạn tổng quát GLE (Fredlund & Krahn, 1977): Tích hợp cả Spencer và Morgenstern - Price bằng cách giải lặp hai hệ số an toàn cân bằng mômen ($F_m$) và cân bằng lực ($F_f$) dựa trên hàm $f(x) = \sin(x)$. Dù phổ biến trong các bộ phần mềm thương mại (GeoStudio Slope/W), GLE đòi hỏi quá trình thử dần phức tạp và tiềm ẩn nguy cơ không hội tụ nghiệm thực nếu hàm $f(x)$ giả định không tương thích với trường ứng suất biến dạng thực tế.
                         TIẾN TRÌNH LÝ THUYẾT PHÂN THỎI

Trong lĩnh vực đất có cốt, kể từ phát minh của Henri Vidal (1963) và các thí nghiệm ba trục xác định ứng xử đất cốt đồng nhất dị hướng của Schlosser & Long (1972), thế giới đã phân nhánh ứng dụng rõ rệt:

  1. Dạng đất trộn cốt rời: Sử dụng các đoạn sợi polyme, xơ dừa hoặc mảnh lưới địa kỹ thuật phân tán đẳng hướng (được Giroud phát triển tiếp thành công nghệ túi đất có đuôi "con nòng nọc đất" thi công chống xói lở lấp hàm ếch dưới nước).
  2. Dạng kết cấu đất đặt cốt định hình dạng màng/lưới: Tiêu biểu với các hướng dẫn thiết kế tường chắn và mái dốc của Schmertmann (1987), Broms (1980), và Tiêu chuẩn Cục Đường bộ Liên bang Hoa Kỳ (FHWA).

So sánh với các nghiên cứu công trình biển quốc tế tiêu biểu:

  • Đê biển Deep Bay (Hồng Kông): Chiều dài 3,5 km thi công trên bùn sét siêu yếu, sử dụng VĐKT để phân cách, ngăn lún chìm và gia cố trượt tổng thể trong thời gian ngắn.
  • Dự án lấn biển cải tạo vùng nước sâu Laem Chabang (Thái Lan): Kết hợp VĐKT cường độ cao với lưới tre tạo đệm cát gia cố nền sét bùn chịu nén trước thủy lực.
  • Đập đất Ait Chouarit (Ma-rốc) & Đê quai Tam Hiệp (Sanxia - Trung Quốc cao 90m): Ứng dụng hàng chục nghìn $m^2$ VĐKT làm tầng lọc và cốt gia cường chịu áp lực thủy động cực lớn tồn tại an toàn suốt hơn 10 năm thi công.
  • Tuyến đê bảo vệ bờ New Orleans (Hoa Kỳ): Trải thảm địa kỹ thuật dệt Nicolon để giảm tiết diện đắp, ngăn chặn hiện tượng mất đất lún chìm vào nền bùn lỏng.

Tại Việt Nam, các công trình của GS. Phan Trường Phiệt (1996, 1999), PGS. Vũ Đình Hùng (2001, 2005), PGS. Nguyễn Quốc Dũng (2006), TS. Phạm Văn Long (2004) đã bước đầu đặt nền móng ứng dụng VĐKT làm tầng lọc, chống xói mái, hoặc tường chắn đất. Tuy nhiên, chưa có công trình nào giải quyết thấu đáo bài toán tương thích tĩnh học phân thỏi khi trong lòng từng lát cắt có sự hiện diện của cốt neo chịu kéo làm biến đổi cục bộ hướng và độ lớn của cặp véc-tơ nội lực $(\Delta E_i, \Delta X_i)$. Luận án của Phan Tiến An định vị chính xác tại điểm giao thoa này, chuyển dịch mô hình tính từ việc "áp đặt thông số hình học giả định" sang "tự cân bằng tương thích theo nguyên lý cơ học cực trị".

Đó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à hoàn thiện lý thuyết phân thỏi cổ điển thông qua việc xác lập cơ chế tương thích tĩnh học cho hệ lực tương tác giữa các thỏi đất:

  1. Khắc phục triệt để tính siêu tĩnh phi vật lý: Thay vì cưỡng bức gán góc nghiêng lực $\theta$ hoặc đường tác dụng lực $h_{pi}$, luận án xem hai phần khối trượt phân chia bởi một lát cắt thẳng đứng bất kỳ như một hệ thống đẩy - chống tương hỗ, tương tự hệ thống "Tường chắn - Khối đất sau tường" ở trạng thái cân bằng giới hạn.
  2. Vận dụng định lý Gvozdev và nguyên lý cực trị Coulomb:
    • Khi phần khối đất phía trên đẩy phần khối đất phía dưới: Lực đẩy chủ động đạt cực đại ($E_a^{\max}$) tương ứng với trạng thái giới hạn chủ động.
    • Khi khối đất đóng vai trò cản trượt: Lực cản thụ động là lực đẩy nhỏ nhất ($E_p^{\min}$) cần thiết để duy trì trạng thái cân bằng giới hạn bị động.
    • Dạng phá hoại trượt thực tế của hệ thống tương ứng với trị số nhỏ nhất của tải trọng phụ phá hoại.
  3. Mô hình hóa số gia lực tương tác $(\Delta E_i, \Delta X_i)$ có xét cốt VĐKT: Lực kéo neo của cốt vải ($V_i$) tại điểm giao cắt với đáy thỏi trượt hoặc thân thỏi trượt được phân tích giải tích thành các thành phần tham gia trực tiếp vào đa giác lực khép kín của từng phân tố: $$\sum F_x = 0 \iff E_{i-1} - E_i + T_i \cos\alpha_i - N_i \sin\alpha_i + V_i \cos\beta_i = 0$$ $$\sum F_y = 0 \iff W_i + X_{i-1} - X_i - N_i \cos\alpha_i - T_i \sin\alpha_i - V_i \sin\beta_i = 0$$ với điều kiện bền Mohr - Coulomb tiếp xúc đáy thỏi: $$T_i = \frac{1}{F_s} \left( N_i - u_i l_i \right) \tan\varphi^* + \frac{c^* l_i}{F_s}$$
                ĐA GIÁC LỰC TÁC DỤNG LÊN MỘT THỎI ĐẤT (CÓ CỐT VĐKT)
                                 Wi (Trọng lượng)
                         Ni (Phản lực pháp tuyến)
                       Ti (Lực chống trượt đáy thỏi)

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

Khung phân tích của luận án tích hợp đa tầng ba trụ cột lý thuyết:

  • Trụ cột 1: Cơ học tiếp xúc đất - cốt màng mỏng: Khảo sát đồng thời hai cơ chế phá hoại: (i) Sự trượt của đất trên bề mặt cốt vải (xác định qua hệ số ma sát tiếp xúc $f_{ds}$); (ii) Sự kéo tuột cốt vải ra khỏi khối đất kẹp (xác định qua hệ số kháng nhổ $f_{po}$).
  • Trụ cột 2: Tối ưu hóa điều khiển hệ thống: Thiết lập phiếm hàm mục tiêu tìm hệ số an toàn nhỏ nhất $F_s^{\min}$ trên tập hợp các cung trượt tròn khả dĩ tâm $O(x, y)$, bán kính $R$, thỏa mãn đồng thời các phương trình cân bằng hình chiếu lực và mômen của toàn hệ thống $n$ thỏi mà không làm vi phạm điều kiện tương thích lực biên: $$\Delta(F_s, \lambda) = 0$$
  • Trụ cột 3: Điều kiện biên giới hạn (Boundary Conditions):
    • Tại đỉnh mái dốc (thỏi số 1): $E_0 = 0, X_0 = 0$.
    • Tại chân mái dốc (thỏi số $n$): $E_n = 0, X_n = 0$.
    • Ứng suất kéo cho phép và biến dạng giới hạn của VĐKT: $V_i \le [T_{\text{allow}}]$, với $[T_{\text{allow}}] = T_{\text{ult}} / (RF_{\text{ID}} \times RF_{\text{CR}} \times RF_{\text{CD}} \times RF_{\text{BD}})$.

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ thế giới quan thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phê phán (Critical Realism) trong cơ học địa kỹ thuật. Phương pháp nghiên cứu là sự kết hợp đa phương thức (Multi-method Triangulation):

  • Phân tích giải tích toán học: Thiết lập hệ phương trình vi phân sai phân cân bằng lực tương tác có điều kiện tương thích.
  • Thực nghiệm trong phòng: Thí nghiệm cắt phẳng quy mô hộp cắt cải tiến và thí nghiệm kéo tuột để xác định lực dính đơn vị $c$, góc ma sát trong $\varphi$, lực dính tiếp xúc $c_{ds}$, và góc ma sát tiếp xúc $\varphi_{ds}$.
  • Mô hình vật lý tỷ lệ lớn 1:1: Thử nghiệm hiện trường với kích thước thật nhằm loại bỏ hoàn toàn hiệu ứng tỷ lệ (Scale Effects) thường gặp ở mô hình ly tâm hay mô hình thu nhỏ.
  • Mô phỏng số học (Computational Modeling): Tự viết chương trình phần mềm chuyên dụng và đối chiếu chéo (Cross-verification) với phần mềm chuẩn quốc tế GeoStudio (Slope/W của GeoSlope International, Canada).
                      SƠ ĐỒ PHƯƠNG PHÁP NGHIÊN CỨU ĐA HỆ

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được kiểm soát nghiêm ngặt theo các tiêu chuẩn kỹ thuật thủy lợi và địa kỹ thuật (TCVN, ASTM, ISO):

  1. Quy trình lấy mẫu đất: Thu thập mẫu nguyên dạng và không nguyên dạng tại nền đê và mỏ vật liệu đắp ven biển Kim Sơn (Ninh Bình). Đất thuộc nhóm sét pha màu xám nâu, xám xanh và cát pha lơ lửng có hàm lượng mùn hữu cơ, đường kính hạt biến thiên từ $0,005\text{ mm}$ đến $0,5\text{ mm}$.
  2. Giao thức thí nghiệm cắt và kéo rút cốt vải: Tiến hành trên thiết bị hộp cắt lớn cải tiến chuyên dụng cho vật liệu tổng hợp địa kỹ thuật:
    • Mô hình cắt trượt mẫu đất trên mặt phẳng VĐKT dưới các cấp áp lực nén pháp tuyến $\sigma_n$.
    • Mô hình kéo rút cốt vải nằm giữa hai lớp đất đầm nén chuẩn để đo ứng suất tiếp $\tau$ và chuyển vị tương đối.
  3. Giao thức thử nghiệm mô hình vật lý 1:1: Đắp khối đê thử nghiệm theo từng lớp có trải cốt VĐKT phân tầng; lắp đặt hệ thống cọc mốc và cảm biến chuyển vị ngang $S_n$ tại các điểm trọng yếu (D2, D4). Đo đạc liên tục tốc độ tăng biến dạng ngang theo từng cấp tải trọng đắp sau 1 ngày và 3 ngày để đánh giá tốc độ tiêu tán áp lực nước lỗ rỗng và độ ổn định cố kết.

Data và phân tích

  • Đặc trưng mẫu và chỉ tiêu cơ lý: Khối lượng thể tích tự nhiên $\gamma = 16,5 \div 18,2\text{ kN/m}^3$; lực dính $c = 8 \div 16\text{ kPa}$; góc ma sát trong $\varphi = 6^\circ \div 14^\circ$; hệ số ma sát tương tác $f_{ds} = \tan\delta / \tan\varphi \approx 0,6 \div 0,8$; lực dính tiếp xúc $c^* \approx (0,8 \div 1,0)c$.
  • Kỹ thuật phân tích số: Lập trình giải thuật ma trận lặp phi tuyến tìm nghiệm hội tụ kép của vector lực tương tác $(\Delta E, \Delta X)$ và hệ số an toàn $F_s$.
  • Kiểm định Robustness Check: So sánh kết quả $F_s$ của phương pháp Luận án (PP Tương thích) với phương pháp Bishop đơn giản và phần mềm thương mại GeoStudio Slope/W qua ba kịch bản (TH1: $F_s = 1,0$; TH2: $F_s > 1,0$ không cốt; TH3: $F_s > 1,0$ có gia cố VĐKT). Kết quả cho thấy khi $F_s = 1,0$ (trạng thái cân bằng giới hạn chuẩn), phân bố hợp lực kháng trượt giữa hai phương pháp hoàn toàn trùng khít; khi $F_s \ne 1,0$, phương pháp tương thích phản ánh chính xác sự giảm ứng suất cắt đáy thỏi nhờ lực neo $V_i$, loại bỏ hiện tượng tập trung ứng suất giả tạo ở chân mái dốc.

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

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

  1. Giải quyết triệt để tính bất định của bài toán phân thỏi: Xác lập công thức giải tích giải $(2n - 2)$ ẩn số lực tương tác dựa trên điều kiện cân bằng nội tại của hệ thống đẩy - chống, chứng minh rằng hướng tác dụng của lực liên thỏi không phải là một hàm gán ngẫu nhiên $\lambda f(x)$ mà phụ thuộc trực tiếp vào độ dốc mái, cấp tải trọng ngoài và ma sát tương đối giữa các thỏi.
  2. Cơ chế tái phân bố ứng suất của cốt VĐKT: Khi đưa các lớp cốt VĐKT vào thân đê biển, lực kéo $V_i$ được huy động từng phần ($0% \to 50% \to 100%$) đã làm thay đổi căn bản quy luật phân bố số gia lực đẩy ngang $\Delta E_i \sin\alpha \tan\varphi$ và số gia lực cắt $\Delta X_i \cos\alpha \tan\varphi$. Hệ số an toàn ổn định tổng thể $F_s$ tăng tỷ lệ thuận phi tuyến với mức độ huy động lực kéo của cốt vải, làm chuyển cung trượt nguy hiểm nhất từ dạng trượt sâu trong thân đê ra sát mép ngoài mái đê.
  3. Hiệu ứng kìm hãm biến dạng ngang ($S_n$): Dữ liệu quan trắc từ mô hình vật lý 1:1 cho thấy tại điểm đo D4 và D2, khi có cốt VĐKT, tốc độ tăng biến dạng ngang sau 1 ngày và 3 ngày đắp giảm từ 40% đến 65% so với mô hình không cốt dưới cùng cấp tải trọng đắp. Điều này minh chứng cốt VĐKT phát huy tác dụng "khóa" biến dạng trượt hông của đất yếu ngay trong quá trình thi công.
  4. Tính khả thi vượt trội khi sử dụng vật liệu đất đắp tại chỗ: Ứng dụng thực tế tại đoạn đê biển KCT2 Bình Minh 3 (Kim Sơn, Ninh Bình) chứng minh công nghệ đất có cốt VĐKT cho phép tận dụng triệt để bùn cát pha, á sét cục bộ có chất lượng cơ lý kém ($c, \varphi$ thấp) để đắp đê với mái dốc thiết kế đứng hơn ($m = 1,5 \div 2,0$ so với truyền thống $m = 3,0 \div 4,5$), giúp giảm hơn 30% khối lượng đất đào đắp và diện tích chiếm đất vĩnh viễn.
  5. Khả năng chịu đựng thiên tai thực tế: Đoạn đê thử nghiệm KCT2 tại Ninh Bình sau khi thi công (tháng 8/2005) đã chống chịu an toàn tuyệt đối trước cơn bão số 7 lịch sử năm 2005 với triều cường và sóng biển dâng cao tràn qua đê mà không hề bị sạt lở hay trượt thân đê, tồn tại ổn định kéo dài suốt quá trình quan trắc đến năm 2009.
                  SO SÁNH CƠ CHẾ MÁI ĐÊ TRUYỀN THỐNG VÀ ĐÊ CÓ CỐT VĐKT

Implications đa chiều

  • Tiến bộ lý thuyết: Làm phong phú thêm cơ học đất tính toán (Computational Soil Mechanics), cung cấp một công cụ giải tích tất định giải bài toán siêu tĩnh phân thỏi, chấm dứt sự phụ thuộc vào các giả thiết kinh nghiệm chủ quan về góc lực liên thỏi.
  • Đổi mới phương pháp luận: Thiết lập quy trình chuẩn tích hợp giữa mô hình toán học giải tích, thí nghiệm xác định thông số tương tác giao diện ($f_{ds}, f_{po}$), và kiểm soát biến dạng hiện trường tỷ lệ 1:1.
  • Ứng dụng thực tiễn kỹ thuật: Cung cấp giải pháp kỹ thuật mang tính đột phá cho các vùng ven biển khan hiếm đất đắp đạt chuẩn (vùng Đồng bằng sông Cửu Long, đồng bằng sông Hồng), cho phép rút ngắn tiến độ đắp đê nhờ khả năng thoát nước lỗ rỗng nhanh qua các lớp VĐKT đóng vai trò tầng lọc và tiêu nước nội bộ.
  • Khuyến nghị chính sách đê điều: Định hình lại Tiêu chuẩn thiết kế đê biển Việt Nam trong bối cảnh thích ứng biến đổi khí hậu; chuyển đổi tư duy từ "đắp đê khổng lồ tiêu tốn đất" sang "kết cấu đê mềm gia cường công nghệ cao", tối ưu hóa nguồn vốn ngân sách đầu tư công trình thủy lợi quốc gia.

Limitations và Future Research

  1. Giới hạn hình học mặt trượt: Mô hình giải tích trong luận án hiện tập trung chủ yếu vào mặt trượt trụ tròn (Circular Slip Surface). Trong thực tế nền đất có phân tầng phức tạp hoặc lớp đệm cát mỏng xen kẹp bùn hữu cơ, mặt trượt thực tế có thể phát triển thành mặt trượt gãy khúc, mặt trượt hỗn hợp hoặc trượt phẳng dọc theo mặt tiếp xúc cốt vải.
  2. Ảnh hưởng của từ biến (Creep) và độ bền dài hạn của polyme: Luận án giả định cốt VĐKT làm việc trong giới hạn đàn hồi tuyến tính với tuổi thọ thiết kế công trình bán vĩnh cửu $30 \div 80\text{ năm}$. Tuy nhiên, hiện tượng suy giảm cường độ kéo do từ biến dưới tác dụng của tải trọng tĩnh duy trì liên tục và tác động ăn mòn vi sinh, oxy hóa trong môi trường nước biển mặn phèn chưa được mô hình hóa thành hàm thời gian động $V_i(t)$.
  3. Tác động tương hỗ động lực học sóng - đê: Nghiên cứu xét tải trọng tĩnh và gia tải tĩnh tương đương; các hiệu ứng thủy động lực học phức tạp như áp lực sóng va đập lặp (Cyclic Wave Loading) gây hóa lỏng cục bộ đất cát bão hòa nước trên mái đê cần tiếp tục được lượng hóa chi tiết.
  4. Chương trình nghiên cứu tương lai:
    • Mở rộng thuật toán tương thích lực tương tác cho mặt trượt có hình dạng bất kỳ (Non-circular Slip Surfaces).
    • Phát triển mô hình liên kết phần tử hữu hạn phi tuyến (FEM) ghép nối cơ chế thủy lực - biến dạng (Hydro-mechanical Coupled Modeling) phân tích quá trình thấm bão hòa và biến dạng cố kết đồng thời.
    • Ứng dụng công nghệ sợi quang thông minh (Fiber Optic Sensors) tích hợp trực tiếp vào cấu trúc VĐKT để quan trắc thời gian thực biến dạng và nội lực kéo $V_i$ trong thân đê biển phục vụ cảnh báo sớm sạt trượt.

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

  • Ảnh hưởng học thuật: Luận án mở ra một nhánh tiếp cận toán - cơ địa kỹ thuật mới trong nước và khu vực; là tài liệu tham khảo then chốt trong các chương trình đào tạo sau đại học chuyên ngành Công trình Thủy, Địa kỹ thuật và Xây dựng Công trình Thủy lợi.
  • Tác động ngành và chuyển giao công nghệ: Công nghệ đất có cốt VĐKT được ứng dụng rộng rãi trong các dự án nâng cấp tuyến đê biển từ Quảng Ninh đến Quảng Nam, các dự án đê bao bờ bao ngăn mặn thích ứng biến đổi khí hậu tại Đồng bằng sông Cửu Long, mái dốc đường cao tốc và taluy chống sạt lở miền núi.
  • Hiệu quả kinh tế - xã hội: Giảm chi phí xây dựng từ 20% đến 35% so với giải pháp truyền thống nhờ giảm cự ly vận chuyển đất đắp xa xôi, tận dụng hàng triệu mét khối đất bùn cát tại chỗ, đồng thời bảo vệ an toàn tính mạng và cơ sở hạ tầng kinh tế ven biển của hàng triệu hộ dân trước các hiểm họa bão lũ cực đoan.
                    MA TRẬN TÁC ĐỘNG VÀ ĐỐI TƯỢNG HƯỞNG LỢI

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận giải tích hệ thống chuẩn xác để xử lý các bài toán siêu tĩnh địa kỹ thuật; tham khảo quy trình thực nghiệm từ phòng đến hiện trường tỷ lệ thật 1:1.
  • Các nhà khoa học và chuyên gia Cơ học đất: Sở hữu một mô hình tính toán tương thích lực tương tác có độ tin cậy cao, giải quyết các tranh luận học thuật kéo dài về giả thiết lực liên thỏi trong lý thuyết phân thỏi cổ điển.
  • Kỹ sư tư vấn thiết kế và doanh nghiệp xây dựng thủy lợi: Ứng dụng trực tiếp thuật toán và chương trình phần mềm vào tính toán thiết kế đê biển, đê sông, đập đất, tường chắn có cốt; tối ưu hóa mặt cắt ngang công trình, rút ngắn tiến độ và kiểm soát chất lượng đầm nén vật liệu địa phương.
  • Bộ Nông nghiệp & PTNT, Chi cục Quản lý Đê điều & PCTT: Có cơ sở khoa học - thực tiễn vững chắc để phê duyệt các dự án đầu tư nâng cấp đê biển quốc gia thích ứng với kịch bản nước biển dâng và biến đổi khí hậu toàn cầu.

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 lý thuyết độc đáo nhất là việc thiết lập điều kiện tương thích cho hệ lực tương tác giữa các thỏi đất ($\Delta E_i, \Delta X_i$), trực tiếp mở rộng lý thuyết cân bằng giới hạn phân thỏi cổ điển và lý thuyết phân tích hệ thống trong địa kỹ thuật của GS. Phan Trường Phiệt. Thay vì giải bài toán siêu tĩnh bằng cách "áp đặt các giả thiết nhân tạo" (như Bishop bỏ qua $X_i=0$, Janbu giả định đường đẩy $h_p$, Spencer/GLE giả định hàm góc $\lambda f(x)$), luận án đã ứng dụng định lý giới hạn Gvozdev và nguyên lý cực trị Coulomb để xác lập trạng thái cân bằng nội tại giữa khối đẩy và khối cản trượt qua từng lát cắt đứng. Khi có sự hiện diện của cốt VĐKT, véc-tơ lực kéo neo $V_i$ được tích hợp trực tiếp vào phương trình tương thích vi phân của thỏi đất, phản ánh đúng bản chất vật lý của sự phân bố lại ứng suất tiếp và ứng suất pháp mà không cần dựa vào bất kỳ tham số hình học giả định nào.

2. Sự đổi mới về phương pháp luận nghiên cứu so với các công trình trước đây thể hiện ở điểm nào?

So với các nghiên cứu trước đây (như nghiên cứu của Spencer 1967, Janbu 1973, hay Schmertmann 1987 vốn chỉ coi lực cốt vải như một mômen giữ thuần túy tác dụng lên toàn khối trượt $M_{V/0} = V_i \cdot R$), phương pháp luận của luận án mang tính đột phá ở quy trình kiểm chứng 4 cấp khép kín:

  1. Cấp độ giải tích vi phân: Xây dựng mô hình toán tương thích lực liên thỏi đa ẩn số.
  2. Cấp độ thực nghiệm tiếp xúc: Sử dụng hộp cắt phẳng cải tiến và kéo tuột xác định chính xác các hệ số tương tác thực $f_{ds}, f_{po}$ thay vì lấy hệ số quy đổi giả định.
  3. Cấp độ mô hình vật lý quy mô 1:1: Đo đạc trường biến dạng ngang thực tế $S_n$ theo thời gian tại các điểm D2, D4 để kiểm chứng động lực học cố kết và tiêu tán áp lực nước lỗ rỗng.
  4. Cấp độ số hóa và kiểm chứng chéo: Tự phát triển phần mềm tính toán chuyên dụng và đối chiếu sai số trực tiếp với chuẩn thương mại quốc tế GeoStudio (Slope/W) trên các kịch bản chuẩn $F_s=1$ và $F_s>1$.

3. Phát hiện nào trong luận án mang tính bất ngờ nhất và được chứng minh bằng dữ liệu thực nghiệm ra sao?

Phát hiện bất ngờ và có ý nghĩa kỹ thuật sâu sắc nhất là: Lực căng cốt vải $V_i$ không làm tăng đồng đều hệ số an toàn trên toàn bộ mặt trượt, mà thực hiện tái cấu trúc trường ứng suất trượt cục bộ, đẩy cung trượt nguy hiểm nhất ($F_s^{\min}$) từ vùng trượt sâu trong thân đê ra sát mép biên mái dốc. Dữ liệu thực nghiệm chứng minh:

  • Khi mức độ huy động lực kéo của VĐKT tăng từ 0% lên 100%, giá trị $F_s$ của mái phía đồng và mái phía biển tăng từ dưới ngưỡng phá hoại ($F_s < 1,0$) lên mức an toàn cao ($F_s = 1,35 \div 1,58$).
  • Phân bố số gia lực kháng trượt $\Delta X_i \cos\alpha \tan\varphi$ và $\Delta E_i \sin\alpha \tan\varphi$ tại các thỏi cắt qua cốt vải triệt tiêu ứng suất cắt tập trung đáy thỏi, làm phẳng biểu đồ ứng suất trượt dọc cung trượt.
  • Mô hình vật lý 1:1 xác nhận biến dạng ngang $S_n$ tại điểm D4 sau 3 ngày đắp tải ở mô hình có cốt giảm hơn một nửa so với mô hình không cốt, ngăn chặn hoàn toàn nguy cơ hình thành mặt trượt cắt sâu liên tục.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp hệ thống giao thức tái lập hoàn chỉnh và minh bạch:

  • Toàn bộ hệ phương trình toán học giải tích xác định $\Delta E, \Delta X$, các ma trận cân bằng tĩnh học và sơ đồ khối thuật toán lặp tìm $F_s$ được trình bày chi tiết trong Chương II và Chương III.
  • Các biểu mẫu thông số đầu vào cơ lý đất ($c, \varphi, \gamma, w$), quy cách vật liệu VĐKT (độ bền kéo, độ giãn dài, khối lượng đơn vị $g/m^2$, kích thước lỗ mở lọc), kích thước hình học mặt cắt đê thử nghiệm KCT2 Bình Minh 3 và vị trí bố trí các lớp cốt (khoảng cách đứng $\Delta h$, chiều dài neo $L_{\text{neo}}$) được công bố cụ thể trong Chương IV và hệ thống Phụ lục (Phụ lục 1 đến Phụ lục 4). Bất kỳ nhà nghiên cứu nào cũng có thể lập trình lại giải thuật hoặc mô phỏng lại trên các phần mềm địa kỹ thuật tiêu chuẩn.

5. Luận án đã phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào?

Luận án định hình chương trình nghiên cứu dài hạn tập trung vào 4 trọng tâm:

  1. Mở rộng trường lý thuyết: Hoàn thiện mô hình tương thích lực cho cung trượt gãy khúc bất kỳ và không gian ba chiều (3D Slope Stability Analysis).
  2. Độ bền dài hạn và tương tác môi trường biển: Nghiên cứu thực nghiệm từ biến dài hạn của sợi polyme dệt và không dệt trong môi trường xâm thực mặn, phèn và chu kỳ khô ướt triều kiệt.
  3. Động lực học thủy văn - địa kỹ thuật: Tích hợp bài toán thấm không dừng khi triều rút nhanh kết hợp sóng tràn đỉnh đê trong các trận bão cấp 12–15.
  4. Thông minh hóa đê biển: Phát triển vật liệu địa kỹ thuật thông minh tích hợp cảm biến đo biến dạng tự động phục vụ công tác giám sát trực tuyến an toàn đê điều quốc gia.

Kết luận

  1. Hoàn thiện lý thuyết phân thỏi địa kỹ thuật: Luận án đã giải quyết thành công bài toán siêu tĩnh bậc cao $(6n - 2)$ ẩn số trong phân tích ổn định mái dốc bằng phương pháp phân thỏi thông qua việc xác lập điều kiện tương thích lực tương tác $(\Delta E, \Delta X)$, loại bỏ các giả thiết chủ quan, thiếu cơ sở cơ học của các phương pháp truyền thống.
  2. Làm sáng tỏ cơ chế làm việc của đất có cốt VĐKT: Mô hình hóa tường minh tác động của lực kéo neo $V_i$ trong việc phân bố lại trường ứng suất nội tại, chứng minh cơ chế nâng cao hệ số ổn định $F_s$ và kìm hãm chuyển vị ngang $S_n$ của thân đê trên nền đất yếu.
  3. Phát triển công cụ tính toán số tin cậy: Xây dựng phần mềm chuyên dụng tính toán ổn định mái dốc có xét tương thích lực tương tác, được kiểm chứng đối chuẩn chính xác với phần mềm quốc tế GeoStudio Slope/W.
  4. Đột phá công nghệ xây dựng đê biển: Mở ra giải pháp kỹ thuật cho phép tận dụng tối đa đất yếu cục bộ tại chỗ (cát pha, bùn sét $0,005 \div 0,5\text{ mm}$) để đắp đê biển bán vĩnh cửu có mái dốc đứng hơn, tiết kiệm trên 30% khối lượng đào đắp và diện tích chiếm đất.
  5. Kiểm chứng thực nghiệm toàn diện: Thành công của đoạn đê thử nghiệm KCT2 Bình Minh 3 (Kim Sơn, Ninh Bình) đứng vững an toàn qua bão số 7 (2005) và nhiều năm vận hành khẳng định tính chuẩn xác, độ tin cậy và giá trị thực tiễn to lớn của công trình.
  6. Giá trị kế thừa và tầm ảnh hưởng lâu dài: Đóng góp của luận án là bước tiến quan trọng trong cơ học đất và công trình thủy tại Việt Nam, cung cấp luận cứ khoa học vững chắc phục vụ chiến lược nâng cấp toàn diện 2.700 km đê biển quốc gia chủ động ứng phó với biến đổi khí hậu và nước biển dâng trong kỷ nguyên mới.