Tổng quan về luận án

Nghiên cứu cơ học đất và kỹ thuật địa kỹ thuật công trình giao thông tại Việt Nam trong giai đoạn công nghiệp hóa, hiện đại hóa đặt ra yêu cầu cấp bách về việc tối ưu hóa giải pháp gia cố nền đắp trên đất yếu và nền đắp cao. Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu ổn định nền đường đắp có gia cường bằng vải địa kỹ thuật trong các công trình xây dựng đường ô tô ở Việt Nam" do Nghiên cứu sinh Huỳnh Ngọc Hào thực hiện (chuyên ngành Xây dựng đường ô tô và đường thành phố, Mã số: 62.01, dưới sự hướng dẫn khoa học của GS.TS. Vũ Đình Phụng và TS. Vũ Đức Sỹ, Trường Đại học Giao thông Vận tải, 2014) là một công trình tiên phong giải quyết toàn diện bài toán tương tác phức tạp giữa đất và cốt mềm polyme tổng hợp.

Thực tiễn xây dựng hạ tầng đường ô tô tại Việt Nam ghi nhận sự mở rộng quy mô ứng dụng vải địa kỹ thuật (VĐKT - Geotextiles) từ dự án cải tạo Quốc lộ 5 năm 1993 (với khối lượng hơn $500.000\text{ m}^2$) đến hàng loạt tuyến cao tốc huyết mạch như TP. Hồ Chí Minh – Trung Lương, Đại lộ Thăng Long, Cầu Giẽ – Ninh Bình, Nội Bài – Lào Cai. Tuy nhiên, một research gap (khoảng trống nghiên cứu) mang tính cốt lõi tồn tại suốt nhiều thập kỷ: các tiêu chuẩn thiết kế hiện hành (22TCN 207-1992, 22TCN 219-1994, BS-6349, OCDI) và các phần mềm thương mại quốc tế (Plaxis 2D, Phase2, GeoStudio) phần lớn dựa trên phương pháp cân bằng giới hạn (Limit Equilibrium Method - LEM) giả thiết mặt trượt tròn (Terzaghi, Bishop, Janbu) hoặc xem quan hệ ứng xử kéo của VĐKT là đàn hồi dẻo lý tưởng tuyến tính theo tiêu chuẩn phá hoại Mohr-Coulomb. Thực tế biến dạng kéo của VĐKT và bề mặt tiếp xúc đất - cốt mang tính chất phi tuyến nhiều giai đoạn, đồng thời dạng mặt trượt thực tế trong nền đắp gia cường không tuân theo cung tròn cổ điển.

Luận án thiết lập hệ thống 04 câu hỏi nghiên cứu và giả thuyết khoa học tương ứng:

  • Câu hỏi 1 (RQ1): Làm thế nào để mô hình hóa chính xác ứng xử phi tuyến đa giai đoạn của VĐKT và phần tử tiếp xúc (interface) trong bài toán biến dạng phẳng bằng phương pháp số?
    • Giả thuyết 1 (H1): Tích hợp mô hình phi tuyến Robert M. Koerner vào phần tử thanh chịu kéo và phần tử tiếp xúc sẽ phản ánh chính xác trạng thái ứng suất - biến dạng thực tế của hệ "đất - cốt".
  • Câu hỏi 2 (RQ2): Hình dạng cung trượt nguy hiểm nhất của mái dốc nền đắp gia cường trên đất yếu thực tế tuân theo quy luật hình học nào?
    • Giả thuyết 2 (H2): Mặt trượt thực tế phát triển theo đường cong Ellipse thay vì cung tròn $Rd$ hoặc khối nêm hai phần gãy khúc.
  • Câu hỏi 3 (RQ3): Mức độ ảnh hưởng tương quan giữa độ cứng dọc trục ($EA_g$), cường độ chịu kéo cực đại ($T_{max}$) của VĐKT và mô đun đàn hồi đất nền ($E_s$) tới hệ số an toàn ($F_s$) diễn ra như thế nào?
    • Giả thuyết 3 (H3): Hệ số an toàn $F_s$ chịu sự chi phối phi tuyến của độ cứng $EA_g$; tồn tại ngưỡng giới hạn mà vượt qua đó, việc tăng $EA_g$ không làm tăng thêm $F_s$.
  • Câu hỏi 4 (RQ4): Có thể xây dựng một công cụ phần mềm chuyên dụng và hệ thống biểu đồ tra cứu chuẩn hóa cho kỹ sư giao thông Việt Nam hay không?
    • Giả thuyết 4 (H4): Thuật toán phần tử hữu hạn kết hợp phương pháp giảm cường độ kháng cắt $c-\phi$ cho phép phát triển phần mềm HNH_RESS V1.00 đạt độ chính xác tương đương các phần mềm quốc tế nhưng tối ưu hóa riêng cho điều kiện địa kỹ thuật Việt Nam.

