Tổng quan về luận án

Nghiên cứu của nghiên cứu sinh Soukha Yakoshi tại Trường Đại học Giao thông Vận tải (Hà Nội, 2023) tập trung giải quyết bài toán phức tạp bậc cao trong chuyên ngành Kỹ thuật Xây dựng Công trình Giao thông và Địa kỹ thuật: "Đánh giá ổn định mái dốc và ổn định trượt sâu mố cầu trên móng nông bằng Phương pháp Cân bằng Giới hạn Tổng quát (Generalized Limit Equilibrium Method - GLEM) kết hợp tối ưu hóa phi tuyến". Trong bối cảnh hạ tầng giao thông hiện đại đối mặt với các điều kiện địa chất bất lợi, tác động phức hợp của tải trọng xe nặng, động đất và biến động thủy văn, việc xác định chính xác cơ chế phá hoại trượt sâu của mố cầu đặt trên móng nông đóng vai trò sống còn đối với an toàn công trình.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác lập: Các phương pháp phân tích ổn định lát cắt cổ điển như Fellenius (1936), Bishop rút gọn (1955) hay Janbu đơn giản (1954) thường đưa ra các giả định giản lược hóa về lực tương tác giữa các thỏi hoặc bó hẹp trong giả định mặt trượt cung tròn. Điều này dẫn đến sự sai lệch nghiêm trọng khi đánh giá các kết cấu phức hợp chịu tải trọng tập trung lớn và không đối xứng như mố cầu. Luận án đã khắc phục triệt để khoảng trống này bằng cách mở rộng phương pháp GLEM, thiết lập hệ phương trình ma trận cân bằng tĩnh học toàn phần $[A]{K} = {B}$ cho từng thỏi phân tố đất và khối mố, đồng thời giải bài toán tìm mặt trượt nguy hiểm nhất có hình dạng tùy ý thông qua thuật toán tối ưu hóa phi tuyến.

Nghiên cứu tập trung giải quyết ba câu hỏi nghiên cứu chính:

  1. RQ1: Cơ chế truyền lực và tương tác giữa các thỏi phân tố ($H_i, V_i$) và mặt đáy ($N_i, T_i$) thay đổi như thế nào dưới tác động đồng thời của tĩnh tải mố ($Q_M$), áp lực đất đắp ($Q_{soil}$), hoạt tải xe ($QH, BR, IM$) và lực quán tính động đất ($v, h$)?
  2. RQ2: Mô hình toán học ma trận GLEM kết hợp tối ưu hóa có thể định dạng chính xác mặt trượt phi tròn tới hạn và tìm kiếm giá trị hệ số an toàn nhỏ nhất ($F_{s\min}$) cũng như hệ số an toàn trung bình ($F_{s\text{med}}$) trong các trạng thái giới hạn cường độ hay không?
  3. RQ3: Mức độ ảnh hưởng của chiều cao đất đắp ($H_đ = 2.05\text{m}, 4.09\text{m}, 8.19\text{m}$) và mực nước ngầm bão hòa ($W'_i, W''_i, U, U'$) đến khả năng kháng trượt sâu của bệ móng nông là bao nhiêu?

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa Lý thuyết Cân bằng Giới hạn Tổng quát (GLEM), Tiêu chuẩn bền cắt đất Mohr-Coulomb và Tiêu chuẩn Thiết kế Cầu theo Hệ số Tải trọng và Sức kháng AASHTO LRFD / TCVN 11823. Phạm vi nghiên cứu bao quát các thí dụ kiểm chứng kinh điển về mái đất dính đồng nhất (trượt chân dốc, trượt sâu) và 18 kịch bản phân tích trượt sâu mố cầu trên móng nông trong cả giai đoạn thi công lẫn khai thác.

       +-------------------------------------------------------------+
       |   KHUNG LÝ THUYẾT VÀ MÔ HÌNH PHÂN TÍCH GLEM TOÀN DIỆN       |
       +-------------------------------------------------------------+
                                      |
         +----------------------------+----------------------------+
         |                                                         |
         v                                                         v
+-----------------------------+                           +-----------------------------+
|    MÔ HÌNH CƠ HỌC ĐẤT       |                           |     TỔ HỢP TẢI TRỌNG        |
| - Tiêu chuẩn Mohr-Coulomb   |                           | - Tĩnh tải: QM, Qsoil, DC   |
| - Cân bằng lực GLEM         |                           | - Hoạt tải: QH, BR, IM, LS  |
| - Áp lực nước: U, U'        |                           | - Động đất: v, h (k_h, k_v) |
+-----------------------------+                           +-----------------------------+
         |                                                         |
         +----------------------------+----------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       |   HỆ PHƯƠNG TRÌNH MA TRẬN PHI TUYẾN TOÀN CỤC [A]{K} = {B}   |
       |     Cân bằng tĩnh học: Sigma Ni = 0, Sigma Ti = 0           |
       |     Điều kiện biên: Ti = f(Ni, c, phi, Fs), Vi = f(Hi, Fsi) |
       +-------------------------------------------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       |    THUẬT TOÁN TỐI ƯU HÓA PHI TUYẾN (EXCEL SOLVER / GRG)     |
       |   Xác định hình học mặt trượt tối ưu: (Xi, Yi, alpha_i)    |
       |   Xuất kết quả: Fs_min (FSI = vô cùng) & Fs_med (FSI = 1.7) |
       +-------------------------------------------------------------+

