Tổng quan luận án

Công trình nghiên cứu thuộc luận án tiến sĩ kỹ thuật chuyên ngành Xây dựng đường ô tô và đường thành phố (Mã số: 62.01), do nghiên cứu sinh Huỳnh Ngọc Hào thực hiện dưới sự hướng dẫn khoa học của GS. TS. Vũ Đình Phụng và TS. Vũ Đức Sỹ tại Trường Đại học Giao thông Vận tải (năm 2014).

Tính cấp thiết và khoảng trống nghiên cứu

Vải địa kỹ thuật (VĐKT - Geotextiles) là vật liệu polyme tổng hợp được phát triển và ứng dụng từ những năm 1970 tại các nước phương Tây và bắt đầu du nhập vào Việt Nam từ năm 1993 trong dự án nâng cấp Quốc lộ 5. Thực tiễn áp dụng trong xây dựng nền đường đắp cao và nền đắp trên đất yếu cho thấy việc sử dụng VĐKT đem lại hiệu quả kinh tế - kỹ thuật rõ rệt (giảm giá thành từ 15% đến 20%, thuận tiện trong thi công, tăng độ ổn định và tuổi thọ công trình).

Tuy nhiên, công tác tính toán thiết kế nền đắp có gia cường VĐKT tại Việt Nam còn tồn tại các khoảng trống kỹ thuật:

  1. Các phương pháp giải tích truyền thống dựa trên trạng thái cân bằng giới hạn chủ yếu giả thiết mặt trượt trụ tròn hoặc mặt trượt gãy khúc, chỉ xét đến cường độ chịu kéo cực đại ($T_{\max}$) của cốt mà bỏ qua độ cứng ($EA$) và mô đun đàn hồi ($E_g$) của vật liệu.
  2. Các phần mềm số thương mại (như Plaxis, Phase2) mặc định quan hệ ứng suất - biến dạng của VĐKT và phần tử tiếp xúc đất - cốt là mô hình đàn hồi dẻo lý tưởng tuyến tính theo tiêu chuẩn Mohr-Coulomb, chưa mô tả đúng bản chất ứng xử phi tuyến qua nhiều giai đoạn của vật liệu polyme. Đồng thời, cấu trúc đóng gói của các phần mềm ngoại nhập không cho phép can thiệp để tùy biến theo điều kiện cơ lý đất độn và vật liệu tại Việt Nam.

Do đó, việc nghiên cứu hoàn thiện phương pháp tính toán bằng phương pháp phần tử hữu hạn (PTHH) mô tả sát ứng xử cơ lý thực tế của hệ kết cấu "đất - cốt" và xây dựng phần mềm chuyên dụng là yêu cầu cấp thiết.

Mục tiêu nghiên cứu

Mục tiêu cốt lõi của đề tài là xây dựng mô hình tính toán nền đường đất đắp gia cường VĐKT bằng phương pháp phần tử hữu hạn, phản ánh chính xác quan hệ phi tuyến thực tế của vật liệu; dự báo chính xác khả năng mất ổn định và trạng thái ứng suất - biến dạng; từ đó thiết lập chương trình tính toán chuyên dụng và hệ thống biểu đồ tiện ích phục vụ thiết kế công trình đường ô tô.

Nhiệm vụ nghiên cứu

Luận án xác lập 3 nội dung nghiên cứu chính:

  1. Tổng kết thực tiễn ứng dụng và các mô hình tính toán nền đắp gia cường VĐKT trên thế giới và tại Việt Nam.
  2. Thiết lập mô hình tính toán và thuật toán bài toán ổn định nền đắp gia cường VĐKT theo phương pháp phần tử hữu hạn có xét đến quan hệ phi tuyến của cốt mềm và phần tử tiếp xúc.
  3. Lập trình xây dựng phần mềm HNH_RESS V1.00, kiểm chứng với các phần mềm quốc tế (Plaxis), tiến hành thực nghiệm số trên các chiều cao đắp khác nhau (6m, 8m, 10m, 12m) trên nền đất tốt và đất yếu, xác lập công thức giải tích tính lực căng cực đại ($T_{\max}$) theo dạng mặt trượt Ellipse và xây dựng các biểu đồ tra cứu thiết kế.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Kết cấu nền đất đắp có sử dụng cốt mềm vải địa kỹ thuật trong các công trình xây dựng nền đường ô tô.
  • Phạm vi nghiên cứu: Mô hình hóa bài toán biến dạng phẳng của nền đường đắp cao và nền đắp trên đất yếu có gia cường VĐKT; xây dựng thuật toán phần tử hữu hạn 2D; khảo sát các trường hợp chiều cao nền đắp từ 6m đến 12m với các thông số độ dốc mái taluy, cường độ kéo và mô đun đàn hồi của đất và cốt.

