Tổng quan luận án

Công trình nghiên cứu "Nghiên cứu trạng thái ứng suất giới hạn trong nền đất tự nhiên dưới tác dụng của tải trọng nền đường đắp và bệ phản áp" được thực hiện bởi nghiên cứu sinh Nguyễn Minh Khoa, dưới sự hướng dẫn khoa học của TS. Hoàng Đình Đạm (cùng sự chỉ dẫn của GS.TSKH Hà Huy Cương), bảo vệ tại Học viện Kỹ thuật Quân sự năm 2013, thuộc chuyên ngành Kỹ thuật xây dựng công trình giao thông (mã số 62 58 02 05).

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

Hạ tầng giao thông đường bộ tại Việt Nam phân bố rộng khắp trên các vùng địa hình có cấu tầng địa chất phức tạp, đặc biệt là các vùng đồng bằng phù sa sông Hồng, đồng bằng sông Cửu Long, dải ven biển miền Trung cũng như các thung lũng, vùng trũng chân núi ở khu vực trung du và miền núi. Đặc trưng nổi bật của các khu vực này là sự hiện diện phổ biến của các tầng đất yếu có cường độ chống cắt thấp, tính nén lún lớn, hệ số rỗng cao và độ bão hòa nước lớn. Hiện tượng mất ổn định nền đường đắp trên đất yếu (dưới các dạng phình trồi hoặc trượt trồi) thường xuyên xảy ra trong giai đoạn thi công và khai thác, gây hư hỏng công trình và gia tăng chi phí xử lý.

Các lý thuyết tính toán cơ học đất truyền thống xác định trạng thái ứng suất và tải trọng giới hạn thường dựa trên các giả thiết:

  1. Đất là môi trường đàn hồi tuyến tính, đàn - dẻo hoặc cứng - dẻo lý tưởng;
  2. Áp dụng lời giải cân bằng giới hạn kinh điển của L. Prandtl vốn chỉ giải cho đất không có trọng lượng ($\gamma = 0$) và dưới tác dụng của móng cứng trơn nhẵn;
  3. Sử dụng các phương pháp gần đúng dựa trên mặt trượt giả định (mặt trượt phẳng, mặt trượt trụ tròn Bishop, Fellenius) hoặc phương pháp cộng tác dụng của Sokolovski, Terzaghi để xét đến trọng lượng bản thân đất.

Tuy nhiên, tải trọng nền đường đắp tác dụng lên nền đất tự nhiên có bản chất là tải trọng móng mềm (áp lực tiếp xúc phụ thuộc trực tiếp vào chiều cao và hình dạng hình học của thân đắp). Việc áp dụng mô hình móng cứng hoặc bỏ qua tương tác ứng suất thực tế làm sai lệch trạng thái làm việc của nền. Mặt khác, phương pháp coi nền đất là môi trường hạt rời ổn định theo điều kiện ứng suất tiếp lớn nhất đạt giá trị nhỏ nhất có khả năng phản ánh tốt hơn trạng thái ứng suất thực và xét được trọng lượng bản thân, nhưng các nghiên cứu trước đó còn hạn chế, chủ yếu mới khảo sát dưới móng cứng.

Giải pháp bệ phản áp là giải pháp tăng cường ổn định và tăng sức chịu tải của nền đất yếu được ứng dụng rộng rãi nhờ ưu điểm thi công đơn giản, tận dụng vật liệu tại chỗ. Mặc dù vậy, việc tính toán kích thước hình học (chiều rộng $L$, chiều cao $h$) và cường độ tải trọng ($q$) của bệ phản áp trong thực tế vẫn phụ thuộc chủ yếu vào kinh nghiệm và các công thức gần đúng, thiếu một phương pháp giải tích - số chuẩn xác xác định trạng thái ứng suất giới hạn dưới tải trọng móng mềm để tối ưu hóa kích thước thiết kế.

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

Luận án xác định hai mục tiêu nghiên cứu cụ thể:

  1. Xác định trạng thái ứng suất giới hạn trong nền đất tự nhiên dưới tác dụng của tải trọng nền đường đắp và bệ phản áp, dựa trên giả thiết nền đất mang tính chất của môi trường hạt rời, ổn định theo điều kiện ứng suất tiếp lớn nhất đạt giá trị nhỏ nhất, chịu tải trọng móng mềm và xét đến ảnh hưởng của trọng lượng bản thân nền đất.
  2. Ứng dụng bài toán trạng thái ứng suất giới hạn để nghiên cứu quy luật bệ phản áp làm tăng tải trọng giới hạn (sức chịu tải) của nền đất yếu, từ đó xác định kích thước và tải trọng hợp lý của bệ phản áp và xây dựng toán đồ phục vụ thiết kế thực hành.