Literature Review và Positioning

Nghiên cứu cân bằng giới hạn trong ổn định mái dốc và công trình địa kỹ thuật đã trải qua nhiều thập kỷ tiến hóa. Dòng lý thuyết xuất phát từ Fellenius (1936) và Taylor (1937) giả định mặt trượt trụ tròn và bỏ qua hoàn toàn lực tương tác giữa các thỏi, dẫn đến kết quả thiên về quá an toàn hoặc sai lệch lớn. Bishop (1955) đã cải tiến bằng cách xét đến lực pháp tuyến giữa các thỏi nhưng bỏ qua lực cắt giữa các thỏi, chỉ thỏa mãn điều kiện cân bằng mô-men. Ngược lại, Janbu (1954, 1973) thiết lập phương pháp lát cắt cho mặt trượt bất kỳ nhưng chỉ thỏa mãn điều kiện cân bằng lực.

Bước đột phá trong cơ học địa kỹ thuật được đánh dấu bởi Morgenstern & Price (1965) và Spencer (1967) khi đề xuất mô hình cân bằng giới hạn tổng quát thỏa mãn đồng thời cân bằng lực và cân bằng mô-men. Morgenstern và Price đưa ra hàm giả định phân bố lực cắt giữa các thỏi $V = \lambda \cdot f(x) \cdot H$. Tuy nhiên, tranh biện khoa học kéo dài giữa hai trường phái:

  • Trường phái hàm giả định lực tiếp tuyến (Morgenstern & Price, 1965; Fredlund & Krahn, 1977): Cho rằng cần áp đặt một hàm số $f(x)$ liên tục để khống chế tỷ số giữa lực cắt và lực pháp tuyến giữa các thỏi.
  • Trường phái ma trận lực độc lập (Chugh, 1982; Zhu, Lee & Jiang, 2005): Cho rằng việc áp đặt hàm giả định nhân tạo có thể gây ra hiện tượng mất hội tụ số học hoặc xuất hiện ứng suất kéo phi thực tế tại các vùng biên của khối trượt.

Luận án của Soukha Yakoshi định vị trực tiếp vào dòng nghiên cứu của Chugh (1982) và phát triển vượt bậc bằng cách thiết lập quan hệ trực tiếp giữa lực cắt và lực pháp tuyến thông qua thông số hướng lực ($k, m$) và hệ số an toàn cục bộ $F_{si}$. Bằng cách đưa điều kiện phá hoại trượt tại đáy $T_i = (k \cdot N_i \tan\phi_1 + c_1 S_i B)/F_s$ và trên mặt bên $V_i = (m_i \cdot H_i \tan\phi_1 + m_i c_1 R_i B)/F_{si}$ vào ma trận toàn cục, nghiên cứu loại bỏ sự phụ thuộc vào các hàm giả định phân bố lực tùy tiện.

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

  1. Nghiên cứu của Duncan & Wright (2005): Sử dụng phương pháp Spencer và Morgenstern-Price để tính toán ổn định mố cầu nhưng chủ yếu áp dụng cho mặt trượt hình cung tròn và elip với tải trọng gộp, chưa phân rã chi tiết từng thành phần tĩnh tải kết cấu mố ($Q_{M1}$ đến $Q_{M5}$) tương tác trực tiếp lên khối trượt $M$.
  2. Nghiên cứu của Griffiths & Lane (1999) và Zhu et al. (2005): Sử dụng phần tử hữu hạn (FEM/SRM) và phương pháp tối ưu hóa không phụ thuộc lưới. Công trình của Soukha Yakoshi tạo ra một cầu nối xuất sắc khi đưa phương pháp giải tích - ma trận GLEM đạt tới độ chính xác tương đương FEM nhưng sở hữu tốc độ hội tụ nhanh hơn nhiều lần, cho phép phân tích hàng loạt kịch bản tải trọng LRFD phức tạp trên công cụ bảng tính tiêu chuẩn.

Đó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 lý thuyết Cân bằng Giới hạn Tổng quát (GLEM) từ bài toán mái dốc đất đơn thuần sang bài toán tương tác phức hợp Địa kỹ thuật - Kết cấu mố cầu - Móng nông. Đóng góp lý thuyết bao gồm việc mô hình hóa khối mố cầu như một thỏi trượt đặc biệt (Khối $M$) có trọng lượng $W_M$, tiếp nhận toàn bộ tổ hợp tải trọng ngoài $Q_M$, tải trọng đất đắp sau lưng $Q_{soil}$, hoạt tải xe chạy và lực quán tính động đất.

Nghiên cứu thiết lập hệ thống phương trình cân bằng hình chiếu tĩnh học cho thỏi thứ $i$: Chiếu lên phương pháp tuyến của mặt đáy khối: $$N_i - W_i \cos\alpha_i - H_i \cos(\theta_i - \alpha_i) + V_i \sin(\theta_i - \alpha_i) + H_{i+1} \cos(\theta_{i+1} - \alpha_i) - V_{i+1} \sin(\theta_{i+1} - \alpha_i) = 0$$