Tổng quan tài liệu và vị trí của luận án

Lịch sử phát triển và thực tiễn ứng dụng vải địa kỹ thuật

VĐKT xuất hiện lần đầu tiên dưới tên thương mại Bidim vào những năm 1960 tại Pháp. Giai đoạn 1965–1973, Hà Lan đã sử dụng 5–6 triệu $\text{m}^2$ VĐKT trong dự án đê biển chống xâm nhập mặn. Từ thập niên 1980, vật liệu này được áp dụng rộng rãi tại Pháp, Hà Lan, Mỹ, Nhật Bản, Trung Quốc, Ấn Độ, Hàn Quốc và các quốc gia Đông Nam Á (Thái Lan, Philippines, Indonesia, Malaysia, Brunei).

Tại Việt Nam, VĐKT được đưa vào sử dụng từ năm 1993 tại dự án cải tạo Quốc lộ 5 do công ty tư vấn KEI (Nhật Bản) thiết kế với khối lượng hơn 500.000 $\text{m}^2$ nhằm xử lý nền đất yếu. Tiếp đó, VĐKT được ứng dụng trong hàng loạt dự án giao thông trọng điểm:

  • Tuyến nâng cấp, cải tạo: Quốc lộ 1, QL10, QL18, QL3, QL51, QL32, QL38, QL39; đoạn tuyến Trới – Vũ Oai (Quảng Ninh, 2000–2003) và tuyến tránh TP. Vinh trên QL1 (2002).
  • Tuyến cao tốc: Cao tốc TP. Hồ Chí Minh – Trung Lương, Đại lộ Thăng Long, Hà Nội – Hải Phòng, Cầu Giẽ – Ninh Bình, Nội Bài – Lào Cai, Hà Nội – Thái Nguyên, Long Thành – Dầu Giây, Bến Lức – Long Thành.
  • Công trình môi trường và bảo vệ mái dốc: Xử lý đáy bãi rác Tam Tân (Củ Chi, TP.HCM, 2011), bãi rác Bố Trạch (Quảng Bình), bãi rác An Nhơn (Bình Định), bãi rác TP. Hưng Yên; chống xói bờ biển Bãi Cháy (Quảng Ninh); kè mái hồ trung tâm TP. Lào Cai (2000), hồ điều hòa Trung tâm Hội nghị Quốc gia Hà Nội; xử lý xói mái dốc mố cầu Phù Đổng, cầu Như Nguyệt, cầu Xương Giang (QL1) do tư vấn PCI thiết kế (1998); tường chắn đất cao 7,2m tại đường dẫn cầu Tân Thịnh – Thị trấn Vôi (Bắc Giang, 1998).

Các hướng nghiên cứu và phương pháp tính toán hiện hữu

Tác giả đã hệ thống hóa các phương pháp tính toán nền đắp gia cường trong và ngoài nước:

Nhóm phương pháp Phương pháp cụ thể & Tác giả / Tiêu chuẩn Đặc điểm và hạn chế theo nhận định của luận án
Giải tích nền đất yếu Phân tích mặt trượt tròn (Bishop, Janbu cải tiến) xác định lực kéo yêu cầu $T_{roj}$ và chiều dài neo bám $L_j$. Giả thiết mặt trượt tròn; xét trạng thái cân bằng giới hạn; chỉ xét lực căng $T_{roj}$, chưa xét độ cứng $EA$ và mô đun biến dạng của cốt; phức tạp khi nền có nhiều lớp đất yếu.
Giải tích nền đất tốt - Khối nêm hai phần (mặt trượt gãy khúc).
- Phương pháp phân mảnh trượt tròn (K. Terzaghi, A.W. Bishop, N. Janbu, A. Nichiprovich, Lý thuyết độ ẩm; các tiêu chuẩn 22TCN 207-1992, 22TCN 219-1994, OCDI Nhật Bản, BS-6349 Part 1-1984).
- Phương pháp ứng suất kết hợp.
- Mặt trượt xoắn ốc logarit (Leschinsky & Boedecker 1989; Bridle & Barr 1990).
- Phương pháp trọng lực dính kết Rankine (Segrestin, P.).		 Phải giả định trước hình dạng mặt trượt (tròn, gãy khúc, xoắn ốc); tìm kiếm mặt trượt nguy hiểm nhất qua thử dần nên khó đạt độ chính xác cao; bỏ qua tương tác biến dạng giữa đất và cốt; không xét đến mô đun đàn hồi $E$ của đất và độ cứng $EA$ của VĐKT.
Phương pháp số - Nguyên lý cực trị Gauss và sai phân hữu hạn (Hoàng Đình Trạm, Hà Huy Cương).
- Phần mềm Geo-Slope (Slope/W).
- Phần mềm Plaxis (Hà Lan) và Phase2 (Canada).		 Phương pháp sai phân khó xử lý bài toán biên phức tạp. Các phần mềm thương mại Plaxis, Phase2 giả thiết quan hệ ứng suất - biến dạng của VĐKT và tiếp xúc đất - cốt là đàn hồi dẻo lý tưởng tuyến tính theo Mohr-Coulomb, không thể hiện được tính phi tuyến thực tế; người dùng không can thiệp được vào thuật toán.