Câu hỏi nghiên cứu

Luận án tập trung giải quyết các câu hỏi khoa học:

  1. Trạng thái ứng suất hữu hiệu và quy luật hình thành, phát triển vùng biến dạng dẻo trong nền đất tự nhiên diễn ra như thế nào khi chịu đồng thời trọng lượng bản thân và tải trọng móng mềm phân bố dạng hình thang của nền đường đắp?
  2. Khi bổ sung bệ phản áp hai bên thân đường, trạng thái ứng suất giới hạn và tải trọng giới hạn $p_{gh}$ của nền đất yếu biến đổi ra sao theo các thông số hình học ($L, h$) và cường độ tải trọng ($q$)?
  3. Tồn tại hay không ngưỡng kích thước và tải trọng hợp lý ($L_{hl}, q_{hl}$) mà tại đó hiệu quả tăng tải trọng giới hạn đạt mức tối ưu, và làm thế nào để số hóa kết quả này thành công cụ toán đồ thiết kế thuận tiện?

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

  • Đối tượng nghiên cứu: Trạng thái ứng suất, trạng thái ứng suất giới hạn, tải trọng giới hạn ($p_{gh}$) của nền đất tự nhiên và các thông số tải trọng, kích thước hình học của bệ phản áp ($L, h, q$).
  • Phạm vi nghiên cứu:
    • Mô hình hóa bài toán phẳng (trạng thái biến dạng phẳng) cho nền đất tự nhiên đồng nhất, đẳng hướng, có mặt thoáng nằm ngang.
    • Trạng thái ứng suất được phân tích là trạng thái ứng suất hữu hiệu ($\sigma'_x, \sigma'z, \sigma'{xz}$).
    • Tải trọng tác dụng gồm trọng lượng bản thân nền đất, tải trọng móng mềm của nền đường đắp và tải trọng bệ phản áp.

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

Luận án đã hệ thống hóa các hướng nghiên cứu trong và ngoài nước liên quan đến bài toán trạng thái ứng suất, sức chịu tải giới hạn của nền đất và giải pháp bệ phản áp:

Hướng nghiên cứu / Nhóm phương pháp Tác giả, công trình tiêu biểu được điểm luận Nội dung chính và đặc điểm
Lý thuyết biến dạng tuyến tính và đàn hồi Boussinesq, Flamant, N. Puzyrevsky, N. Maslov, Iaropolxki, V. Gorbunov – Poxadov, G. Frohlich, Đặng Hữu Xác định ứng suất trong bán không gian đàn hồi; xác định tải trọng an toàn $p_0$ khi vùng dẻo bắt đầu xuất hiện ($z_{\max}=0$ theo Puzyrevsky); khống chế độ sâu vùng dẻo ($z_{\max} = b\tan\varphi$ theo Maslov; $z_{\max}$ lớn nhất theo Iaropolxki); kiểm tra ổn định theo ứng suất cắt hoạt động $\tau_a$ (Đặng Hữu).
Lý thuyết cân bằng giới hạn kinh điển L. Prandtl (1920), Novotortsev (1938), Sokolovski Prandtl đưa ra lời giải giải tích tải trọng giới hạn dưới móng băng cứng trơn nhẵn cho đất không trọng lượng ($\gamma=0$); Novotortsev mở rộng cho tải trọng nghiêng; Sokolovski đề xuất phương pháp số giải hệ phương trình vi phân cân bằng có xét trọng lượng bản thân theo nguyên lý cộng tác dụng.
Lý thuyết đàn - dẻo và mô hình Cam-Clay Schofield & Wroth (1968), Roscoe & Burland (1968), Britto & Gunn (1987), Muir Wood (1991) Mô phỏng ứng xử của đất từ giai đoạn đàn hồi đến phá hoại dẻo; xét quy luật nén logarit và mặt bao trạng thái hóa mềm/thế dẻo.
Phương pháp mặt trượt giả định W. Gerxevanov, K. Florin, Bishop, K. Terzaghi, V. Berezansev, Vesic (1973), P. Ebdokimov (1964) Giả định mặt trượt phẳng, gãy khúc hoặc trụ tròn; Terzaghi thiết lập công thức sức chịu tải 3 thành phần ($N_\gamma, N_q, N_c$); Berezansev đưa ra mô hình nêm đất hình tam giác cân góc đáy $\pi/4$; Ebdokimov kết hợp lý thuyết Novotortsev và nguyên lý Caquot.
Phương pháp phân tích giới hạn Định lý giới hạn dưới (tĩnh học) và định lý giới hạn trên (động học) Xác định cận trên và cận dưới của tải trọng phá hoại mà không cần giải chi tiết phương trình vi phân đàn - dẻo phức tạp.
Phương pháp cực trị ứng suất tiếp nhỏ nhất Ngô Thị Thanh Hương Xem đất là môi trường hạt rời; bổ sung điều kiện phiếm hàm năng lượng biến dạng trượt đạt cực tiểu $\int_V (\tau_{\max}^2 / G) dV \to \min$ kết hợp phương trình cân bằng tĩnh để khép kín hệ phương trình vi phân.

Khoảng trống luận án lựa chọn giải quyết

Mặc dù các nghiên cứu trên đã đóng góp nền tảng cho cơ học đất, vẫn tồn tại các hạn chế kỹ thuật:

  1. Nghiệm giải tích của Prandtl và hầu hết công thức sức chịu tải móng nông áp dụng cho móng cứng, trong khi nền đường đắp truyền tải trọng xuống nền đất dưới dạng móng mềm có phân bố áp lực biến đổi theo mặt cắt hình học;
  2. Việc xét trọng lượng bản thân nền đất bằng phương pháp cộng tác dụng (như Terzaghi, Sokolovski) mang tính gần đúng;
  3. Phương pháp cực trị ứng suất tiếp của Ngô Thị Thanh Hương mới dừng ở móng cứng và chưa nghiên cứu trạng thái ứng suất giới hạn dưới tải trọng móng mềm;
  4. Việc tính toán thiết kế bệ phản áp cho nền đường đắp chưa có công cụ định lượng chuẩn xác dựa trên lý thuyết ứng suất giới hạn móng mềm có xét trọng lượng bản thân. Luận án chọn giải quyết trực tiếp khoảng trống này.

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

Cơ sở lý thuyết và khung phân tích

Luận án tiếp cận bài toán cơ học đất dựa trên mô hình môi trường hạt rời với hệ thống các phương trình và điều kiện biên:

  1. Điều kiện bền Mohr – Coulomb: Khối đất duy trì trạng thái cân bằng bền khi hàm giá trị bền $f(k)$ thỏa mãn: $$f(k) = \tau - \sigma \tan\varphi - c \le 0$$ Biểu diễn qua các thành phần ứng suất hữu hiệu trong bài toán phẳng: $$f(k) = \sqrt{\frac{(\sigma'_z - \sigma'x)^2}{4} + \tau{xz}^2} - \frac{\sigma'_z + \sigma'_x}{2}\sin\varphi - c\cos\varphi \le 0$$ Trạng thái cân bằng giới hạn (chảy dẻo) đạt được tại các điểm có $f(k) = 0$.

  2. Hệ phương trình vi phân cân bằng tĩnh: Trong hệ tọa độ phẳng $x0z$ ($z$ hướng xuống dưới), phương trình cân bằng tĩnh của phân tố đất có xét trọng lượng thể tích $\gamma$ là: $$\frac{\partial \sigma'z}{\partial z} + \frac{\partial \tau{xz}}{\partial x} = \gamma$$ $$\frac{\partial \sigma'x}{\partial x} + \frac{\partial \tau{xz}}{\partial z} = 0$$

  3. Nguyên lý cực trị ứng suất tiếp lớn nhất: Để khép kín hệ phương trình khi đất ở trạng thái hạt rời, luận án sử dụng điều kiện phiếm hàm đạt cực tiểu: $$\Phi = \iint_V \frac{\tau_{\max}^2}{G} dx dz \to \min$$ trong đó $\tau_{\max} = \sqrt{\frac{(\sigma'_z - \sigma'x)^2}{4} + \tau{xz}^2}$, $G$ là mô đun trượt của đất, $V$ là miền tính toán.

  4. Trạng thái ứng suất giới hạn: Xác định thông qua việc tăng dần tải trọng móng mềm $p$ cho đến khi vùng biến dạng dẻo ($f(k)=0$) phát triển liên tục và nối liền, cấu thành cơ chế phá hoại của nền.

Phương pháp nghiên cứu thực tế

  • Phương pháp số: Luận án sử dụng phương pháp sai phân hữu hạn (Finite Difference Method). Bán không gian nền đất được rời rạc hóa thành lưới sai phân chữ nhật với các bước lưới $\Delta x$ và $\Delta z$, số nút lưới theo phương ngang và đứng lần lượt là $n_a$ và $m_a$. Các đạo hàm riêng của các thành phần ứng suất được thay thế bằng các biểu thức sai phân hữu hạn tương ứng tại từng nút lưới $(i, j)$.
  • Thuật toán quy hoạch phi tuyến: Chuyển bài toán tìm nghiệm ứng suất về bài toán tối ưu phi tuyến có điều kiện: hàm mục tiêu là tổng năng lượng biến dạng trượt trên toàn bộ các nút lưới sai phân, ràng buộc là hệ phương trình vi phân cân bằng dạng sai phân, các điều kiện biên ứng suất trên mặt thoáng và điều kiện không vượt quá giới hạn bền Mohr - Coulomb.
  • Công cụ lập trình: Toàn bộ thuật toán sai phân hữu hạn và quy hoạch phi tuyến được lập trình tự động hóa bằng ngôn ngữ Matlab.

Nguồn dữ liệu và kiểm chứng

  • Dữ liệu địa kỹ thuật đất yếu: Sử dụng các dải thông số cơ lý đặc trưng của đất yếu trầm tích ven biển và đồng bằng Việt Nam theo tiêu chuẩn 22 TCN 262:2000 và TCVN 4054:2005 (lực dính không thoát nước $c_u < 15\text{ kPa}$, góc ma sát trong $\varphi = 0^\circ \div 10^\circ$, hệ số rỗng $e > 1$, trọng lượng thể tích $\gamma$).
  • Dữ liệu kiểm chứng mô hình: So sánh kết quả tính toán số với nghiệm giải tích chính xác của L. Prandtl cho trường hợp đất không trọng lượng ($\gamma = 0, \varphi = 0$) để đánh giá mức độ hội tụ và sai số của thuật toán.

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

Chương 1: Tổng quan về trạng thái ứng suất và tải trọng giới hạn của nền đất tự nhiên dưới tác dụng của tải trọng nền đường đắp

Chương 1 tổng thuật các kiến thức cơ bản về nền đường đắp, đặc tính cơ lý đất yếu tại Việt Nam, cơ chế phá hoại nền đắp và các lý thuyết cơ học đất hiện hành.

  • Đặc điểm nền đường đắp và đất yếu: Nền đường ô tô thông thường có bề rộng mặt đắp $B \ge 6,0\text{ m}$, độ dốc ta luy thiết kế chuẩn $1/1,5$, chiều cao đắp tối thiểu từ $1,2 \div 1,5\text{ m}$. Đất yếu tại Việt Nam tập trung nhiều ở đồng bằng sông Hồng, sông Cửu Long và ven biển miền Trung, phân loại thành hai nhóm chính: nhóm khoáng vật (bùn sét, á sét trầm tích) và nhóm hữu cơ (than bùn hàm lượng hữu cơ $20 \div 80%$).
  • Cơ chế mất ổn định: Hai dạng mất ổn định điển hình gồm:
    1. Dạng phình trồi: Vùng dẻo xuất hiện cục bộ dưới đáy và đẩy trồi đất yếu sang hai bên taluy nhưng chưa cắt đứt thân nền đắp;
    2. Dạng trượt trồi: Hình thành mặt trượt liên tục xé rách nền đường và trượt trồi khối đất lớn.
  • Đặc trưng tải trọng móng mềm: Khác với móng cứng có độ lún đều và ứng suất đáy móng tập trung ở mép móng, tải trọng nền đắp là tải trọng móng mềm phân bố dạng hình thang với cường độ lớn nhất ở giữa và vát dần ra mép taluy.
  • Tổng hợp giải pháp bệ phản áp:
Đặc tính giải pháp Đánh giá đối với bệ phản áp
Tăng mức độ ổn định trong thi công Đạt yêu cầu
Tăng mức độ ổn định sau khi đắp xong Đạt yêu cầu
Giảm chuyển vị ngang và lực đẩy ngang của đất yếu Đạt yêu cầu
Mức độ phức tạp công nghệ Đơn giản, tận dụng đất tại chỗ, không cần đầm lèn cao
Yêu cầu mặt bằng chiếm dụng Trung bình đến lớn (cần giải phóng mặt bằng hai bên)
Chi phí và tính khả thi Chi phí trung bình, thời gian thi công nhanh, ít rủi ro

Chương 2: Nghiên cứu trạng thái ứng suất trong nền đất tự nhiên dưới tác dụng của trọng lượng bản thân và tải trọng nền đường đắp

Chương 2 thiết lập bài toán trạng thái ứng suất phẳng trong nền đất tự nhiên chịu đồng thời trọng lượng bản thân và tải trọng móng mềm của nền đường đắp, giải bằng phương pháp sai phân hữu hạn và quy hoạch phi tuyến trên Matlab.

  • Xây dựng bài toán:
    • Hệ trục tọa độ $x0z$ với gốc tại mặt thoáng, trục $z$ hướng xuống.
    • Phân tích trạng thái ứng suất tự nhiên chỉ chịu trọng lượng bản thân: Ứng suất đứng $\sigma'{oz} = \gamma z$, ứng suất ngang $\sigma'{ox} = k\gamma z$ (với $k$ là hệ số áp lực ngang), ứng suất tiếp $\tau'_{oxz} = 0$.
    • Khi chịu tải trọng nền đường đắp cường độ $p$, tải trọng móng mềm hình thang được đưa vào điều kiện biên mặt thoáng $z = 0$.
  • Mô hình sai phân hữu hạn: Miền tính toán bán không gian được giới hạn với biên đủ xa để không ảnh hưởng đến vùng tập trung ứng suất. Lập hệ phương trình đại số phi tuyến từ phương trình cân bằng tĩnh và điều kiện tối tiểu hóa phiếm hàm ứng suất tiếp $\Phi$.
  • Quy luật phát triển vùng biến dạng dẻo:
    • Khảo sát sự hình thành các điểm chảy dẻo ($f(k) \ge 0$) theo cấp tải trọng tăng dần $p$.
    • Biến dạng dẻo xuất hiện đầu tiên ở độ sâu nhất định dưới vai đường/chân ta luy (nơi có gradient ứng suất tiếp lớn nhất), sau đó phát triển mở rộng thành các dải dẻo đối xứng qua trục tim đường.
    • Phân tích chi tiết sự phân bố lại của các trường ứng suất chính $\sigma'_z, \sigma'_x$ và các đường đẳng bền $f(k)$ trong nền đất.

Chương 3: Nghiên cứu trạng thái ứng suất giới hạn trong nền đất tự nhiên dưới tác dụng của tải trọng nền đường đắp và bệ phản áp

Chương 3 tập trung xác định trạng thái ứng suất giới hạn và tải trọng giới hạn ($p_{gh}$) của nền đất khi có sự tham gia của bệ phản áp, đồng thời đánh giá độ nhạy của mô hình số.

  • Thiết lập bài toán trạng thái ứng suất giới hạn: Sử dụng trạng thái ứng suất do trọng lượng bản thân làm hệ so sánh nền tảng. Khi đặt tải trọng nền đắp $p$ và tải trọng bệ phản áp $q$ (chiều rộng $L$, chiều cao $h$), bài toán giải hệ phương trình sai phân hữu hạn để tìm giá trị tải trọng lớn nhất $p_{gh}$ mà hệ vẫn duy trì được trạng thái cân bằng tĩnh học.
  • Khảo sát ảnh hưởng của mạng lưới sai phân:
    • Kích thước ô lưới ($\Delta x, \Delta z$): Khảo sát cho thấy khi giảm kích thước ô lưới, độ mịn tăng lên và giá trị tải trọng giới hạn $p_{gh}$ tiệm cận về nghiệm ổn định.
    • Kích thước biên miền tính toán: Miền lưới phải đủ lớn ($n_a, m_a$ đủ lớn) để biên nhân tạo không cản trở sự phát triển tự nhiên của vùng biến dạng dẻo.
  • Kiểm chứng với lời giải L. Prandtl: Khi cho $\gamma = 0$, $\varphi = 0$, tải trọng bên $q = 0$, nghiệm số từ chương trình Matlab cho giá trị không thứ nguyên: $$\frac{p_{gh}}{c} \approx 5{,}14$$ sai số so với nghiệm giải tích của Prandtl ($p_{gh} = (\pi + 2)c \approx 5{,}1416c$) là dưới $0{,}1%$, khẳng định độ tin cậy của thuật toán quy hoạch phi tuyến và mã nguồn Matlab.
  • Ảnh hưởng của các yếu tố địa kỹ thuật và hình học:
    • Ảnh hưởng của góc ma sát trong $\varphi$: Tải trọng giới hạn $p_{gh}/c$ tăng phi tuyến rất nhanh khi $\varphi$ tăng.
    • Ảnh hưởng của chiều rộng nền đắp $b$: Dưới tải trọng móng mềm có xét trọng lượng bản thân $\gamma$, $p_{gh}$ phụ thuộc rõ rệt vào bề rộng $b$, khác biệt với lời giải không trọng lượng của Prandtl.
    • Ảnh hưởng của bệ phản áp đến vùng dẻo: Tải trọng bệ phản áp $q$ tạo ra áp lực hông giữ chặt chân khối trượt, ngăn cản sự phát triển trồi lên mặt thoáng của vùng dẻo, làm giảm chiều sâu lớn nhất $h_d$ và bề rộng lớn nhất $b_d$ của vùng biến dạng dẻo dưới thân nền đường chính.

Chương 4: Nghiên cứu bệ phản áp để làm tăng tải trọng giới hạn của nền đất yếu dưới tải trọng nền đường đắp

Chương 4 là phần ứng dụng và thực nghiệm số chuyên sâu, đưa ra các kết quả định lượng về mối quan hệ giữa tải trọng giới hạn của nền đất và các thông số bệ phản áp.

  • Quan hệ giữa tải trọng giới hạn $p_{gh}/c$ với chiều rộng ($L$) và cường độ ($q$) bệ phản áp:
    • Khi tăng chiều rộng tải trọng bệ phản áp $L$, tải trọng giới hạn $p_{gh}/c$ tăng nhanh trong giai đoạn đầu, sau đó tốc độ tăng chậm dần và đạt trạng thái tiệm cận bão hòa.
    • Khi tăng cường độ tải trọng bệ phản áp $q$ (tương ứng chiều cao bệ $h$), $p_{gh}/c$ tăng tuyến tính trong một khoảng nhất định, nhưng nếu $q$ quá lớn, bản thân bệ phản áp sẽ gây phá hoại cục bộ nền đất dưới chân bệ.
Mối quan hệ định lượng giữa tải trọng giới hạn và thông số bệ phản áp:
1. Trường hợp nền đất dính thuần túy (φ = 0):
   - Tải trọng giới hạn p_gh/c tăng tỷ lệ thuận với cường độ bệ phản áp q/c.
   - Chiều rộng hiệu quả của bệ phản áp L có giới hạn hữu hạn (L_hl).
2. Trường hợp nền đất có ma sát trong (φ > 0):
   - Sức chịu tải được nhân thêm hệ số khuếch đại: f(φ) > 1,0.
   - Hiệu quả gia cường của bệ phản áp tăng mạnh khi góc ma sát trong φ tăng.
  • Xác định kích thước và tải trọng hợp lý của bệ phản áp:
    • Chiều rộng hợp lý ($L_{hl}$): Là giá trị chiều rộng mà vượt qua nó, tải trọng giới hạn $p_{gh}$ không tăng thêm quá $1 \div 2%$.
    • Cường độ hợp lý ($q_{hl}$) và chiều cao hợp lý ($h_{hl}$): Mức tải trọng bệ phản áp tối ưu đảm bảo sức chịu tải của nền đường chính đạt cực đại mà không gây mất ổn định dưới bệ phản áp ($h_{hl} \le H_{gh}$).
  • Trường hợp bệ phản áp rộng vô hạn ($L \to \infty$): Thiết lập quan hệ cận trên của sức chịu tải $p_{gh}/c$ phụ thuộc vào $q/c$ và góc ma sát trong $\varphi$.
  • Toán đồ thiết kế bệ phản áp: Luận án tổng hợp toàn bộ kết quả số hóa thành các đường cong toán đồ kỹ thuật:
    1. Toán đồ xác định tải trọng giới hạn $p_{gh}/c$ theo kích thước bệ phản áp $(L, q)$;
    2. Toán đồ xác định chiều rộng hợp lý $L_{hl}$ và cường độ hợp lý $q_{hl}$;
    3. Toán đồ hiệu chỉnh theo góc ma sát trong $f(\varphi)$.

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

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

  1. Đã xây dựng hoàn chỉnh mô hình toán và phương pháp xác định trạng thái ứng suất cũng như trạng thái ứng suất giới hạn trong nền đất tự nhiên đồng nhất có xét trọng lượng bản thân dưới tác dụng của tải trọng móng mềm (nền đường đắp và bệ phản áp), dựa trên giả thiết đất là môi trường hạt rời ổn định theo điều kiện ứng suất tiếp lớn nhất đạt giá trị nhỏ nhất.
  2. Ứng dụng thành công phương pháp sai phân hữu hạn kết hợp thuật toán quy hoạch phi tuyến trên ngôn ngữ Matlab để giải bài toán cơ học đất phức tạp, khắc phục được các nhược điểm của giả thiết móng cứng kinh điển và phương pháp cộng tác dụng gần đúng.
  3. Làm sáng tỏ cơ chế cơ học về sự phát triển của các dải biến dạng dẻo và quy luật phân bố lại ứng suất hữu hiệu trong nền đất yếu khi có bệ phản áp trợ lực.

Đóng góp mới về thực tiễn và giải pháp đề xuất

  1. Đưa ra phương pháp định lượng xác định tải trọng giới hạn $p_{gh}$ của nền đường đắp trên đất yếu có bệ phản áp, thay thế cho các phương pháp tra bảng kinh nghiệm hoặc tính toán gần đúng mặt trượt trụ tròn trước đây.
  2. Xác định được các tiêu chí và giá trị cụ thể của chiều rộng hợp lý ($L_{hl}$) và cường độ hợp lý ($q_{hl}$, $h_{hl}$) của bệ phản áp, giúp các kỹ sư thiết kế tối ưu hóa hình học bệ phản áp, tiết kiệm khối lượng đất đắp và giảm thiểu diện tích giải phóng mặt bằng.
  3. Xây dựng hệ thống toán đồ thiết kế bệ phản áp trực quan, cho phép tra cứu nhanh sức chịu tải và thông số bệ phản áp tương thích với điều kiện địa chất công trình cụ thể ($c, \varphi, \gamma$).

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

Hạn chế của luận án

  • Mô hình nghiên cứu dừng lại ở bài toán phẳng hai chiều (2D), chưa mô phỏng được các yếu tố không gian ba chiều (3D) tại các vị trí đầu mố cầu hoặc đoạn chuyển tiếp.
  • Giả thiết nền đất đồng nhất, đẳng hướng, có mặt thoáng nằm ngang; chưa xét trường hợp nền đất có nhiều lớp xen kẹp với các chỉ tiêu cơ lý biến đổi phức tạp theo chiều sâu hoặc mặt đất tự nhiên nghiêng dốc.
  • Phân tích trạng thái ứng suất hữu hiệu tĩnh, chưa tích hợp quá trình cố kết thấm biến đổi theo thời gian ($c_v$, áp lực nước lỗ rỗng thặng dư tiêu tán dần) và chưa xét tải trọng động trùng phục của phương tiện giao thông.

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

  1. Mở rộng bài toán cho mô hình nền đất yếu nhiều lớp không đồng nhất và có xét đến mái dốc tự nhiên của địa hình.
  2. Phát triển bài toán liên kết thủy - cơ (Hydro-Mechanical coupling) xét đến sự tiêu tán áp lực nước lỗ rỗng theo thời gian trong quá trình đắp nền theo từng giai đoạn kết hợp bệ phản áp.
  3. Nghiên cứu bài toán không gian 3 chiều và tích hợp tác dụng của tải trọng động từ đoàn xe chạy trên mặt đường.

Giá trị tham khảo

Luận án mang lại giá trị tham khảo khoa học và ứng dụng thực tiễn:

  • Đối với nghiên cứu sinh, học viên cao học và giảng viên:
    • Tài liệu tham khảo hữu ích cho các chuyên ngành Kỹ thuật xây dựng công trình giao thông, Địa kỹ thuật công trình, Cơ học đất và Nền móng.
    • Phần xây dựng thuật toán sai phân hữu hạn và kỹ thuật lập trình quy hoạch phi tuyến trên Matlab (trình bày tại Chương 2 và Chương 3) là tài liệu tham khảo tốt cho việc giải các bài toán trường ứng suất và cơ học môi trường liên tục.
  • Đối với kỹ sư tư vấn thiết kế và các cơ quan quản lý dự án:
    • Cung cấp cơ sở khoa học định lượng để tính toán thiết kế bệ phản áp xử lý nền đường đắp qua vùng đất yếu (Chương 4).
    • Hệ thống toán đồ thiết kế có thể ứng dụng trực tiếp trong giai đoạn lập dự án và thiết kế kỹ thuật thi công nền đường ô tô, đường cao tốc qua vùng đồng bằng sông Hồng, sông Cửu Long nhằm tối ưu hóa chi phí xây dựng.

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

1. Tại sao luận án khẳng định tải trọng nền đường đắp là tải trọng móng mềm thay vì móng cứng?

Móng cứng (như móng bê tông cốt thép của công trình dân dụng) có độ cứng chống uốn vô cùng lớn, độ lún tại đáy móng là phẳng đều và ứng suất tiếp xúc đáy móng tập trung cao ở các mép móng. Ngược lại, nền đường đắp là công trình bằng đất đắp có độ cứng hữu hạn, biến dạng mềm uốn theo biến dạng của nền đất tự nhiên bên dưới. Áp lực tiếp xúc truyền xuống nền đất tỷ lệ trực tiếp với chiều cao cột đất đắp tại từng vị trí trắc ngang (phân bố hình thang), do đó mang bản chất cơ học của móng mềm.

2. Điều kiện ổn định bổ sung theo nguyên lý cực trị ứng suất tiếp được thiết lập như thế nào?

Để khép kín hệ phương trình vi phân tĩnh định cho môi trường hạt rời mà không cần các giả thiết đàn hồi tuyến tính cổ điển, luận án áp dụng điều kiện phiếm hàm năng lượng biến dạng trượt đạt cực tiểu: $$\Phi = \iint_V \frac{\tau_{\max}^2}{G} dx dz \to \min$$ trong đó ứng suất tiếp lớn nhất $\tau_{\max} = \sqrt{\frac{(\sigma'_z - \sigma'x)^2}{4} + \tau{xz}^2}$. Điều kiện này kết hợp cùng hai phương trình vi phân cân bằng tĩnh tạo thành hệ 3 phương trình xác định duy nhất 3 thành phần ứng suất hữu hiệu ($\sigma'_x, \sigma'z, \tau{xz}$).

3. Nghiệm số của luận án được kiểm chứng với lý thuyết cổ điển nào để khẳng định độ tin cậy?

Thuật toán và chương trình lập trình Matlab được kiểm chứng trực tiếp bằng cách so sánh với lời giải giải tích chính xác của L. Prandtl (1920) cho bài toán móng băng trên nền đất không trọng lượng ($\gamma = 0, \varphi = 0, q = 0$). Kết quả chạy số của luận án cho hệ số sức chịu tải $p_{gh}/c \approx 5{,}14$, hoàn toàn trùng khớp với giá trị lý thuyết của Prandtl là $(\pi + 2) \approx 5{,}1416$.

4. Khái niệm "kích thước hợp lý" ($L_{hl}, q_{hl}$) của bệ phản áp được xác định dựa trên nguyên tắc nào?

Kích thước hợp lý của bệ phản áp được xác định tại điểm uốn tiệm cận của đường cong quan hệ giữa tải trọng giới hạn $p_{gh}$ và kích thước bệ:

  • Chiều rộng hợp lý ($L_{hl}$): Là ngưỡng chiều rộng mà nếu tiếp tục đắp rộng thêm ($L > L_{hl}$), tải trọng giới hạn $p_{gh}$ của nền đường chính tăng không đáng kể (dưới $1 \div 2%$).
  • Cường độ hợp lý ($q_{hl}$ hay chiều cao $h_{hl}$): Là giá trị tải trọng bệ phản áp tối đa giúp kìm hãm vùng dẻo dưới thân đường chính mà không gây ra hiện tượng phá hoại mất ổn định dẻo cục bộ ngay dưới chân bản thân bệ phản áp.

5. Góc ma sát trong $\varphi$ ảnh hưởng như thế nào đến hiệu quả của bệ phản áp?

Khi xét góc ma sát trong $\varphi > 0$, sức chịu tải giới hạn của nền đất yếu tăng lên đáng kể so với trường hợp đất thuần dính ($\varphi = 0$). Luận án đưa ra hệ số hiệu chỉnh $f(\varphi) > 1{,}0$. Khi góc ma sát trong $\varphi$ tăng, bệ phản áp phát huy hiệu quả giữ chân khối trượt và gia tăng áp lực hông chủ động mạnh hơn, làm tăng đáng kể tải trọng giới hạn $p_{gh}/c$ của hệ thống nền đắp.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Minh Khoa đã giải quyết bài toán cơ học nền móng về xác định trạng thái ứng suất giới hạn trong nền đất tự nhiên dưới tác dụng của tải trọng móng mềm nền đường đắp và bệ phản áp có xét trọng lượng bản thân đất. Bằng việc kết hợp lý thuyết môi trường hạt rời theo nguyên lý cực trị ứng suất tiếp với phương pháp sai phân hữu hạn và quy hoạch phi tuyến trên Matlab, công trình đã làm sáng tỏ quy luật phát triển vùng dẻo và cơ chế tăng cường sức chịu tải khi có bệ phản áp. Các kết quả nghiên cứu định lượng cùng hệ thống toán đồ thiết kế bệ phản áp hợp lý là đóng góp có giá trị khoa học và thực tiễn kỹ thuật trong lĩnh vực xây dựng công trình giao thông qua vùng đất yếu.