Chiếu lên phương tiếp tuyến của mặt đáy khối: $$T_i - W_i \sin\alpha_i + H_i \sin(\theta_i - \alpha_i) + V_i \cos(\theta_i - \alpha_i) - H_{i+1} \sin(\theta_{i+1} - \alpha_i) - V_{i+1} \cos(\theta_{i+1} - \alpha_i) = 0$$

Chuyển đổi thành dạng tường minh để thiết lập ma trận: $$-N_i + H_i \cos(\theta_i - \alpha_i) - V_i \sin(\theta_i - \alpha_i) - H_{i+1} \cos(\theta_{i+1} - \alpha_i) + V_{i+1} \sin(\theta_{i+1} - \alpha_i) = -W_i \cos\alpha_i \quad (1)$$ $$-T_i - H_i \sin(\theta_i - \alpha_i) - V_i \cos(\theta_i - \alpha_i) + H_{i+1} \sin(\theta_{i+1} - \alpha_i) + V_{i+1} \cos(\theta_{i+1} - \alpha_i) = -W_i \sin\alpha_i \quad (2)$$

Hệ thống điều kiện biên tiếp xúc tại đáy và mặt bên: $$T_i - N_i \frac{k_i \tan\phi_i}{F_s} = \frac{c_i S_i B}{F_s} \quad (7)$$ $$V_i - H_i \frac{m_i \tan\phi_i}{F_{si}} = \frac{m_i c_i R_i B}{F_{si}} \quad (8)$$

Hệ phương trình (1)-(15) được đóng gói hoàn chỉnh thành mô hình ma trận đại số tuyến tính $[A]{K} = {B}$ kích thước $(4n \times 4n)$, với véc-tơ ẩn số ${K} = [N_1, T_1, H_2, V_2, \dots, N_n, T_n, H_{n+1}, V_{n+1}]^T$.

       +-------------------------------------------------------------+
       |             MA TRẬN HỆ SỐ CÂN BẰNG TỔNG QUÁT [A]            |
       |  [ -1    0  -cos(d_1)  sin(d_1)   0    0   ...   0    0  ]  |
       |  [  0   -1   sin(d_1)  cos(d_1)   0    0   ...   0    0  ]  |
       |  [-a1    1     0          0       0    0   ...   0    0  ]  |
       |  [  0    0   -b2          1       0    0   ...   0    0  ]  |
       |  [  0    0   cos(d_2) -sin(d_2)  -1    0   ...   0    0  ]  |
       +-------------------------------------------------------------+
                                      *
       +-------------------------------------------------------------+
       |                    VÉC-TƠ ẨN SỐ LỰC {K}                     |
       |    {K} = [ N_1, T_1, H_2, V_2, N_2, T_2, ..., H_7, V_7 ]^T   |
       +-------------------------------------------------------------+
                                      =
       +-------------------------------------------------------------+
       |                    VÉC-TƠ TẢI TRỌNG {B}                     |
       |  {B} = [ -W_1*cos(a_1), -W_1*sin(a_1), c1*S1*B/Fs, ... ]^T   |
       +-------------------------------------------------------------+

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

Khung phân tích của luận án kết hợp đồng thời ba cấu phần:

  1. Lý thuyết phân tố lát cắt GLEM: Tự do hóa hình học mặt trượt thông qua tập hợp tọa độ đỉnh $(X_i, Y_i)$, bán kính/chiều dài mặt phân chia $R_i$, góc nghiêng đáy $\alpha_i$ và góc nghiêng mặt bên $\theta_i$.
  2. Lý thuyết trạng thái giới hạn LRFD: Tích hợp trực tiếp các hệ số tải trọng $\gamma_P$ vào các thành phần lực tác dụng lên khối mố $M$, bao gồm hệ số tải trọng tĩnh $\gamma_{P,QM} = 1.25$ và hệ số áp lực đất $\gamma_{P,Qsoil} = 1.35$.
  3. Thuật toán tối ưu hóa phi tuyến Solver (GRG - Generalized Reduced Gradient): Sử dụng hàm mục tiêu tối thiểu hóa hệ số an toàn: $$\min F_s = f(X_i, Y_i, \theta_i, R_i)$$ với các ràng buộc hình học không xuyên cắt, liên tục và điều kiện vật lý $F_{si} \ge F_s$.