Vị trí và giải pháp của luận án

Khắc phục các nhược điểm trên, luận án lựa chọn phương pháp phần tử hữu hạn với mô hình tiếp xúc và ứng xử kéo phi tuyến theo mô hình Robert M. Koerner (2005), xác định mặt trượt thực tế dựa trên trường chuyển vị nút (không cần giả định trước cung trượt tròn) và phát triển chương trình tính HNH_RESS V1.00 độc lập.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết biến dạng phẳng và mô hình vật liệu

  1. Lý thuyết đàn hồi bài toán biến dạng phẳng: Xét khối đất nền và thân đường đắp có kích thước theo phương dọc trục (trục $Z$) rất lớn so với mặt cắt ngang, biến dạng dọc trục $\varepsilon_z = 0$ và các biến dạng trượt $\gamma_{yz} = \gamma_{zx} = 0$. Mối quan hệ ứng suất – biến dạng trong mặt phẳng $XY$ tuân theo định luật Hooke tổng quát thông qua mô đun đàn hồi $E$ và hệ số Poisson $\nu$: $$\varepsilon_x = \frac{1-\nu^2}{E}\left(\sigma_x - \frac{\nu}{1-\nu}\sigma_y\right), \quad \varepsilon_y = \frac{1-\nu^2}{E}\left(\sigma_y - \frac{\nu}{1-\nu}\sigma_x\right), \quad \gamma_{xy} = \frac{2(1+\nu)}{E}\tau_{xy}$$

  2. Mô hình đất nền Mohr-Coulomb: Đất đắp và đất nền tự nhiên được mô hình hóa là môi trường đàn hồi dẻo lý tưởng, đặc trưng bởi các chỉ tiêu: trọng lượng thể tích $\gamma$, lực dính đơn vị $c$, góc ma sát trong $\varphi$, góc giãn nở $\psi$, mô đun đàn hồi $E_s$ và hệ số Poisson $\nu$.

  3. Mô hình ứng xử kéo phi tuyến của VĐKT (Robert M. Koerner, 2005): VĐKT được phân loại gồm hai nhóm chính:

  • VĐKT dệt (Woven): Robusta Nicolon (Hà Lan), Krafter (Nhật Bản), Amoco (Anh), cường độ chịu kéo cao theo phương cuộn, tính thấm nhỏ.
  • VĐKT không dệt (Nonwoven): Fibertex (Đan Mạch), Polyfelt TS (Malaysia/Úc), Terrafix (Canada), TS, ART, VNT, PH, HD (Việt Nam), cường độ đẳng hướng, tính thấm và khả năng lọc cao.