Phạm vi nghiên cứu bao quát các chiều cao nền đắp điển hình $H = 6\text{ m}, 8\text{ m}, 10\text{ m}, 12\text{ m}$ với hệ số mái dốc $m = 1:1; 1:1,25; 1:1,5$, cường độ chịu kéo VĐKT biến thiên $T_{max} = 12 \div 28\text{ kN/m}$. Công trình tạo ra bước đột phá định lượng khi chứng minh việc ứng dụng VĐKT đúng kỹ thuật giúp "giảm giá thành từ 15 - 20%, tăng chất lượng sử dụng, tăng tuổi thọ của công trình" và kiểm soát hoàn toàn nguy cơ trượt trồi, mất ổn định tổng thể.

Literature Review và Positioning

Lịch sử phát triển của vải địa kỹ thuật ghi nhận sự ra đời từ thập niên 1960 tại Pháp dưới tên thương mại Bidim. Bước ngoặt ứng dụng quy mô lớn diễn ra từ năm 1965 đến 1973 tại Hà Lan, khi công trình đê biển Delta Works sử dụng từ $5 \div 6\text{ triệu m}^2$ VĐKT nhằm ngăn mặn và gia cố đê. Tại Đông Nam Á, VĐKT được triển khai thành công tại các dự án mái dốc Bukit Panggal Mosque (Tutong, Brunei), tường chắn Arca Budaya (Kuala Lumpur, Malaysia), tường chắn khách sạn Novotel Patong (Phuket, Thái Lan), và đường đê Punarbhava (Tây Bengal, Ấn Độ).

Tổng quan y văn xác định hai trường phái tính toán và tranh luận học thuật đối lập:

  1. Trường phái Cân bằng Giới hạn Cổ điển (LEM): Đại diện bởi K. Terzaghi (1943), A.W. Bishop (1955), N. Janbu (1956) với giả thiết mặt trượt trụ tròn phân mảnh; Leschinsky & Boedecker (1989), Bridle & Barr (1990) phát triển mặt trượt xoắn ốc logarit; Segrestin & Bastick (1988) với phương pháp trọng lực dính kết (Rankine). Nhược điểm chí mạng của trường phái này là giả định khối trượt là vật thể cứng tuyệt đối, bỏ qua hoàn toàn mô đun đàn hồi $E_s, E_g$ và độ cứng tương đối $EA_g$, dẫn đến việc xác định lực căng yêu cầu $T_{roj}$ và chiều dài neo bám $L_{ej}$ mang tính gần đúng thô sơ: $$L_{ej} \ge \frac{f_n \cdot f_p \cdot T_{roj}}{2 \left( \frac{f_{fs} \gamma h + f_q w_s}{f_{ms}} \cdot \alpha' \tan\varphi'p + \frac{\alpha'{bc} c'}{f_{ms}} \right)}$$
  2. Trường phái Phương pháp Phần tử Hữu hạn (FEM) và Cơ học Môi trường Liên tục: Được ứng dụng trong Plaxis 2D (Hà Lan) và Phase2 (Canada, Rocscience). Tuy nhiên, các phần mềm này "trong tính toán đều xem quan hệ ứng xử kéo của VĐKT và phần tử tiếp xúc VĐKT với đất nền là đàn hồi dẻo lý tưởng, quan hệ tuyến tính theo tiêu chuẩn phá hoại Mohr-Coulomb".

Vị thế nghiên cứu (positioning) của luận án được xác lập rõ nét: bắc cầu nối giữa cơ học tiếp xúc phi tuyến thực nghiệm của Robert M. Koerner (2005) với giải thuật phần tử hữu hạn 2D, thay thế mặt trượt trụ tròn truyền thống bằng mô hình mặt trượt Ellipse chính xác hóa, đồng thời tạo ra bộ công cụ tính toán tự chủ đóng gói hoàn chỉnh mang tên HNH_RESS.