Điều kiện biên xác định rõ ràng: Tại mặt trượt ban đầu $H_1 = 0, V_1 = 0$; tại điểm thoát trượt cuối cùng $H_{n+1} = 0, V_{n+1} = 0$. Hai trạng thái hệ số an toàn được xác lập mạch lạc:

  • $F_{s\min}$: Hệ số an toàn tới hạn khi mặt trượt giữa các thỏi đạt trạng thái cân bằng tuyệt đối ($FSI = \infty$).
  • $F_{s\text{med}}$: Hệ số an toàn trung bình khi giả định trạng thái ứng suất cắt giữa các thỏi đạt giá trị giới hạn thực tế ($FSI = 1.700$).

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi hệ hình nhận thức thực chứng duy lý (Positivist Paradigm) và lập trường bản thể luận khách quan (Objective Epistemology), sử dụng mô hình toán học giải tích kết hợp phương pháp số tối ưu hóa. Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng 1 - Kiểm chứng mô hình lý thuyết: Phân tích mái dốc đất dính đồng nhất trong hai trường hợp chuẩn tắc: trượt qua chân dốc và trượt sâu, nhằm đối chiếu với nghiệm giải tích kinh điển.
  • Tầng 2 - Phân tích trạng thái thi công: Khảo sát mố cầu móng nông chịu áp lực đất đắp với ba mức chiều cao $H_đ = 2.05\text{m}, 4.09\text{m}, 8.19\text{m}$; kết hợp tải trọng cơ giới thi công (xe 2 trục, xe 3 trục, xe lu rung) và tải trọng động đất giả tĩnh ($k_h, k_v$).
  • Tầng 3 - Phân tích trạng thái khai thác: Đánh giá ổn định trượt sâu mố cầu dưới tác dụng của tổ hợp hoạt tải xe thiết kế, tác động thủy văn dòng thấm (áp lực đẩy nổi $U, U'$, trọng lượng riêng bão hòa $W'_i$ và không bão hòa $W''_i$) và tải trọng động đất theo các Trạng thái Giới hạn Cường độ II và Cường độ III.
+---------------------------------------------------------------------------------------------------+
|                              QUY TRÌNH PHÂN TÍCH VÀ THIẾT KẾ NGHIÊN CỨU                           |
+---------------------------------------------------------------------------------------------------+
|                                                                                                   |
|  [BƯỚC 1: KHỞI TẠO HÌNH HỌC]                                                                      |
|  Thiết lập tọa độ mái đất (XPi, YPi), mặt trượt giả định (Xi, Yi), góc đáy alpha_i, mặt bên theta_i|
|                                                |                                                  |
|                                                v                                                  |
|  [BƯỚC 2: XÁC ĐỊNH TẢI TRỌNG VÀ TỔ HỢP LRFD]                                                      |
|  - Tĩnh tải mố: QM = QM1 + QM2 + QM3 + QM4 + QM5                                                 |
|  - Tải trọng đất đắp: Qsoil (theo Hđ = 2.05m, 4.09m, 8.19m)                                      |
|  - Hoạt tải & Động lực: QH, BR, IM, LS, PL, WL, WS, Gia tốc chấn động v, h                      |
|                                                |                                                  |
|                                                v                                                  |
|  [BƯỚC 3: THÀNH LẬP PHƯƠNG TRÌNH MA TRẬN GLEM]                                                    |
|  Xây dựng ma trận hệ số [A] kích thước 24x24 hoặc 28x28 và véc-tơ tải trọng {B}                  |
|                                                |                                                  |
|                                                v                                                  |
|  [BƯỚC 4: GIẢI HỆ VÀ TỐI ƯU HÓA HÌNH HỌC]                                                         |
|  Sử dụng thuật toán phi tuyến tìm min Fs qua việc biến thiên tọa độ mặt trượt                     |
|  Xuất kết quả: Fs_min (tại FSI = vô cùng) và Fs_med (tại FSI = 1.700)                             |
|                                                                                                   |
+---------------------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình tính toán phân rã miền trượt thành $n = 6$ đến $7$ thỏi phân tố. Đối với bài toán mố cầu, khối mố và khối đất đắp trong lòng mố được gộp thành phần tử tập trung có diện tích $A_M$, khối lượng $M_M$, chịu tác dụng của hệ lực phức hợp:

Ký hiệu tải trọng Đơn vị Giá trị thành phần tiêu chuẩn Hệ số tải trọng $\gamma_P$ Tải trọng tính toán (kN)
Bệ mố ($Q_{M1}$) kN 2829.75 1.25 3537.19
Tường thân ($Q_{M2}$) kN 2263.80 1.25 2829.75
Tường cánh ($Q_{M3}$) kN 902.40 1.25 1128.00
Tường đỉnh ($Q_{M4}$) kN 294.43 1.25 368.04
Vai kê ($Q_{M5}$) kN 36.38 1.25 45.48
Tổng tĩnh tải mố ($Q_M$) kN 7908.45 1.25 7908.45 (gộp hệ số: 9908.46)
Đất đắp trong lòng mố ($Q_{soil}$) kN 6096.87 ($H_đ=2.05\text{m}$) 1.35 8230.77
Đất đắp trong lòng mố ($Q_{soil}$) kN 3223.94 ($H_đ=4.09\text{m}$) 1.35 4352.31

Quy trình giải bài toán ma trận lặp được bảo đảm tính hội tụ nhờ thuật toán tối ưu hóa phi tuyến Solver, kiểm soát sai số lực khép kín ở mức $\varepsilon < 10^{-6}\text{ kN}$.

Data và phân tích

Dữ liệu hình học và cơ lý của bài toán kiểm chứng mái đất dính đồng nhất (trường hợp trượt qua chân dốc và trượt sâu) được chi tiết hóa trong các bảng tính ma trận:

Trọng lượng riêng của đất $\gamma = 20.0\text{ kN/m}^3$, gia tốc trọng trường $g = 10\text{ m/s}^2$, trọng lượng riêng bê tông $\gamma_c = 25.0\text{ kN/m}^3$. Các thông số hình học và tải trọng phân tố thỏi sau khi tối ưu hóa được định lượng chính xác:

Bảng thông số phân tố sau tối ưu hóa - Mái dốc đồng nhất trượt qua chân dốc ($F_{s\min} = 1.380$):

  • Thỏi 1: $\theta_1 = 33.69^\circ$, $R_1 = 2.42\text{m}$, $X_1 = -9.15$, $Y_1 = -6.10$, $\alpha_1 = -10.71^\circ$, $A_1 = 1.24\text{m}^2$, $W_1 = 24.74\text{ kN}$, $S_1 = 1.46\text{m}$.
  • Thỏi 2: $\theta_2 = 70.47^\circ$, $R_2 = 1.71\text{m}$, $X_2 = -7.71$, $Y_2 = -6.37$, $\alpha_2 = 10.81^\circ$, $A_2 = 3.74\text{m}^2$, $W_2 = 74.71\text{ kN}$, $S_2 = 2.23\text{m}$.
  • Thỏi 3: $\theta_3 = 84.74^\circ$, $R_3 = 2.43\text{m}$, $X_3 = -5.52$, $Y_3 = -5.95$, $\alpha_3 = 28.04^\circ$, $A_3 = 6.66\text{m}^2$, $W_3 = 133.12\text{ kN}$, $S_3 = 3.65\text{m}$.
  • Thỏi 4: $\theta_4 = 110.77^\circ$, $R_4 = 2.31\text{m}$, $X_4 = -2.30$, $Y_4 = -4.24$, $\alpha_4 = 41.55^\circ$, $A_4 = 5.31\text{m}^2$, $W_4 = 106.17\text{ kN}$, $S_4 = 3.00\text{m}$.
  • Thỏi 5: $\theta_5 = 130.26^\circ$, $R_5 = 1.85\text{m}$, $X_5 = -0.05$, $Y_5 = -2.25$, $\alpha_5 = 63.98^\circ$, $A_5 = 2.31\text{m}^2$, $W_5 = 46.22\text{ kN}$, $S_5 = 1.98\text{m}$.
  • Thỏi 6: $\theta_6 = 150.38^\circ$, $R_6 = 0.94\text{m}$, $X_6 = 0.82$, $Y_6 = -0.46$, $\alpha_6 = 90.00^\circ$, $A_6 = 0.19\text{m}^2$, $W_6 = 3.79\text{ kN}$, $S_6 = 0.46\text{m}$.

Bảng thông số phân tố sau tối ưu hóa - Mái dốc đồng nhất trượt sâu ($F_{s\min} = 1.469$):

  • Thỏi 1: $\theta_1 = 0.00^\circ$, $R_1 = 1.85\text{m}$, $X_1 = -11.00$, $Y_1 = -6.10$, $\alpha_1 = -20.06^\circ$, $A_1 = 0.52\text{m}^2$, $W_1 = 10.35\text{ kN}$, $S_1 = 1.63\text{m}$.
  • Thỏi 2: $\theta_2 = 60.45^\circ$, $R_2 = 0.64\text{m}$, $X_2 = -9.47$, $Y_2 = -6.66$, $\alpha_2 = -0.09^\circ$, $A_2 = 2.74\text{m}^2$, $W_2 = 54.76\text{ kN}$, $S_2 = 2.51\text{m}$.
  • Thỏi 3: $\theta_3 = 95.46^\circ$, $R_3 = 1.91\text{m}$, $X_3 = -6.96$, $Y_3 = -6.66$, $\alpha_3 = 15.72^\circ$, $A_3 = 5.29\text{m}^2$, $W_3 = 105.73\text{ kN}$, $S_3 = 2.67\text{m}$.
  • Thỏi 4: $\theta_4 = 110.76^\circ$, $R_4 = 2.57\text{m}$, $X_4 = -4.39$, $Y_4 = -5.94$, $\alpha_4 = 29.92^\circ$, $A_4 = 8.06\text{m}^2$, $W_4 = 161.20\text{ kN}$, $S_4 = 3.52\text{m}$.
  • Thỏi 5: $\theta_5 = 130.25^\circ$, $R_5 = 2.76\text{m}$, $X_5 = -1.34$, $Y_5 = -4.18$, $\alpha_5 = 49.96^\circ$, $A_5 = 8.78\text{m}^2$, $W_5 = 175.65\text{ kN}$, $S_5 = 4.37\text{m}$.
  • Thỏi 6: $\theta_6 = 150.38^\circ$, $R_6 = 1.70\text{m}$, $X_6 = 1.47$, $Y_6 = -0.84$, $\alpha_6 = 49.96^\circ$, $A_6 = 0.91\text{m}^2$, $W_6 = 18.26\text{ kN}$, $S_6 = 1.09\text{m}$.

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

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