Ứng xử kéo của VĐKT được mô tả qua các thông số cơ lý (theo ASTM D5199 cho độ dày $d = 0{,}25 \div 3{,}5\text{ mm}$; ASTM D1388 cho độ cứng uốn):

  • Cường độ chịu kéo $T_{\max}\text{ (kN/m)}$ và ứng suất kéo lớn nhất $\sigma_{\max} = T_{\max}/d\text{ (kN/m}^2\text{)}$.
  • Mô đun tiếp tuyến ban đầu ($E_i$), mô đun tiếp tuyến trục tung dời ngang ($E_{OT}$ áp dụng cho vải không dệt khâu lỗ kim), và mô đun cát tuyến ($E_{S10}, E_{S35}$).
  • Độ cứng kéo của phần tử thanh: $EA = E_g \cdot A\text{ (kN)}$, với $A$ là diện tích tiết diện danh định theo chiều dày vải.
  1. Mô hình phi tuyến tiếp xúc Đất – VĐKT: Ứng xử trượt tiếp xúc trong thí nghiệm kéo trượt tuân theo 3 giai đoạn:
  • Giai đoạn phi tuyến (0–1): Ứng suất tiếp tăng nhanh, biến dạng trượt tăng chậm.
  • Giai đoạn tái bền (Hardening, 1–2): Ứng suất tiếp và biến dạng đồng thời tăng đến giá trị đỉnh.
  • Giai đoạn hóa mềm (Softening, 2–3): Ứng suất tiếp suy giảm khi biến dạng tiếp tục tăng.
  1. Phương pháp giảm tham số cường độ chống cắt ($c-\varphi$ reduction): Hệ số an toàn ổn định ($F_s$) được xác định thông qua việc giảm dần các chỉ tiêu lực dính $c$ và góc ma sát trong $\varphi$ của đất cho đến khi hệ thống đạt trạng thái giới hạn phá hoại: $$F_s = \frac{\tan\varphi_i}{\tan\varphi_r} = \frac{c_i}{c_r} = \frac{s_{u,i}}{s_{u,r}}$$ $$\varphi_r = \arctan\left(\frac{\tan\varphi_i}{F_s}\right); \quad c_r = \frac{c_i}{F_s}; \quad s_{u,r} = \frac{s_{u,i}}{F_s}$$

Kỹ thuật phần tử hữu hạn và thuật toán giải phi tuyến

  • Phần tử môi trường liên tục (Đất): Sử dụng phần tử tam giác 3 nút và tam giác 6 nút đẳng tham số với phép tích phân số trên miền tam giác.
  • Phần tử cốt mềm VĐKT: Mô hình hóa bằng phần tử thanh 1D chịu lực dọc trục phi tuyến (chỉ chịu kéo, độ cứng chịu nén bằng 0).
  • Phần tử tiếp xúc (Interface element): Mô phỏng tương tác dính bám - trượt giữa bề mặt VĐKT và các phần tử đất liền kề.
  • Thuật toán phi tuyến: Sử dụng thuật toán lặp Newton-Raphson (N-R) và Newton-Raphson cải tiến để giải hệ phương trình cân bằng phi tuyến từng bước tải trọng.

Nội dung chính theo từng chương

Chương 1: Tổng quan về tình hình sử dụng và phương pháp tính toán nền đắp đường vải địa kỹ thuật

Chương 1 tổng thuật lịch sử phát triển, phân loại công nghệ chế tạo (dệt và không dệt), tiêu chuẩn đánh giá (ASTM, ISO) và 7 chức năng cơ bản của VĐKT trong công trình giao thông: phân cách, gia cường nền đất yếu, cốt tường chắn đất, lọc - tiêu thoát nước lưng tường, hạ mực nước ngầm, chống xói mòn bảo vệ mái dốc bờ sông/biển và làm ống địa kỹ thuật (Geotextile tubes). Chương này phân tích chi tiết các dự án áp dụng tại Việt Nam (QL5, QL1, các tuyến cao tốc, bãi rác xử lý chất thải).