Đó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à thách thức các giới hạn của lý thuyết cân bằng giới hạn cổ điển và lý thuyết đàn hồi dẻo tuyến tính áp dụng cho vật liệu tổng hợp địa kỹ thuật:

  • Thách thức mô hình Mohr-Coulomb tuyến tính đối với VĐKT: Luận án chứng minh rằng biến dạng của VĐKT polyme (dệt sợi đơn, dệt dạng tấm, không dệt xuyên kim, không dệt nhiệt hóa) tuân theo quan hệ ứng suất - biến dạng ($\sigma - \varepsilon$) phi tuyến trải qua 3 pha rõ rệt: $$\sigma = f(\varepsilon) \implies E_g = E_{initial} \to E_{secant} (E_{S10}, E_{S35})$$
  • Mở rộng lý thuyết tương tác tiếp xúc đất - cốt (Soil-Geosynthetic Interface Mechanics): Khẳng định cơ chế trượt tiếp xúc giữa đất và VĐKT gồm 3 pha liên tục theo mô hình Koerner: pha phi tuyến ban đầu (đoạn 0-1, ứng suất tiếp tăng nhanh, biến dạng trượt nhỏ), pha tái bền (strain hardening, đoạn 1-2, ứng suất và biến dạng tăng đồng thời), và pha hóa mềm (strain softening, đoạn 2-3, ứng suất trượt suy giảm dần về giá trị dư khi biến dạng tiếp tục tăng).
  • Xác lập mô hình mặt trượt Ellipse thay thế mặt trượt tròn: Bằng cách phân tích trường chuyển vị nút và vector biến dạng trượt lớn nhất từ PTHH, luận án chứng minh hình dạng mặt trượt thực tế của nền đắp gia cường VĐKT trên đất yếu là đường Ellipse có phương trình tham số: $$\frac{(x - x_0)^2}{a^2} + \frac{(y - y_0)^2}{b^2} = 1$$ Từ đó, xây dựng công thức giải tích xác định phân bố lực căng cực đại $T_{max}$ dọc theo cung trượt Ellipse cho từng lớp cốt $j$.

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

Khung phân tích của đề tài tích hợp chặt chẽ 03 trụ cột lý thuyết: Lý thuyết biến dạng phẳng đàn hồi - dẻo môi trường liên tục (Định luật Hooke mở rộng), Tiêu chuẩn chảy dẻo Mohr-Coulomb cho vật liệu đất đắp/đất nền, và Cơ học vật liệu phân tử Polyme phi tuyến Robert M. Koerner (2005).

Điều kiện biên và phạm vi áp dụng (Boundary Conditions):

  • Trạng thái biến dạng phẳng: $\varepsilon_z = 0, \gamma_{yz} = 0, \gamma_{zx} = 0$.
  • Cốt VĐKT phân bố nằm ngang, chỉ chịu lực dọc trục kéo ($N \ge 0$), triệt tiêu khả năng chịu nén ($N < 0 \implies [K_g] = 0$) và độ cứng uốn bỏ qua ($EJ = 0$).
  • Giới hạn bền kéo đơn vị xác định theo độ dày danh định $t = 0,25 \div 3,5\text{ mm}$ ở áp suất chuẩn $2\text{ kPa}$ (ASTM D5199).

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường thực chứng (positivism) kết hợp thuyết hiện thực phản biện (critical realism), sử dụng phương pháp giải tích số phần tử hữu hạn phi tuyến 2D làm trọng tâm, đối sánh trực tiếp với nghiệm giải tích giải tích cổ điển và kiểm chứng chéo (cross-verification) với phần mềm thương mại Plaxis 2D (Hà Lan).

Mô hình hóa hình học bài toán gồm:

  1. Nền đất (Soil mass): Rời rạc hóa bằng phần tử tam giác 3 nút và phần tử tấm tam giác đẳng tham số 6 nút với ma trận chuyển đổi tọa độ tự nhiên $(\xi, \eta)$.
  2. Cốt VĐKT (Reinforcement): Mô phỏng bằng phần tử thanh chịu lực dọc trục 2 nút (truss element), hàm dạng tuyến tính bậc nhất.
  3. Lớp tiếp xúc (Interface element): Phần tử tiếp xúc 4 nút vô cùng mỏng mô tả trượt tương đối giữa các nút đất và nút cốt liền kề.

Quy trình nghiên cứu rigorous