Quá trình mô phỏng số học và tối ưu hóa ma trận GLEM đã dẫn đến 5 phát hiện mang tính bước ngoặt:

  1. Sự phân kỳ rõ rệt giữa $F_{s\min}$ và $F_{s\text{med}}$ phản ánh mức độ huy động sức kháng cắt:

    • Đối với mái đất đồng nhất trượt qua chân dốc: Hệ số an toàn tối thiểu đạt $F_{s\min} = 1.380$ (khi $FSI = \infty$), trong khi hệ số an toàn trung bình đạt $F_{s\text{med}} = 1.510$ (khi $FSI = 1.700$), chênh lệch $9.42%$.
    • Đối với trường hợp trượt sâu: $F_{s\min} = 1.469$, trong khi $F_{s\text{med}} = 1.654$, chênh lệch $12.59%$. Điều này chứng minh rằng việc áp đặt độ bền trượt giữa các thỏi ở mức thực tế ($FSI = 1.7$) giúp giải phóng ứng suất cục bộ, làm tăng hệ số an toàn tính toán của công trình.

  2. Quy luật dịch chuyển mặt trượt nguy hiểm theo chiều cao đất đắp $H_đ$:

    • Khi chiều cao đất đắp sau mố tăng từ $H_đ = 2.05\text{m}$ lên $H_đ = 4.09\text{m}$ và đạt cực đại $H_đ = 8.19\text{m}$, tổng tải trọng thẳng đứng tính toán tăng từ $12260.76\text{ kN}$ lên $16139.22\text{ kN}$.
    • Áp lực chủ động gia tăng mạnh đẩy cung trượt ăn sâu vào tầng đất tự nhiên bên dưới bệ móng, làm tăng bán kính trượt và kéo dài đáy cung trượt từ $X = -15.19\text{m}$ ($H_đ = 4.09\text{m}$) ra xa đến $X = -16.69\text{m}$ ($H_đ = 2.05\text{m}$ với tải trọng bản thân lớn).

  3. Hiện tượng bất đối xứng ứng suất tại khối mố $M$: Khối trượt $M$ chứa mố cầu gánh chịu lực nén pháp tuyến đáy cực lớn (ví dụ trong trường hợp $H_đ = 4.09\text{m}$, lực pháp tuyến $N_4$ và lực tiếp tuyến $T_4$ đạt giá trị áp đảo trong véc-tơ ${B}$ với $-(W_M + Q_M)\cos\alpha_M = -17272\text{ kN}$ và $-(W_M + Q_M)\sin\alpha_M = -5412.5\text{ kN}$). Sự tập trung ứng suất này tạo ra hiệu ứng "nêm chèn", làm thay đổi toàn bộ hướng của véc-tơ phản lực giữa các thỏi lân cận $H_4, V_4$ và $H_5, V_5$.

  4. Độ nhạy cảm cao đối với tải trọng thi công động và động đất: Hoạt tải xe 3 trục và xe lu rung trong giai đoạn thi công gây suy giảm $F_s$ từ $15%$ đến $22%$ so với trạng thái tĩnh. Khi có lực quán tính động đất ($v, h$), mặt trượt có xu hướng phẳng hóa ở phần đáy và mở rộng về phía hạ lưu mái dốc.

  5. Hiệu ứng giảm áp lực hữu hiệu do mực nước ngầm: Áp lực đẩy nổi $U, U'$ làm triệt tiêu đáng kể thành phần lực pháp tuyến hữu hiệu $N'_i = N_i - U_i$, dẫn đến sức kháng cắt tiếp tuyến $T_i$ giảm mạnh, đưa hệ số an toàn trong giai đoạn khai thác về sát ngưỡng giới hạn ổn định ($F_s \approx 1.05 - 1.15$).

+---------------------------------------------------------------------------------------------------+
|                        SO SÁNH HỆ SỐ AN TOÀN VÀ TỔ HỢP TẢI TRỌNG CỐT LÕI                          |
+---------------------------------------------------------------------------------------------------+
| Kịch bản phân tích                    | Fs_min (FSI = vô cùng) | Fs_med (FSI = 1.700) | Tải trọng (kN)|
|---------------------------------------|------------------------|----------------------|---------------|
| Mái dốc đồng nhất (Trượt chân dốc)    |         1.380          |        1.510         |       --      |
| Mái dốc đồng nhất (Trượt sâu)         |         1.469          |        1.654         |       --      |
| Mố móng nông (Thi công, Hđ = 2.05m)   |      Tối ưu hóa        |      Tối ưu hóa      |   16139.22    |
| Mố móng nông (Thi công, Hđ = 4.09m)   |      Tối ưu hóa        |      Tối ưu hóa      |   12260.76    |
| Mố móng nông (Thi công, Hđ = 8.19m)   |      Tối ưu hóa        |      Tối ưu hóa      |   16139.22    |
+---------------------------------------------------------------------------------------------------+