Về mặt lý thuyết, chương 1 đối chiếu chi tiết các công thức giải tích cân bằng giới hạn:

  • Công thức kiểm tra ổn định trượt cục bộ: $\frac{H}{L_s} \le \frac{\tan\varphi'{cv}}{f{ms}}$
  • Phương pháp phân tích mô men giữ của cốt: $M_{RR} = \sum T_{roj} \cdot Y_j$
  • Chiều dài neo bám tối thiểu chống tuột cốt trên nền đất yếu: $$T_{roj} \le \frac{1}{f_n f_p}\left(2h\gamma \tan\varphi'{cv}\frac{\alpha'}{f{ms}} + 2c_u\alpha'_{bc}L_j\right)$$
  • Phương pháp khối nêm hai phần trên nền đất tốt với lực gây trượt tổng hợp: $R_h = \frac{1}{2}K f_{fs}\gamma H^2$ và khoảng cách bố trí cốt theo phương đứng: $S_{vj} \le \frac{T_j}{f_{fs}\gamma h_j K + f_q w_s K}$.

Đồng thời, chương 1 chỉ ra các bất cập của phần mềm Plaxis và Phase2 khi tuyến tính hóa mô hình VĐKT theo Mohr-Coulomb, làm tiền đề xác lập mục tiêu nghiên cứu của luận án.

Chương 2: Mô hình tính bài toán nền đất đắp đường bằng cốt mềm vải địa kỹ thuật

Chương 2 thiết lập hệ thống giả thiết cơ bản cho bài toán nền đắp nhiều lớp đất trên nền tự nhiên đồng nhất hoặc phân tầng. Cốt VĐKT được đặt nằm ngang tại các cao độ thiết kế.

Chương này trình bày chi tiết các ví dụ định lượng xác định thông số cơ lý từ đường cong quan hệ ứng suất – biến dạng của Robert M. Koerner:

  • Ví dụ mẫu VĐKT không dệt nhiệt hóa (Loại D): Chiều dày $d = 0{,}33\text{ mm}$, trọng lượng $135\text{ g/m}^2$, lực căng phá hoại $T_{\max} = 23\text{ kN/m}$, ứng suất kéo cực đại $\sigma_{\max} = 69.700\text{ kPa}$, độ giãn dài khi đứt $\varepsilon = 69%$, độ dẻo $U_g = 24.000\text{ kPa}$, mô đun đàn hồi ban đầu $E_g = 139.390\text{ kPa}$.
  • Ví dụ mẫu VĐKT không dệt khâu lỗ kim (Loại E): Chiều dày $d = 0{,}63\text{ mm}$, trọng lượng $200\text{ g/m}^2$, xác định mô đun cát tuyến tại các mức biến dạng $10%$ và $35%$ ($E_{S10}, E_{S35}$) và mô đun tiếp tuyến trục tung dời ngang ($E_{OT}$).

Chương 2 hoàn thiện hệ phương trình vi phân đàn hồi – dẻo trong điều kiện biến dạng phẳng và tích hợp phương pháp suy giảm $c-\varphi$ vào quy trình phân tích PTHH.

Chương 3: Xây dựng thuật toán và chương trình tính ổn định nền đất đắp đường vải địa kỹ thuật bằng phương pháp PTHH

Chương 3 chi tiết hóa ma trận độ cứng cho phần tử tam giác 3 nút và tam giác 6 nút đẳng tham số, thiết lập ma trận phần tử thanh 1D cho cốt mềm VĐKT và phần tử tiếp xúc phi tuyến mô tả 3 giai đoạn (phi tuyến, tái bền, hóa mềm).

Nội dung trọng tâm là việc xây dựng hoàn chỉnh phần mềm chuyên dụng HNH_RESS V1.00 (Reinforced Embankment Stability Software) với các module chức năng:

  • Module nhập dữ liệu hình học, phân lớp địa chất, sơ đồ xếp tải xe tiêu chuẩn.
  • Module khai báo phi tuyến đường cong ứng suất – biến dạng và độ cứng cát tuyến $EA_g$ của VĐKT.
  • Module giải thuật phi tuyến Newton-Raphson tự động phân chia bước lặp.
  • Module xử lý hậu kỳ: vẽ trường biến dạng, tự động nhận dạng cung trượt nguy hiểm nhất, so sánh độ khớp sai số giữa mặt trượt dạng cung tròn cổ điển và dạng hình học Ellipse.

Chương 4: Thực nghiệm tính toán nền đắp đường VĐKT trong xây dựng đường ô tô

Chương 4 thực hiện hàng loạt bài toán mô phỏng thực nghiệm số trên nền đất tốt và nền đất yếu với chiều cao đắp $H = 6\text{m}, 8\text{m}, 10\text{m}, 12\text{m}$, tải trọng xe thiết kế và các góc mái dốc $1:1$, $1:1{,}25$, $1:1{,}5$:

Sơ đồ mặt cắt nền đường đắp gia cường nhiều lớp VĐKT và cung trượt xấp xỉ Ellipse:
+-------------------------------------------------------------+
|                     Mặt đường / Tải trọng xe                |
|               ====================================          |
|              /                                    \         |
|             /   Lớp VĐKT n (Tmax, EAg)             \        |
|            /----------------------------------------\       |
|           /     Lớp VĐKT 2 (khoảng cách Sv)          \      |
|          /--------------------------------------------\     |
|         /       Lớp VĐKT 1 (đáy nền đắp)               \    |
|  ======/------------------------------------------------\===|  <-- Mặt đất tự nhiên
|       (...... Cung trượt thực tế dạng Ellipse ..........)   |
|        (...............................................)    |
+-------------------------------------------------------------+

Các kết quả và quy luật rút ra từ thực nghiệm số bao gồm:

  1. Phân tích ổn định nền đắp trên đất tốt:

    • Khi tăng chiều cao đắp từ 6m đến 12m, việc bố trí VĐKT làm tăng đáng kể hệ số an toàn $F_s$.
    • Bố trí các lớp VĐKT với khoảng cách dày ở phần dưới chân ta-luy mang lại hiệu quả giữ ổn định cao hơn so với rải đều ở phần đỉnh. Khảo sát khoảng cách bố trí giữa các lớp ($S_v = 0{,}5\text{m}$ so với $S_v = 1{,}5\text{m}$) cho thấy cự ly $0{,}5\text{m}$ khống chế biến dạng trượt cục bộ tốt hơn rõ rệt.
    • Luận án đã lập các bảng tính chi tiết về lực căng phát sinh trong từng lớp VĐKT khi mái dốc đạt trạng thái giới hạn phá hoại tương ứng với các cấp cường độ ($T = 12, 14, 16, 18, 20, 22, 24, 26, 28\text{ kN/m}$).

  2. Dạng mặt trượt thực tế và công thức tính lực căng $T_{\max}$ trên nền đất yếu:

    • Phân tích trường chuyển vị nút từ HNH_RESS V1.00 chứng minh rằng đối với nền đắp có cốt trên nền đất yếu, mặt trượt thực tế phát triển sâu xuống các lớp đất yếu bên dưới và có dạng hình học cung Ellipse, sai khác rõ rệt so với giả định cung tròn Bishop cổ điển.
    • Luận án đã xây dựng công thức giải tích giải bài toán phân mảnh xác định lực căng lớn nhất ($T_{\max}$) cho từng lớp VĐKT tương ứng với cung trượt Ellipse.

  3. Ảnh hưởng tương hỗ của độ cứng $EA_g$, cường độ $T_{\max}$ và mô đun đàn hồi đất $E_s$:

    • Hệ số an toàn $F_s$ tăng tỷ lệ thuận với độ cứng $EA_g$ và cường độ $T_{\max}$ của VĐKT, đồng thời chịu ảnh hưởng mạnh bởi mô đun đàn hồi $E_s$ của đất đắp.
    • Luận án lập hệ thống đồ thị quan hệ $E_s - EA_g - T_{\max} - F_s$ phục vụ tra cứu nhanh.

  4. Kiểm chứng độ tin cậy giữa HNH_RESS V1.00 và Plaxis: Kết quả so sánh hệ số an toàn ổn định $F_s$ trên cùng điều kiện bài toán cho thấy sự tương thích rất cao giữa hai phần mềm:

Trường hợp tính toán (Chiều cao nền đắp $H$) Hệ số an toàn $F_s$ (Phần mềm Plaxis) Hệ số an toàn $F_s$ (Phần mềm HNH_RESS V1.00)
Nền đắp cao $H = 6\text{m}$ (không gia cường VĐKT) 1,18 1,19
Nền đắp cao $H = 6\text{m}$ (có gia cường VĐKT) 1,42 1,44
Nền đắp cao $H = 8\text{m}$ (có gia cường VĐKT) 1,31 1,33
Nền đắp cao $H = 10\text{m}$ (có gia cường VĐKT) 1,24 1,25
Nền đắp cao $H = 12\text{m}$ (có gia cường VĐKT) 1,16 1,18

Sự sai lệch giữa HNH_RESS và Plaxis nằm trong khoảng dưới $2%$, khẳng định độ tin cậy của thuật toán và chương trình phần mềm do tác giả tự xây dựng.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và khoa học

  1. Đã 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 cho bài toán ổn định nền đường đắp có cốt mềm VĐKT, trong đó tích hợp mô hình phi tuyến của Robert M. Koerner cho ứng xử kéo của VĐKT và phần tử tiếp xúc đất – cốt (thay thế cho giả thiết đàn hồi dẻo tuyến tính Mohr-Coulomb thông thường).
  2. Phát hiện và chứng minh dạng hình học mặt trượt nguy hiểm nhất trong nền đắp gia cường VĐKT trên nền đất yếu có dạng cung Ellipse thay vì dạng cung tròn cổ điển.
  3. Thiết lập công thức giải tích mới xác định lực căng lớn nhất ($T_{\max}$) trong từng lớp VĐKT dựa trên mặt trượt hình Ellipse, giải quyết hạn chế của phương pháp phân mảnh truyền thống.
  4. Làm rõ quy luật định lượng về sự tương tác đồng thời giữa ba yếu tố: độ cứng VĐKT ($EA_g$), cường độ chịu kéo ($T_{\max}$) và mô đun đàn hồi của đất ($E_s$) đối với hệ số an toàn ổn định $F_s$.

Đóng góp mới về mặt thực tiễn và ứng dụng

  1. Xây dựng hoàn chỉnh bộ phần mềm HNH_RESS V1.00 có giao diện trực quan, cho phép kỹ sư ngành đường tại Việt Nam chủ động tính toán, phân tích trạng thái ứng suất – biến dạng và độ ổn định của nền đắp có gia cường VĐKT mà không phụ thuộc vào các phần mềm thương mại đắt tiền của nước ngoài.
  2. Thiết lập hệ thống biểu đồ quan hệ tiện ích giữa chiều cao nền đắp ($H = 6\text{m} \div 12\text{m}$), độ dốc mái taluy ($1:1 \div 1:1{,}5$), số lớp vải, khoảng cách bố trí và cấp cường độ $T_{\max}$, hỗ trợ tra cứu nhanh trong công tác thiết kế sơ bộ và tối ưu hóa kết cấu nền đắp đường ô tô.

Kiến nghị và giải pháp đề xuất

  • Khi thiết kế nền đắp gia cường VĐKT, bắt buộc phải tính toán đồng thời điều kiện bền kéo đứt và điều kiện chống tuột cốt (chiều dài neo bám $L_e$).
  • Ưu tiên bố trí các lớp cốt tập trung với cự ly hẹp ($S_v \le 0{,}5\text{m}$) ở nửa dưới thân nền đắp để triệt tiêu biến dạng trượt cung tròn và trượt trồi.
  • Đưa chương trình HNH_RESS V1.00 vào ứng dụng rộng rãi tại các đơn vị tư vấn thiết kế công trình giao thông trong nước.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của đề tài

  1. Đề tài giới hạn phạm vi khảo sát ở 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 của các đoạn nền đắp chuyển tiếp mố cầu hoặc nền đường cong có bán kính nhỏ.
  2. Tải trọng xe cộ trong mô hình được quy đổi tương đương về tải trọng phân bố tĩnh, chưa mô phỏng trực tiếp tác động động học trùng phục và hiện tượng tích lũy biến dạng dẻo dưới tác dụng của dòng xe tải nặng.
  3. Chưa đi sâu xét đến hiện tượng từ biến (creep) kéo dài theo thời gian của vật liệu polyme trong môi trường đất ngập nước bão hòa nhiều năm.

Hướng nghiên cứu tiếp theo

  • Phát triển module tính toán 3D mở rộng cho phần mềm HNH_RESS để giải quyết các bài toán mố cầu, cống chui và nền đường đắp trên sườn dốc nghiêng phức tạp.
  • Nghiên cứu tích hợp bài toán thấm cố kết phi tuyến đa pha (nước - khí - hạt đất) đồng thời với quá trình gia tải theo từng giai đoạn thi công nền đắp trên đất yếu.
  • Bổ sung quy luật suy giảm cường độ và từ biến của VĐKT theo thời gian vào thuật toán phần tử hữu hạn.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị ứng dụng và tham khảo cao cho nhiều đối tượng trong ngành xây dựng công trình giao thông và địa kỹ thuật:

  • Nghiên cứu sinh, học viên cao học và giảng viên: Tham khảo phương pháp xây dựng ma trận độ cứng phần tử tiếp xúc phi tuyến, kỹ thuật lập trình giải bài toán phần tử hữu hạn địa kỹ thuật bằng phương pháp suy giảm $c-\varphi$, và thuật toán nhận dạng mặt trượt phi tròn trong Chương 2 và Chương 3.
  • Kỹ sư tư vấn thiết kế đường ô tô: Sử dụng trực tiếp hệ thống biểu đồ tra cứu thiết kế ở Chương 4 để lựa chọn nhanh chủng loại VĐKT (dệt hoặc không dệt), số lớp cốt, khoảng cách rải và chiều dài neo bám cho các bài toán nền đắp cao từ 6m đến 12m.
  • Các cơ quan quản lý và chủ đầu tư dự án giao thông: Tham khảo giải pháp kỹ thuật gia cường VĐKT nhằm thay thế một phần giải pháp đào thay đất truyền thống hoặc xử lý nền đắp cao, giúp rút ngắn thời gian thi công và tiết kiệm từ 15% đến 20% chi phí xây dựng.

Câu hỏi thường gặp

1. Vải địa kỹ thuật được phân loại và đánh giá những chỉ tiêu cơ lý cốt lõi nào trong luận án?

Theo luận án, VĐKT được chia thành hai loại chính dựa trên công nghệ chế tạo: VĐKT dệt (Woven) có cường độ kéo cao nhưng tính thấm nhỏ (như Robusta Nicolon, Krafter, Amoco) và VĐKT không dệt (Nonwoven) có tính thấm và lọc tốt, chịu lực đẳng hướng (như Polyfelt TS, Fibertex, ART, HD). Các chỉ tiêu cơ lý chủ yếu gồm: chiều dày danh định ($d$ đo theo ASTM D5199), cường độ chịu kéo cực đại ($T_{\max}\text{ kN/m}$), độ giãn dài khi đứt ($\varepsilon%$), mô đun đàn hồi ban đầu ($E_i$), mô đun cát tuyến ($E_{S10}, E_{S35}$), mô đun tiếp tuyến dời trục ($E_{OT}$) và độ cứng kéo của phần tử ($EA$).

2. Điểm khác biệt căn bản giữa mô hình tính trong luận án và các phần mềm Plaxis, Phase2 là gì?

Các phần mềm thương mại Plaxis và Phase2 mặc định xem quan hệ ứng suất – biến dạng của VĐKT và phần tử tiếp xúc đất – cốt là mô hình đàn hồi dẻo lý tưởng tuyến tính theo tiêu chuẩn Mohr-Coulomb (khi đạt giới hạn kéo thì bị phá hoại ngay). Ngược lại, mô hình của luận án áp dụng lý thuyết phi tuyến của Robert M. Koerner (2005), mô tả chi tiết 3 giai đoạn làm việc thực tế của tiếp xúc đất – cốt (phi tuyến ban đầu, tái bền - hardening, và hóa mềm - softening), giúp kết quả phản ánh sát hơn trạng thái làm việc thực tế của công trình.

3. Tại sao luận án khẳng định mặt trượt của nền đắp có cốt trên đất yếu có dạng hình học Ellipse?

Qua phân tích trường chuyển vị và biến dạng trượt của các nút phần tử hữu hạn trên phần mềm HNH_RESS V1.00, kết quả cho thấy khi có cốt VĐKT gia cường ở đáy nền, cốt đã ngăn cản sự dịch chuyển ngang và triệt tiêu cung trượt tròn nông cục bộ. Vùng phá hoại bị đẩy sâu xuống tầng đất yếu bên dưới và kéo dài theo phương ngang, tạo thành biên trượt giới hạn có dạng hình học đường cong Ellipse. Luận án đã thiết lập giải thuật xấp xỉ sai số để chứng minh mặt trượt Ellipse khớp với thực tế biến dạng hơn mặt trượt tròn cổ điển.

4. Phần mềm HNH_RESS V1.00 do tác giả xây dựng có những chức năng nổi bật gì?

HNH_RESS V1.00 là phần mềm chuyên dụng phân tích ổn định nền đắp gia cường bằng phương pháp phần tử hữu hạn 2D. Phần mềm cho phép mô phỏng nền nhiều lớp đất, khai báo đặc tính phi tuyến của VĐKT và phần tử tiếp xúc; giải hệ phi tuyến bằng thuật toán Newton-Raphson cải tiến; tự động tính toán hệ số an toàn ổn định $F_s$ theo phương pháp giảm $c-\varphi$; tự động tìm kiếm và hiển thị mặt trượt nguy hiểm nhất (dạng tròn hoặc Ellipse) và tính toán lực căng lớn nhất trong từng lớp cốt. Kết quả của phần mềm đã được kiểm chứng tương thích với Plaxis (sai số $< 2%$).

5. Yếu tố độ cứng ($EA_g$) của VĐKT ảnh hưởng như thế nào đến độ ổn định nền đắp?

Phương pháp giải tích cổ điển chỉ xét đến cường độ kéo phá hoại $T_{\max}$ mà bỏ qua độ cứng $EA_g$. Nghiên cứu thực nghiệm số của luận án chứng minh rằng với cùng một giá trị lực căng $T_{\max}$, nếu VĐKT có độ cứng $EA_g$ lớn hơn thì khả năng khống chế biến dạng trượt ban đầu tốt hơn, làm giảm ứng suất cắt phát sinh trong nền đất và giúp nâng cao hệ số an toàn ổn định tổng thể $F_s$ của mái dốc nền đắp.

Kết luận

Luận án tiến sĩ của tác giả Huỳnh Ngọc Hào đã giải quyết thành công bài toán tính toán ổn định nền đường ô tô đắp cao và nền đắp trên đất yếu có gia cường vải địa kỹ thuật bằng phương pháp phần tử hữu hạn. Công trình đã đưa mô hình phi tuyến thực tế của vật liệu polyme và mặt trượt dạng cung Ellipse vào quy trình tính toán cơ học đất, khắc phục các giả định đơn giản hóa của phương pháp giải tích truyền thống. Việc phát triển thành công phần mềm chuyên dụng HNH_RESS V1.00 và thiết lập hệ thống biểu đồ tra cứu thiết kế mang lại đóng góp thiết thực cho công tác tư vấn thiết kế và thi công nền đường giao thông tại Việt Nam.