Quy trình tính toán tuân thủ thuật toán lặp số giải hệ phương trình phi tuyến đại số giải tích ma trận độ cứng: $$[K(u)] {u} = {P_{ext}}$$

Để tìm hệ số an toàn ổn định $F_s$, luận án áp dụng nguyên lý suy giảm độ bền chống cắt (c - $\varphi$ reduction technique): $$\frac{\tan\varphi_i}{\tan\varphi_r} = \frac{c_i}{c_r} = \frac{s_{u,i}}{s_{u,r}} = F_s \implies \varphi_r = \arctan\left(\frac{\tan\varphi_i}{F_s}\right); \quad c_r = \frac{c_i}{F_s}$$

Hệ phương trình phi tuyến tại từng bước tăng tải trọng và suy giảm độ bền được giải bằng thuật toán Newton-Raphson tiêu chuẩn và Newton-Raphson cải tiến (Modified Newton-Raphson - MNR) nhằm tối ưu hóa tốc độ hội tụ và triệt tiêu dao động số khi đất bước vào trạng thái dẻo chảy.

                    Hội tụ              Không hội tụ (Sụp đổ dẻo)
                 Tăng tiếp Fs              Fs_limit = Fs_critical

Data và phân tích

Toàn bộ thuật toán được tác giả độc lập lập trình và đóng gói thành phần mềm chuyên dụng HNH_RESS V1.00 (Reinforced Embankment Stability Software). Ma trận kiểm tra thực nghiệm số bao gồm hàng trăm kịch bản tính toán chi tiết:

  • Chiều cao nền đắp: $H = 6\text{ m}, 8\text{ m}, 10\text{ m}, 12\text{ m}$.
  • Cường độ chịu kéo của VĐKT: $T_{max} = 12, 14, 16, 18, 20, 22, 24, 26, 28\text{ kN/m}$.
  • Khoảng cách bố trí lớp cốt: $S_v = 0,5\text{ m}$ và $1,5\text{ m}$.
  • Số lượng lớp vải gia cường: $n = 1, 2, 3, 4\text{ lớp}$.
  • Tải trọng xe cộ thiết kế: Mô hình hóa đoàn xe tiêu chuẩn quy đổi thành tải trọng phân bố tương đương $q$.

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

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

  1. Phát hiện 1: Dạng hình học mặt trượt thực tế là Cung Ellipse thay vì Cung tròn. Phân tích trường biến dạng trượt trong HNH_RESS chỉ ra rằng, đối với nền đắp gia cường VĐKT trên đất yếu, sự hiện diện của các lớp cốt mềm có tính cản trượt và truyền ứng suất cắt sang phương ngang đã kéo giãn vùng phá hoại dẻo. Mặt trượt nguy hiểm nhất không bao giờ là cung tròn $Rd$ mà biến dạng thành dạng Ellipse lệch tâm rõ rệt. Sai số góc trượt giữa cung tròn giả thiết và cung Ellipse thực tế có thể dẫn đến việc đánh giá sai hệ số an toàn từ $12% \div 18%$.

  2. Phát hiện 2: Quy luật phân bố lực căng lớn nhất $T_{max}$ dọc theo thân cốt. Lực căng trong các lớp VĐKT không phân bố đều mà đạt cực đại tại giao điểm giữa lớp vải với mặt trượt Ellipse nguy hiểm. Công thức giải tích mới do tác giả thiết lập phản ánh chính xác giá trị $T_{max}$ huy động tại từng cao độ lớp cốt $j$, khắc phục hoàn toàn nhược điểm tính đồng đều $T_{roj}$ của phương pháp phân mảnh cổ điển.

  3. Phát hiện 3: Tác động phi tuyến của độ cứng $EA_g$ và hiện tượng "bão hòa hiệu quả gia cường". Khi tăng độ cứng kéo $EA_g$ của VĐKT từ $100\text{ kN/m}$ lên $1.000\text{ kN/m}$, hệ số an toàn $F_s$ của mái dốc tăng trưởng rõ rệt (từ $15% \div 28%$). Tuy nhiên, khi $EA_g > 1.500\text{ kN/m}$, tốc độ gia tăng của $F_s$ tiệm cận về 0 (đạt ngưỡng bão hòa). Điều này chứng minh rằng việc tăng độ cứng quá mức không đem lại hiệu quả kinh tế - kỹ thuật do sự khống chế của cường độ chống cắt nội tại của đất đắp ($c, \varphi$).

  4. Phát hiện 4: Phân tích so sánh giữa nguy cơ đứt cốt ($T > T_{max}$) và tuột cốt ($L_e < L_{ej}$). Đối với các nền đắp cao ($H \ge 10\text{ m}$), sự phá hoại chủ yếu bị khống chế bởi hiện tượng đứt cốt tại chân mái taluy; trong khi đối với nền đắp thấp trên đất yếu ($H \le 6\text{ m}$), cơ chế phá hoại chủ đạo lại là tuột neo bám ở đuôi cốt do ứng suất pháp hiệu quả $\sigma'_v$ không đủ lớn để tạo ma sát bám dính.

  5. Phát hiện 5: Kiểm chứng tương quan giữa HNH_RESS V1.00 và Plaxis 2D. Kết quả so sánh hệ số an toàn ổn định cho nền đắp cao $6\text{ m}$ và $8\text{ m}$ giữa chương trình HNH_RESS và Plaxis cho thấy độ sai lệch cực kỳ nhỏ, chỉ dao động trong khoảng $1,2% \div 3,5%$. Điều này khẳng định độ tin cậy tuyệt đối của giải thuật phần tử hữu hạn và mã nguồn do tác giả tự xây dựng.

                    BẢNG SO SÁNH NĂNG LỰC TÍNH TOÁN

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp khung cơ học hoàn chỉnh mô tả hệ tương tác "Đất - Kết cấu cốt mềm Polyme", giải quyết thỏa đáng mối liên hệ giữa biến dạng trượt cục bộ và sự phát triển lực kéo trong cốt.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình phân tích ổn định mái dốc bằng phần tử hữu hạn kết hợp thuật toán giảm $c-\phi$ cho các kết cấu địa kỹ thuật có cốt mềm tại Việt Nam.
  • Về mặt thực tiễn và thiết kế: Xây dựng hệ thống biểu đồ tra cứu quan hệ giữa chiều cao đắp $H$, độ dốc mái $m$, số lớp VĐKT $n$, cường độ $T_{max}$ và mô đun đất $E_s$. Kỹ sư thiết kế có thể tra cứu trực tiếp thông số vải địa mà không cần chạy lại các mô hình PTHH phức tạp trong giai đoạn lập báo cáo nghiên cứu khả thi (FS).
  • Về mặt chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học thực chứng để sửa đổi, bổ sung Tiêu chuẩn Quốc gia TCVN về thiết kế nền đường có gia cường vải địa kỹ thuật, chuyển dịch từ tiếp cận hệ số an toàn tổng thể sang tiếp cận trạng thái giới hạn (Limit States Design).

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận 04 giới hạn nghiên cứu:

  1. Mô hình bài toán không gian 2D: Nghiên cứu tập trung vào bài toán biến dạng phẳng 2D, chưa xét đến hiệu ứng không gian 3 chiều tại các vị trí đầu cầu, mố cầu hoặc đoạn chuyển tiếp hình học phức tạp.
  2. Yếu tố từ biến theo thời gian (Creep behavior): Đặc tính suy giảm cường độ kéo và dão từ biến của polyme dưới tác dụng của tải trọng duy trì dài hạn chưa được tích hợp trực tiếp vào hàm quan hệ ứng suất - biến dạng tức thời.
  3. Áp lực nước lỗ rỗng động (Dynamic excess pore pressure): Mô hình tính toán chủ yếu xét tải trọng xe cộ quy đổi tĩnh tương đương, chưa mô phỏng trực tiếp tải trọng động trùng phục (cyclic loading) gây hóa lỏng hoặc gia tăng áp lực nước lỗ rỗng tức thời trong nền đất yếu bão hòa nước.
  4. Vấn đề ăn mòn sinh hóa: Chưa định lượng ảnh hưởng của môi trường hóa chất, pH đất và vi sinh vật đến độ bền kéo dài hạn của VĐKT.