Implications đa chiều

  • Về mặt học thuật: Luận án cung cấp một bộ khung thuật toán ma trận mẫu mực cho bài toán tương tác phức hợp giữa kết cấu cầu và khối trượt địa kỹ thuật.
  • Về mặt phương pháp luận: Chứng minh tính khả thi của việc tích hợp giải thuật tối ưu hóa phi tuyến Solver vào hệ ma trận đại số GLEM, thay thế hoàn toàn các phần mềm thương mại đóng kín (black-box software) như GeoStudio/SLOPE/W hay PLAXIS trong công tác thẩm tra độc lập.
  • Về mặt thực tiễn kỹ thuật: Đưa ra khuyến cáo thiết kế quan trọng: Đối với mố cầu trên móng nông có chiều cao đất đắp $H_đ > 4.0\text{m}$, không được phép chỉ kiểm tra trượt phẳng đáy móng đơn thuần theo quy trình cơ học kết cấu cổ điển mà bắt buộc phải kiểm toán trượt sâu khối trượt tổng thể đi qua nền đất dưới bệ móng.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án có ba giới hạn khoa học cần được nhìn nhận khách quan:

  1. Giới hạn bài toán biến dạng phẳng 2D: Nghiên cứu giả định điều kiện bài toán biến dạng phẳng ($B = 1\text{m}$ bề rộng), chưa xét đến hiệu ứng không gian ba chiều (3D end effects) và ma sát thành bên của mố cầu có cánh chữ U.
  2. Mô hình vật liệu đàn - dẻo lý tưởng: Tiêu chuẩn Mohr-Coulomb được giả định đẳng hướng, chưa phản ánh hiện tượng suy giảm độ bền theo thời gian (softening behavior) và từ biến của đất sét yếu dưới đáy móng.
  3. Phương pháp động đất giả tĩnh: Việc sử dụng hệ số gia tốc quán tính tương đương ($k_h, k_v$) chưa mô phỏng được sự tích tụ áp lực nước lỗ rỗng động và hiện tượng hóa lỏng cục bộ của cát bão hòa trong động đất.

Chương trình nghiên cứu mở rộng 5-10 năm tới bao gồm:

  • Mở rộng thuật toán ma trận GLEM sang không gian 3 chiều (3D GLEM Matrix Formulation) cho khối trượt mố cầu phức tạp.
  • Tích hợp mô hình trạng thái tới hạn (Critical State Soil Mechanics - Cam-Clay) vào điều kiện trượt tiếp xúc đáy và mặt bên thỏi phân tố.
  • Phát triển module phân tích động phi tuyến theo lịch sử thời gian (Time-history dynamic response) kết hợp lan truyền sóng chấn động qua nền đất móng nông.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning / Deep Neural Operators) để xấp xỉ tức thời mặt trượt tới hạn dựa trên tập dữ liệu ma trận GLEM đã huấn luyện.

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

Nghiên cứu của Soukha Yakoshi tạo ra sức ảnh hưởng sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Định hình lại phương pháp giảng dạy và nghiên cứu chuyên sâu về ổn định công trình ngầm và mố móng cầu tại các trường đại học kỹ thuật. Các bài báo khoa học phát triển từ mô hình ma trận GLEM của luận án có tiềm năng trích dẫn cao trong các tạp chí ISI/Scopus chuyên ngành như Computers and Geotechnics, Geotechnique, hay ASCE Journal of Geotechnical and Geoenvironmental Engineering.
  • Chuyển đổi công nghiệp xây dựng: Cung cấp công cụ tính toán mã nguồn mở, minh bạch trên Excel/VBA cho các tập đoàn tư vấn thiết kế giao thông, giúp tối ưu hóa kích thước bệ mố nông, giảm thiểu chi phí xử lý nền móng từ $15%$ đến $30%$ so với phương án cọc khoan nhồi trong điều kiện địa chất phù hợp.
  • Tác động chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học định lượng để bổ sung, sửa đổi các quy định kiểm toán ổn định tổng thể mố cầu trong Tiêu chuẩn Thiết kế Cầu quốc gia TCVN 11823 và AASHTO LRFD.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận hệ phương pháp toán học ma trận chuẩn xác, tường minh từ khâu thiết lập phương trình cân bằng hình chiếu lực đến thuật toán giải số tối ưu.
  • Kỹ sư thiết kế Cầu - Đường và Địa kỹ thuật: Nắm bắt công cụ lập trình tính toán linh hoạt để kiểm soát an toàn mố móng nông dưới các tổ hợp tải trọng bất lợi nhất (thi công, khai thác, ngập nước, động đất).
  • Các cơ quan quản lý và Chủ đầu tư dự án hạ tầng: Có cơ sở khoa học tin cậy để thẩm tra biện pháp thi công, đánh giá rủi ro sạt trượt và nghiệm thu chất lượng công trình giao thông trọng điểm.

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ì?

Luận án đã mở rộng thành công lý thuyết Cân bằng Giới hạn Tổng quát (GLEM) sang hệ tương tác kết cấu - nền đất bằng cách thiết lập ma trận toàn cục $[A]{K} = {B}$ giải đồng thời $4n$ ẩn số lực pháp tuyến và lực tiếp tuyến ($N_i, T_i, H_i, V_i$). Điểm độc đáo nhất là mô hình hóa khối mố $M$ thành phần tử tương tác đặc biệt tiếp nhận đầy đủ tĩnh tải kết cấu ($\sum Q_M = 7908.45\text{ kN}$), tải trọng đất đắp lòng mố ($Q_{soil} = 6096.87\text{ kN}$), lực hãm xe ($BR$), xung kích ($IM$) và động đất ($v, h$) mà không cần bất kỳ giả định nhân tạo nào về hàm phân bố lực cắt giữa các thỏi.

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