Chương trình nghiên cứu tương lai (Future Research Directions):

  • Phát triển mô hình phần tử hữu hạn 3D phân tích ổn định nền đắp gia cố VĐKT tại các vị trí tiếp giáp mố cầu.
  • Tích hợp mô hình lưu biến dẻo nhiệt (visco-elasto-plastic) để dự báo chính xác hiện tượng từ biến của vải địa sau 20-50 năm khai thác.
  • Phân tích tương tác động lực học đất - VĐKT dưới tác động của tải trọng xe tải nặng và động đất bằng phương pháp sai phân ngầm phi tuyến.
  • Mở rộng thuật toán HNH_RESS sang các dạng vật liệu địa kỹ thuật mới như lưới địa kỹ thuật (Geogrids), ô địa kỹ thuật (Geocells), và bấc thấm thẳng đứng (PVD).

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

Luận án mang lại những tác động sâu rộng trên cả phương diện học thuật và kinh tế - xã hội:

  • Tác động học thuật (Academic Impact): Tạo tiền đề cho hàng loạt nghiên cứu chuyên sâu về địa kỹ thuật số tại các trường đại học khối công trình (Trường ĐH Giao thông Vận tải, ĐH Xây dựng, ĐH Thủy lợi, Học viện Kỹ thuật Quân sự). Đóng góp các công thức giải tích mặt trượt Ellipse có tiềm năng trích dẫn cao trong các tạp chí chuyên ngành Geotextiles and Geomembranes, Geosynthetics International.
  • Chuyển đổi ngành công nghiệp xây dựng giao thông: Ứng dụng kết quả tính toán vào thiết kế thực tế giúp các nhà thầu và tư vấn (như TEDI, VEC) tối ưu hóa số lớp vải địa, giảm chiều sâu đào bóc thay đất, rút ngắn tiến độ thi công từ $20% \div 30%$ và giảm thiểu giá thành xây lắp từ $15% \div 20%$.
  • Lợi ích xã hội và môi trường: Hạn chế việc khai thác vật liệu đất đắp chọn lọc quý hiếm, giảm khối lượng đất đào đổ thải gây ô nhiễm môi trường, đảm bảo an toàn giao thông và nâng cao tuổi thọ khai thác của hệ thống hạ tầng đường bộ quốc gia.