So với phương pháp Spencer (1967) áp đặt góc nghiêng lực liên thỏi không đổi và phương pháp Morgenstern-Price (1965) sử dụng hàm $f(x)$ giả định, mô hình của luận án thiết lập điều kiện bền trượt thực tế độc lập tại từng mặt phân chia $V_i = (m_i H_i \tan\phi + m_i c R_i B)/F_{si}$. Quy trình này kết hợp trực tiếp công cụ tối ưu hóa phi tuyến Solver GRG để tự động biến thiên tọa độ mặt trượt $(X_i, Y_i)$, đạt độ chính xác tương đương phương pháp phần tử hữu hạn (FEM) nhưng minh bạch và nhanh chóng hơn.

3. Phát hiện bất ngờ và có ý nghĩa thực tiễn lớn nhất từ dữ liệu?

Phát hiện cho thấy hệ số an toàn trượt sâu $F_s$ không tỷ lệ nghịch đơn thuần với tải trọng thẳng đứng. Khi chiều cao đất đắp tăng từ $H_đ = 2.05\text{m} \to 4.09\text{m} \to 8.19\text{m}$, dù tổng tĩnh tải tăng mạnh lên $16139.22\text{ kN}$, nhưng do sự dịch chuyển của tâm xoay và sự thay đổi góc dốc mặt trượt đáy $\alpha_M$, ứng suất nén hữu hiệu dưới đáy mố tăng lên làm tăng sức kháng cắt ma sát $N_M \tan\phi$. Điều này chứng minh rằng việc đánh giá ổn định mố móng nông phải dựa trên sự cân bằng động lực học giữa mô-men gây trượt và sức kháng ma sát tăng cường dưới đáy bệ móng.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) minh bạch không?

Có. Luận án cung cấp phụ lục chương trình tính toán cực kỳ chi tiết với cấu trúc ma trận số học $[A]$, véc-tơ tải trọng ${B}$, bảng thông số tọa độ $(X_i, Y_i, XP_i, YP_i)$, diện tích $A_i$, trọng lượng $W_i$, chiều dài đáy $S_i$, chiều dài vách $R_i$ và các hệ số chỉ hướng $k, m$. Mọi nhà nghiên cứu đều có thể tái lập hoàn toàn thuật toán trên Microsoft Excel, MATLAB hoặc Python.

5. Lộ trình nghiên cứu 10 năm tới được phác thảo ra sao?

Lộ trình hướng tới việc số hóa toàn diện: (1) Phát triển phần mềm GLEM-3D Bridge Abutment; (2) Tích hợp quan trắc áp lực lỗ rỗng thời gian thực (BIM-IoT) vào mô hình ma trận để cảnh báo sớm trượt mố cầu; (3) Hoàn thiện chỉ dẫn kỹ thuật thiết kế móng nông mố cầu chịu tải trọng động đất cho ngành giao thông vận tải.

Kết luận

  1. Luận án của Soukha Yakoshi (2023) đã hoàn thiện xuất sắc khung lý thuyết toán học ma trận cho Phương pháp Cân bằng Giới hạn Tổng quát (GLEM) trong phân tích ổn định mái dốc và trượt sâu mố cầu trên móng nông.
  2. Thiết lập thành công hệ phương trình đại số tuyến tính $[A]{K} = {B}$ kích thước lớn, giải chính xác phân bố nội lực $N_i, T_i, H_i, V_i$ trên từng thỏi trượt và khối mố kết cấu.
  3. Tích hợp hiệu quả thuật toán tối ưu hóa phi tuyến để xác định tự động mặt trượt nguy hiểm nhất có dạng hình học bất kỳ, lượng hóa chính xác $F_{s\min} = 1.380$ và $F_{s\text{med}} = 1.510$ (trượt chân dốc) cũng như $F_{s\min} = 1.469$ và $F_{s\text{med}} = 1.654$ (trượt sâu).
  4. Định lượng toàn diện tác động của 18 kịch bản tải trọng phức hợp bao gồm biến thiên chiều cao đất đắp ($H_đ = 2.05\text{m}, 4.09\text{m}, 8.19\text{m}$), hoạt tải cơ giới thi công (xe 2 trục, 3 trục, xe lu), mực nước ngầm bão hòa và lực chấn động động đất theo tiêu chuẩn LRFD.
  5. Cung cấp quy trình tính toán minh bạch, khả thi cao trên các công cụ tính toán phổ thông, mở ra hướng ứng dụng thực tiễn to lớn cho công tác tư vấn thiết kế và thẩm tra an toàn công trình giao thông.
  6. Đặt nền móng vững chắc cho các nghiên cứu tiếp theo về ổn định không gian 3 chiều và cơ chế phá hoại động lực học của hệ thống công trình mố trụ cầu trong tương lai.