Đố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 giải thuật mã nguồn, phương pháp lập trình phần tử hữu hạn phi tuyến cho kết cấu tương tác đất - cốt, mở rộng các hướng nghiên cứu vật liệu mới.
  • Kỹ sư tư vấn thiết kế cầu đường: Sử dụng hệ thống biểu đồ nomograph tiện ích để chọn nhanh loại VĐKT, xác định chính xác số lớp, khoảng cách rải và chiều dài neo bám an toàn.
  • Ban quản lý dự án và Cơ quan quản lý nhà nước (Bộ GTVT, Cục Đường bộ): Có cơ sở khoa học để thẩm định phương án kinh tế - kỹ thuật, phê duyệt dự toán tối ưu hóa chi phí xử lý nền đất yếu.
  • Doanh nghiệp sản xuất vật liệu Geosynthetics: Định hướng sản xuất các dòng sản phẩm VĐKT dệt và không dệt có các chỉ tiêu cơ lý ($EA_g, T_{max}$, độ giãn dài $\varepsilon%$) phù hợp nhất với điều kiện địa chất công trình Việt Nam.

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 mô hình cơ học tiếp xúc phi tuyến đa giai đoạn cho hệ "đất - VĐKT" dựa trên mô hình Robert M. Koerner (2005) và chứng minh quy luật mặt trượt phá hoại dạng Ellipse. Luận án đã vượt qua giới hạn của lý thuyết Mohr-Coulomb đàn hồi dẻo lý tưởng tuyến tính (vốn được đóng gói cố định trong Plaxis và Phase2), phản ánh trung thực bản chất vật liệu dẻo gia cường có biến dạng dẻo hóa mềm và tái bền.

2. Đổi mới phương pháp luận so với các nghiên cứu tiền nhiệm?

So với phương pháp giải tích phân mảnh cân bằng giới hạn truyền thống (Terzaghi, Bishop, Janbu) vốn giả định mặt trượt trụ tròn cứng tuyệt đối và bỏ qua độ cứng $EA_g$, luận án đã:

  • Phát triển thành công thuật toán PTHH 2D kết hợp phương pháp giảm cường độ kháng cắt $c-\phi$.
  • Mô hình hóa thành công phần tử thanh chịu kéo 1D phi tuyến không chịu nén và phần tử tiếp xúc 4 nút đặc thù.
  • Đóng gói hoàn chỉnh thành phần mềm HNH_RESS V1.00, cho phép xác định duy nhất một hệ số an toàn $F_s$ hội tụ mà không cần giả định trước vị trí tâm trượt.

3. Phát hiện nào gây bất ngờ nhất về mặt số liệu thực nghiệm?

Phát hiện bất ngờ nhất là sự tồn tại của ngưỡng bão hòa độ cứng $EA_g \approx 1.500\text{ kN/m}$. Trong thực tế kỹ thuật, các kỹ sư thường kỳ vọng việc tăng độ cứng vải địa kỹ thuật sẽ luôn làm tăng tuyến tính hệ số an toàn mái dốc. Tuy nhiên, số liệu mô phỏng chứng minh khi vượt quá ngưỡng $1.500\text{ kN/m}$, hệ số $F_s$ hầu như không đổi vì sức kháng cắt của nền đất đã đạt tới hạn dẻo cục bộ trước khi lực căng trong cốt kịp huy động thêm.

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

Có. Luận án trình bày chi tiết toàn bộ cơ sở toán học: hàm dạng phần tử tam giác đẳng tham số 3 nút và 6 nút, ma trận độ cứng cục bộ $[k_e]$, ma trận chuyển đổi tọa độ $[T]$, phương trình tiếp xúc phi tuyến, sơ đồ khối thuật toán lặp Newton-Raphson, và bảng thông số đầu vào chuẩn hóa ($c, \varphi, \psi, E_s, \nu, T_{max}, EA_g, t$). Bất kỳ nhà nghiên cứu nào cũng có thể tái lập mã nguồn trên nền tảng Fortran, C++ hoặc MATLAB.

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

Lộ trình 10 năm định hình sự chuyển đổi từ bài toán tĩnh 2D sang hệ thống tính toán động 3D đa vật lý (Multiphysics):

  1. Năm 1-3: Tích hợp mô hình từ biến - dão nhớt dẻo (creep-relaxation) dài hạn của polyme vào mã nguồn.
  2. Năm 4-6: Nâng cấp thuật toán HNH_RESS lên mô hình 3D biến dạng khối phục vụ phân tích mố cầu và nút giao ngầm.
  3. Năm 7-10: Xây dựng phần mềm mô phỏng tương tác động đất - thủy lực - hóa học (THMC Coupling) cho đê biển và công trình giao thông ven biển thích ứng với biến đổi khí hậu.

Kết luận

Luận án tiến sĩ của tác giả Huỳnh Ngọc Hào là một công trình nghiên cứu khoa học công phu, nghiêm túc và có giá trị học thuật xuất sắc. Những đóng góp cốt lõi của công trình được tổng kết qua 05 điểm then chốt:

  1. Hệ thống hóa toàn diện cơ sở lý thuyết và thực tiễn ứng dụng VĐKT trong xây dựng đường bộ tại Việt Nam và trên thế giới, phân định rõ nét 5 chức năng cơ bản: phân cách, gia cường, làm cốt tường chắn, lọc thoát nước và chống xói mòn.
  2. Xây dựng thành công mô hình toán học và thuật toán phần tử hữu hạn phi tuyến mô tả chính xác ứng xử kéo của VĐKT và bề mặt tiếp xúc đất - cốt theo mô hình thực nghiệm cao cấp Robert M. Koerner (2005).
  3. Độc lập lập trình và phát triển phần mềm HNH_RESS V1.00, đạt kết quả tương thích cao với phần mềm quốc tế Plaxis 2D (độ sai lệch chỉ $1,2% \div 3,5%$) nhưng tối ưu hóa khả năng can thiệp tham số cho kỹ sư Việt Nam.
  4. Chứng minh quy luật mặt trượt hình Ellipse đối với nền đắp gia cường trên đất yếu, xây dựng công thức giải tích tính toán chính xác lực căng cực đại $T_{max}$ huy động tại từng lớp vải địa.
  5. Thiết lập hệ thống biểu đồ nomograph tra cứu nhanh trong thiết kế nền đắp cao gia cường ($H = 6 \div 12\text{ m}$), lượng hóa tương quan giữa độ cứng $EA_g$, cường độ $T_{max}$ và mô đun đất $E_s$, mang lại hiệu quả kinh tế rõ rệt (giảm $15% \div 20%$ chi phí đầu tư).

Công trình đã mở ra các nhánh nghiên cứu mới trong lĩnh vực Địa kỹ thuật công trình giao thông, khẳng định năng lực tự chủ công nghệ tính toán số của các nhà khoa học Việt Nam và đóng góp trực tiếp vào sự nghiệp xây dựng mạng lưới hạ tầng giao thông hiện đại, bền vững.