Tổng quan về luận án

Tính ổn định tổng thể của công trình mố cầu trên móng nông đặt trên sườn dốc hoặc nền đất đắp là bài toán cơ học đất - địa kỹ thuật công trình phức tạp, đòi hỏi sự kết hợp đồng thời giữa cơ học môi trường rời và tương tác kết cấu - nền đất (Soil-Structure Interaction - SSI). Trong thực tiễn kỹ thuật cầu đường, sự phá hoại ổn định trượt sâu (deep-seated sliding failure) của cụm công trình mố cầu thường dẫn đến các biến dạng trượt trồi, lún sụt nghiêm trọng mố trụ và làm gián đoạn khai thác tuyến giao thông huyết mạch.

Khoảng trống nghiên cứu then chốt (specific research gap) nằm ở chỗ: các phương pháp phân chia thỏi cổ điển (Bishop, Janbu, Spencer, Morgenstern-Price) thường dựa trên các giả thiết hình học mặt trượt cố định (cung tròn hoặc đa giác đơn giản) và chưa mô hình hóa đầy đủ ma trận cân bằng tĩnh học tương tác đa diện giữa khối mố cứng có tải trọng tập trung cao với các khối đất đắp sau lưng mố chịu áp lực chủ động/bị động biến thiên. Hơn nữa, việc xác định đồng thời hệ số an toàn nhỏ nhất ($F_{s,\min}$) và quỹ đạo mặt trượt nguy hiểm nhất khi thay đổi đồng thời chiều cao đất đắp ($H_đ$), mực nước ngầm, hoạt tải xe cơ giới (xe 3 trục, xe 2 trục, xe lu) và tải trọng động đất theo các trạng thái giới hạn (TTGH) cường độ vẫn còn thiếu một công cụ phân tích ma trận tường minh có thể tự động hóa tối ưu hóa phi tuyến.

Nghiên cứu giải quyết hệ thống câu hỏi và giả thuyết khoa học được chuẩn hóa:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập hệ phương trình ma trận cân bằng tĩnh học lực pháp tuyến ($N_i$), lực tiếp tuyến ($T_i$) tại đáy khối trượt và lực tương tác giữa các phân tố ($H_i, V_i$) theo Phương pháp Cân bằng Giới hạn Tổng quát (GLEM - General Limit Equilibrium Method) khi có sự hiện diện của khối kết cấu mố cầu ($W_M, Q_M$)?
  • Câu hỏi nghiên cứu 2 (RQ2): Ảnh hưởng của các tham số điều kiện biên (chiều cao đất đắp $H_đ = 2{,}05\text{ m}; 4{,}09\text{ m}; 8{,}19\text{ m}$, hoạt tải xe lu, xe tải 2-3 trục, mực nước ngầm và lực quán tính động đất $v, h$) đến hệ số an toàn trượt sâu $F_s$ của mố cầu trên móng nông diễn ra theo quy luật giải tích nào?
  • Giả thuyết nghiên cứu 1 (H1): Việc chuyển đổi hệ phương trình vi tích phân phi tuyến cân bằng giới hạn thành hệ phương trình đại số ma trận dạng $[A]{K} = {B}$ kích thước $24 \times 24$ cho phép giải chính xác đồng thời các ứng suất tiếp xúc đáy và mặt phân chia phân tố mà không cần giả định trước hàm phân bố lực liên thỏi.
  • Giả thuyết nghiên cứu 2 (H2): Bằng cách áp dụng thuật toán tối ưu hóa phi tuyến trên nền tảng hàm Solver để hiệu chỉnh tọa độ $(X_i, Y_i)$ và góc nghiêng đáy $(\alpha_i)$, mặt trượt nguy hiểm nhất sẽ hội tụ về cấu hình có hệ số an toàn thực nhỏ hơn đáng kể so với mặt trượt tròn giả định ban đầu.

Khung lý thuyết của công trình tích hợp Lý thuyết Trạng thái Giới hạn Cơ học Đất Mohr-Coulomb, Lý thuyết Cân bằng Giới hạn Tổng quát GLEM, 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. Đóng góp đột phá của nghiên cứu được định lượng qua việc giải quyết bài toán ma trận tương tác hệ $24 \times 24$, tối ưu hóa hệ số an toàn mái dốc đồng nhất từ trạng thái ban đầu về $F_{s,\min} = 1{,}380$ (trượt chân dốc) và $F_{s,\min} = 1{,}469$ (trượt sâu), cũng như đánh giá mố cầu với tổng tĩnh tải tính toán giai đoạn thi công lên tới $16.139{,}22\text{ kN}$ ($H_đ = 2{,}05\text{ m}$) và $12.260{,}76\text{ kN}$ ($H_đ = 4{,}09\text{ m}$).

Literature Review và Positioning

Lịch sử phát triển các phương pháp tính toán ổn định mái dốc và kết cấu chắn giữ chứng kiến sự tiến hóa liên tục từ các mô hình đơn giản đến các hệ thống cân bằng giới hạn đa biến:

  • Trường phái Cân bằng Cổ điển (Limit Equilibrium Methods - LEM): Khởi nguồn từ phương pháp thỏi kinh điển của Fellenius (1936) và Bishop (1955) giản lược, chỉ xét cân bằng mô-men và giả định lực liên thỏi nằm ngang. Janbu (1957) mở rộng cho mặt trượt phi tròn dựa trên cân bằng lực nhưng hệ số điều chỉnh phụ thuộc kinh nghiệm. Morgenstern & Price (1965) cùng Spencer (1967) hoàn thiện mô hình cân bằng toàn phần (cân bằng đồng thời lực và mô-men) bằng cách đưa vào hàm liên thỏi giả định $f(x)$ và tham số vô hướng $\lambda$. Sarma (1973) phát triển phương pháp tính gia tốc động đất tới hạn cho mái dốc nhiều lớp.
  • Trường phái Cân bằng Tổng quát (GLEM) và Phân tích Ma trận: Fredlund & Krahn (1977) tổng quát hóa các phương pháp thỏi thành khung phân tích GLEM thống nhất. Tuy nhiên, việc áp dụng GLEM vào phân tích ổn định trượt sâu của mố cầu chịu tải trọng kết cấu phức tạp gặp trở ngại do sự mất cân bằng cục bộ tại ranh giới tiếp xúc giữa bê tông và đất. Zhu et al. (2003) và Chen et al. (2005) đưa ra các thuật toán tìm kiếm mặt trượt tới hạn nhưng phần lớn áp dụng cho mái dốc địa hình tự nhiên mà chưa mô đun hóa tổ hợp tải trọng mố cầu.
  • Trường phái Phần tử Hữu hạn Giảm Độ Bền (FEM-SRM): Griffiths & Lane (1999) chứng minh FEM-SRM không cần giả định mặt trượt nhưng gặp hạn chế lớn về thời gian tính toán và khó giải thích trực tiếp cơ chế truyền tải trọng cục bộ theo các tiêu chuẩn thiết kế cầu hiện hành.
       TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT ỔN ĐỊNH MÁI DỐC VÀ MỐ CẦU
 [Mặt trượt tròn]      [Bỏ qua lực đứng]  [Cân bằng toàn phần, giả định f(x)]
                       LUẬN ÁN HIỆN TẠI (2023): GLEM Ma Trận Hóa + Solver
                       [Tích hợp tương tác Mố - Nền - Tải trọng AASHTO LRFD]

Tranh luận học thuật cốt lõi tồn tại giữa quan điểm giả định hàm lực liên thỏi cố định (Spencer, Morgenstern-Price) và quan điểm tự do hóa ứng suất liên thỏi thông qua điều kiện bền cục bộ. Luận án định vị nghiên cứu tại điểm giao thoa giữa độ chính xác cơ học của GLEM và khả năng ứng dụng thực hành kỹ thuật cao thông qua kỹ thuật ma trận hóa tường minh. So sánh với các nghiên cứu quốc tế:

  1. So với nghiên cứu của Duncan & Wright (2005) trong sổ tay Soil Strength and Slope Stability (chỉ dừng lại ở phân tích mái đắp giao thông tiêu chuẩn), nghiên cứu này mở rộng phân tích khối mố phức hợp gồm 5 thành phần kết cấu ($Q_{M1}$ đến $Q_{M5}$) chịu tác động phối hợp của áp lực đất đắp $Q_{\text{soil}}$ và tải trọng khai thác.
  2. So với khung phân tích mố cầu của Cục Đường bộ Liên bang Hoa Kỳ (FHWA-NHI-06-088: Geotechnical Aspects of Pavements and Bridge Foundations), mô hình đề xuất trong luận án thiết lập liên kết đại số trực tiếp giữa trạng thái giới hạn cường độ (Strength Limit States) và ma trận phân tố trượt sâu, loại bỏ sự phụ thuộc vào các biểu đồ tra cứu đơn giản hóa.

Đó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 bằng cách tích hợp trực tiếp các điều kiện bền cắt trên mặt tiếp xúc giữa các khối trượt và mặt trượt đáy vào hệ phương trình giải tích tường minh.

                               MÔ HÌNH HỆ KHỐI TRƯỢT PHÂN TỐ
     (Lực tiếp tuyến/pháp      (Trọng lượng Wi)         (Lực tiếp tuyến/pháp
                             Ti (Lực cắt đáy khối)
                             Ni (Lực pháp tuyến đáy)
                             [Góc nghiêng đáy αi]

Mô hình thiết lập cân bằng lực trên từng phân tố thứ $i$:

  1. Phương trình cân bằng lực theo phương pháp tuyến mặt đáy phân tố: $$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$$ hay chuyển về dạng chuẩn tắc: $$-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$$

  2. Phương trình cân bằng lực theo phương tiếp tuyến mặt đáy phân tố: $$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$$ hay chuyển về dạng chuẩn tắc: $$-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$$

  3. Các điều kiện bền Mohr-Coulomb trên mặt đáy và mặt liên thỏi được tuyến tính hóa theo hệ số an toàn $F_s$ và $F_{si}$: $$T_i = \frac{N_i \tan \varphi_1 + c_1 S_i B}{F_s} \Longleftrightarrow -k_i \frac{\tan \varphi_1}{F_s} N_i + T_i = \frac{c_1 S_i B}{F_s}$$ $$V_i = \frac{m_i (H_i \tan \varphi_1 + c_1 R_i B)}{F_{si}} \Longleftrightarrow -m_i \frac{\tan \varphi_1}{F_{si}} H_i + V_i = \frac{m_i c_1 R_i B}{F_{si}}$$

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

Khung phân tích hợp nhất ba nền tảng lý thuyết cơ sở: (1) Lý thuyết Cơ học Môi trường Rời Cổ điển Mohr-Coulomb; (2) Lý thuyết Tối ưu hóa Phi tuyến Cổ điển với thuật toán Generalized Reduced Gradient (GRG2); và (3) Lý thuyết Thiết kế Kết cấu Nhịp và Mố Trụ Cầu theo Trạng thái Giới hạn (AASHTO LRFD).

Hệ ma trận phương trình tổng quát biểu diễn dưới dạng: $$[A]{24 \times 24} \cdot {K}{24 \times 1} = {B}_{24 \times 1}$$ Trong đó véc-tơ ẩn số chứa toàn bộ các nội lực tiếp xúc của 6 phân tố đất và 1 khối mố phức hợp: $${K} = [N_1, T_1, H_2, V_2, N_2, T_2, H_3, V_3, N_3, T_3, H_M, V_M, N_M, T_M, H_5, V_5, N_5, T_5, H_6, V_6, N_6, T_6, H_7, V_7]^T$$

Điều kiện biên được xác định nghiêm ngặt: $H_1 = 0$, $V_1 = 0$ tại đỉnh cung trượt và $H_7 = 0, V_7 = 0$ tại vị trí thoát của mặt trượt ngoài chân dốc.

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

Thiết kế nghiên cứu

Nghiên cứu tuân theo thế giới quan thực chứng (positivism) kết hợp mô hình hóa tất định (deterministic computational modeling). Hệ thống phân tích được thiết kế đa cấp độ (multi-level design):

  • Cấp độ 1 (Level 1): Khảo sát chuẩn hóa tính ổn định của mái dốc đất dính đồng nhất ($c, \varphi, \gamma$) với hai kịch bản trượt qua chân dốc và trượt sâu.
  • Cấp độ 2 (Level 2): Khảo sát tương tác kết cấu mố cầu móng nông trong giai đoạn thi công chịu áp lực đất đắp với các chiều cao thay đổi ($H_đ = 2{,}05\text{ m}; 4{,}09\text{ m}; 8{,}19\text{ m}$).
  • Cấp độ 3 (Level 3): Khảo sát tổng hợp trong giai đoạn thi công và khai thác khi chịu tác động đồng thời của hoạt tải giao thông trục nặng, xe lu rung nén tĩnh, áp lực thủy tĩnh nước ngầm và tải trọng động đất theo các TTGH Cường độ 1, 2 và 3.
                      SƠ ĐỒ THIẾT KẾ NGHIÊN CỨU ĐA CẤP ĐỘ

Quy trình nghiên cứu rigorous

Quy trình tính toán số học được kiểm soát sai số thông qua công cụ Microsoft Excel Solver với độ chính xác số dư ma trận đạt mức hội tụ cao:

  1. Thiết lập hình học ban đầu: Khởi tạo tọa độ nút của mái dốc $(XP_i, YP_i)$ và tọa độ ban đầu của mặt trượt $(X_i, Y_i)$. Phân chia hệ trượt thành 6 thỏi đất và 1 thỏi mố ($M$).
  2. Tính toán đặc trưng hình học và trọng lượng: Tính diện tích phân tố $A_i$, chiều dài đáy $S_i = \sqrt{\Delta X_i^2 + \Delta Y_i^2}$, chiều dài mặt phân chia $R_i$, góc nghiêng đáy $\alpha_i = \arctan(\Delta Y_i / \Delta X_i)$, và góc nghiêng mặt biên $\theta_i$. Trọng lượng phân tố được xác định: $W_i = \gamma \cdot A_i$.
  3. Tổ hợp tải trọng kết cấu mố cầu: Xác định chính xác các thành phần tĩnh tải mố:
    • Tĩnh tải bệ mố: $Q_{M1} = 2.829{,}75\text{ kN}$
    • Tĩnh tải tường thân: $Q_{M2} = 2.263{,}80\text{ kN}$
    • Tĩnh tải tường cánh: $Q_{M3} = 902{,}40\text{ kN}$
    • Tĩnh tải tường đỉnh: $Q_{M4} = 294{,}43\text{ kN}$
    • Tĩnh tải vai kê: $Q_{M5} = 36{,}38\text{ kN}$
    • Tổng tĩnh tải kết cấu mố: $W_M = 12.423{,}63\text{ kN}$ (khi $H_đ = 2{,}05\text{ m}$ và $8{,}19\text{ m}$) hoặc $9.550{,}69\text{ kN}$ (khi $H_đ = 4{,}09\text{ m}$).
    • Tĩnh tải đất đắp trong lòng mố: $Q_{\text{soil}} = 6.096{,}87\text{ kN}$ ($H_đ = 2{,}05\text{ m}$) hoặc $3.223{,}94\text{ kN}$ ($H_đ = 4{,}09\text{ m}$).
  4. Áp dụng hệ số tải trọng theo TTGH Cường độ I: Hệ số $\gamma_P = 1{,}25$ cho tĩnh tải kết cấu mố và $\gamma_P = 1{,}35$ cho đất đắp trong lòng mố, mang lại tổng tải trọng tính toán thẳng đứng $V = 16.139{,}22\text{ kN}$ đối với $H_đ = 2{,}05\text{ m}$ và $V = 12.260{,}76\text{ kN}$ đối với $H_đ = 4{,}09\text{ m}$.

Data và phân tích

Ma trận $[A]$ được lắp ghép tự động từ các hệ số lượng giác của góc nghiêng và tham số cường độ chống cắt. Dữ liệu trích xuất từ các bảng tính toán tối ưu hóa cho thấy cấu trúc dữ liệu tường minh:

Phân tố ($i$) Góc đáy $\alpha_i$ (rad) Góc $\alpha_i$ ($^\circ$) Chiều dài đáy $S_i$ (m) Diện tích $A_i$ (m$^2$) Trọng lượng $W_i$ (kN) Bán kính $R_i$ (m)
Thỏi 1 $-0{,}19$ $-10{,}71$ $1{,}46$ $1{,}24$ $24{,}74$ $2{,}42$
Thỏi 2 $0{,}19$ $10{,}81$ $2{,}23$ $3{,}74$ $74{,}71$ $1{,}71$
Thỏi 3 $0{,}49$ $28{,}04$ $3{,}65$ $6{,}66$ $133{,}12$ $2{,}43$
Thỏi 4 $0{,}73$ $41{,}55$ $3{,}00$ $5{,}31$ $106{,}17$ $2{,}31$
Thỏi 5 $1{,}12$ $63{,}98$ $1{,}98$ $2{,}31$ $46{,}22$ $1{,}85$
Thỏi 6 $1{,}57$ $90{,}00$ $0{,}46$ $0{,}19$ $3{,}79$ $0{,}94$

Bảng dữ liệu tối ưu hóa hình học phân tố cho mái dốc đồng nhất trượt qua chân dốc.

Quá trình kiểm tra độ vững chắc (robustness checks) được thực hiện bằng cách so sánh nghiệm của ma trận khi gán điều kiện an toàn liên thỏi $FSI = \infty$ (tương đương hệ cân bằng lực độc lập) và khi áp đặt $FSI = 1{,}700$ (tương đương trạng thái phân bố ứng suất tiếp hạn chế trên biên thỏi).

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

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

  1. Sự phụ thuộc của hệ số an toàn thực vào giả thiết liên thỏi: Đối với mái dốc đất dính đồng nhất trong trường hợp trượt qua chân dốc, khi giả thiết $FSI = \infty$, giá trị $F_{s,\min} = 1{,}380$. Khi thiết lập điều kiện ràng buộc liên thỏi $FSI(i) = 1{,}700$, hệ số an toàn đáy hội tụ về $FS = 1{,}510$ (tăng $9{,}42%$). Tương tự, đối với trường hợp trượt sâu, khi $FSI = \infty$, $F_{s,\min} = 1{,}469$; nhưng khi áp đặt $FSI = 1{,}700$, $FS = 1{,}654$ (tăng $12{,}59%$).
  2. Quy luật phi tuyến của ổn định trượt sâu theo chiều cao đất đắp: Khi chiều cao đất đắp sau mố tăng từ $H_đ = 2{,}05\text{ m}$ lên $H_đ = 4{,}09\text{ m}$ và $H_đ = 8{,}19\text{ m}$, lực đẩy ngang chủ động tác dụng vào khối mố tăng phi tuyến, làm dịch chuyển tâm cung trượt nguy hiểm sâu hơn vào tầng đất tự nhiên phía dưới móng mố, làm giảm đáng kể hệ số an toàn tổng thể trong giai đoạn thi công.
  3. Hiện tượng suy giảm hệ số an toàn dưới tổ hợp tải trọng thi công: Sự hiện diện của tải trọng xe lu rung ép và xe đắp 3 trục trực tiếp sau lưng mố tạo ra các xung lực cục bộ cực đại tại phân tố số 4 và phân tố mố $M$, làm gia tăng lực cắt tiếp xúc đáy $T_M$ lên mức $2.900{,}3\text{ kN}$ ($H_đ = 2{,}05\text{ m}$) và $5.412{,}5\text{ kN}$ ($H_đ = 4{,}09\text{ m}$).
  4. Ảnh hưởng bất lợi kép của mực nước ngầm và động đất: Dưới tác động của áp lực đẩy nổi $U, U'$ làm giảm ứng suất pháp hiệu quả $N_i'$ cùng lực quán tính nằm ngang $h = a_h \cdot M_i$, hệ số an toàn trượt sâu của mố cầu suy giảm nghiêm trọng nhất ở TTGH Cường độ 3 và TTGH Đặc biệt (Động đất), đòi hỏi phải mở rộng kích thước bệ mố $B$ hoặc xử lý gia cố nền đất dưới đáy móng nông.
       BIẾN THIÊN HỆ SỐ AN TOÀN THEO ĐIỀU KIỆN BIÊN LIÊN THỎI (FSI)
  FS 
                                FSI = ∞                       FSI = 1.7

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh rằng các phương pháp cổ điển bỏ qua lực cắt liên thỏi hoặc cố định tỷ số lực liên thỏi có thể đánh giá thừa hoặc thiếu độ ổn định của hệ thống mố - nền từ $9{,}4%$ đến $12{,}6%$.
  • Về mặt phương pháp luận: Khẳng định tính khả thi và độ chính xác của giải pháp ma trận hóa $24 \times 24$ kết hợp thuật toán Solver, mở ra hướng ứng dụng tự động hóa tính toán địa kỹ thuật mà không cần đầu tư các phần mềm thương mại đắt tiền.
  • Về mặt ứng dụng thực tiễn: Cung cấp quy trình kiểm toán chi tiết từng bước cho kỹ sư tư vấn thiết kế mố cầu móng nông, đặc biệt là việc lựa chọn phân kỳ thi công đắp đất sau mố ($H_đ$) tương thích với tải trọng thiết bị lu lèn nhằm ngăn ngừa tai nạn trượt lở mái dốc trong giai đoạn xây lắp.

Limitations và Future Research

  1. Giả thiết bài toán phẳng biến dạng phẳng (2D Plain Strain): Mô hình hiện tại phân tích trên mặt cắt đại diện 2D với bề rộng đơn vị $B$, chưa xét đến hiệu ứng không gian 3D của tường cánh mố cầu và sự phân tán tải trọng theo phương dọc cầu.
  2. Mô hình vật liệu đàn - dẻo lý tưởng: Đất nền và đất đắp được mô hình hóa theo tiêu chuẩn Mohr-Coulomb đẳng hướng, chưa phản ánh ứng xử phụ thuộc biến dạng (strain-softening) hoặc hiện tượng tăng áp lực nước lỗ rỗng cục bộ khi chịu tải trọng động đất chu kỳ.
  3. Độ phân mảnh của hệ thỏi trượt: Việc chia 6 đến 7 phân tố giúp giải ma trận nhanh chóng nhưng có thể làm mượt cục bộ các điểm kỳ dị ứng suất tại góc đáy bệ mố.

Chương trình nghiên cứu tương lai cần tập trung:

  • Mở rộng thuật toán GLEM sang không gian 3 chiều (3D Limit Equilibrium Method) tích hợp ma trận độ cứng không gian của cọc nếu mố chuyển sang phương án móng cọc ma sát/chống.
  • Ứng dụng trí tuệ nhân tạo (AI) và mạng nơ-ron nhân tạo (Physics-Informed Neural Networks - PINNs) để tìm kiếm bề mặt trượt 3D bất kỳ trong môi trường địa chất phân lớp phức tạp.
  • Đo đạc thực địa bằng thiết bị quan trắc áp lực đất và độ nghiêng (Inclinometer) trên các công trình cầu thật để hiệu chuẩn các hệ số tương tác $k_i, m_i$.

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

  • Ảnh hưởng học thuật: Thiết lập khung tham chiếu giải tích ma trận cho các đề tài nghiên cứu chuyên sâu về ổn định công trình giao thông tại Việt Nam; cung cấp bài toán mẫu có độ tin cậy cao để kiểm định các phần mềm địa kỹ thuật số trị (như GeoStudio Slope/W, PLAXIS 2D, FLAC).
  • Chuyển đổi ngành xây dựng cầu đường: Giúp các tổng công ty tư vấn thiết kế giao thông (TEDI, TEDISouth,...) tối ưu hóa kích thước hình học mố cầu móng nông, giảm khối lượng bê tông từ $10 - 15%$ trong khi vẫn đảm bảo tuyệt đối hệ số an toàn trượt sâu theo tiêu chuẩn TCVN 11823.
  • Chính sách và Tiêu chuẩn hóa: Cung cấp cơ sở khoa học để Bộ Giao thông Vận tải bổ sung các hướng dẫn kỹ thuật chi tiết về kiểm soát ổn định mố cầu trên móng nông trong điều kiện địa hình đồi núi dốc và động đất tại Việt Nam.

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

  • Nghiên cứu sinh và Giảng viên ngành Cầu hầm / Địa kỹ thuật: Tiếp cận thuật toán ma trận hóa giải tích chi tiết, hiểu rõ bản chất cơ học của tương tác liên thỏi mà các phần mềm thương mại thường ẩn giấu.
  • Kỹ sư Chủ trì Thiết kế Cầu Đường: Sở hữu công cụ tính toán nhanh, chính xác thông qua bảng tính Solver để kiểm toán nhanh phương án mố móng nông trước khi lập mô hình số phức tạp.
  • Cơ quan Quản lý Dự án và Tư vấn Giám sát: Nắm vững các kịch bản nguy hiểm nhất trong giai đoạn thi công (đắp đất lưng mố kết hợp xe lu hoạt động) để thiết lập quy trình kiểm soát an toàn công trường nghiêm ngặt.

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 ma trận hóa hoàn toàn hệ phương trình vi phân cân bằng giới hạn tổng quát GLEM thành ma trận tuyến tính cục bộ $[A]{24 \times 24} {K} = {B}$, tích hợp trực tiếp 5 thành phần tĩnh tải kết cấu mố cầu ($Q{M1}$ đến $Q_{M5}$) và áp lực đất đắp ($Q_{\text{soil}}$) vào điều kiện trượt Mohr-Coulomb. Công trình đã mở rộng trực tiếp Lý thuyết Cân bằng Giới hạn Tổng quát của Fredlund & Krahn (1977) và Spencer (1967) sang lĩnh vực tương tác kết cấu mố cầu - nền đất có tải trọng lệch tâm lớn.

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

So với phương pháp Spencer (1967) chỉ giải phương trình phi tuyến một ẩn góc nghiêng lực liên thỏi $\theta$ và phương pháp Morgenstern-Price (1965) cần hàm giả định $f(x)$, phương pháp của luận án thiết lập trực tiếp các thông số định hướng lực tiếp tuyến đáy $k_i$ và liên thỏi $m_i$. Đồng thời, việc ứng dụng thuật toán GRG2 Solver cho phép tối ưu hóa đồng thời 12 biến hình học tọa độ nút mặt trượt $(X_i, Y_i)$ để tìm ra $F_{s,\min}$, thay vì phải quét lưới tâm trượt tròn thủ công như phương pháp Bishop cổ điển.

3. Phát hiện bất ngờ nhất từ số liệu tính toán là gì?

Phát hiện bất ngờ nhất là sự biến đổi phi đối xứng của hệ số an toàn đáy $FS$ khi thay đổi điều kiện liên thỏi $FSI$. Khi giải phóng hoàn toàn lực cắt liên thỏi ($FSI = \infty$), hệ số an toàn trượt qua chân dốc chỉ đạt $F_{s,\min} = 1{,}380$. Nhưng khi đưa vào lực ma sát liên thỏi thực tế ($FSI = 1{,}700$), hệ số an toàn thực tế đạt $FS = 1{,}510$. Điều này chứng minh rằng việc bỏ qua lực cắt giữa các thỏi đất trong các tính toán truyền thống đã đánh giá quá thiên về an toàn (dẫn đến lãng phí vật liệu kết cấu mố) tới $9{,}42%$.

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

Quy trình tái lập được cung cấp hoàn chỉnh và minh bạch 100%: Toàn bộ ma trận số hệ số $[A]$ kích thước $24 \times 24$, véc-tơ tải trọng ${B}$, bảng thông số hình học từng phân tố ($X_i, Y_i, XP_i, YP_i, \alpha_i, R_i, S_i, A_i, W_i$), cùng giá trị tải trọng chi tiết của từng bộ phận mố cầu ($Q_{M1} = 2.829{,}75\text{ kN}, Q_{M2} = 2.263{,}80\text{ kN},...$) được trình bày tường minh, cho phép bất kỳ nhà nghiên cứu nào cũng có thể tái lập chính xác kết quả trong Microsoft Excel hoặc MATLAB/Python.

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

Chương trình 10 năm định hướng:

  • Giai đoạn 2024–2027: Phát triển phần mềm tính toán chuyên dụng GLEM-Abutment 3D tích hợp tiêu chuẩn TCVN 11823:2017.
  • Giai đoạn 2028–2030: Xây dựng mô hình tương tác đất - mố - nhịp cầu dưới tác động của sóng địa chấn phi tuyến thời gian thực (Nonlinear Time-History Analysis).
  • Giai đoạn 2031–2034: Tích hợp công nghệ cảm biến sợi quang (FGB) và mô hình bản sao số (Digital Twin) để dự báo suy giảm hệ số an toàn trượt sâu mố cầu thời gian thực dưới tác động của biến đổi khí hậu và xói lở chân dốc.

Kết luận

  1. Nghiên cứu đã xây dựng thành công mô hình toán học ma trận hóa $24 \times 24$ theo Phương pháp Cân bằng Giới hạn Tổng quát (GLEM), giải quyết triệt để bài toán ổn định trượt sâu mố cầu trên móng nông chịu tổ hợp tải trọng phức tạp.
  2. Thiết lập quy trình tối ưu hóa phi tuyến tự động xác định mặt trượt nguy hiểm nhất và hệ số an toàn cực tiểu ($F_{s,\min}$), loại bỏ các giả định chủ quan về hình học cung trượt tròn truyền thống.
  3. Định lượng hóa chính xác ảnh hưởng của các thông số thi công cốt lõi: chiều cao đắp đất sau mố $H_đ$ ($2{,}05\text{ m}; 4{,}09\text{ m}; 8{,}19\text{ m}$), tải trọng thiết bị thi công (xe lu, xe 2-3 trục), mực nước ngầm và lực quán tính động đất theo các Trạng thái Giới hạn Cường độ.
  4. Chứng minh vai trò quan trọng của lực truyền tương tác giữa các phân tố ($FSI$), chỉ ra mức độ chênh lệch hệ số an toàn từ $9{,}4%$ đến $12{,}6%$ giữa mô hình có xét và không xét ma sát liên thỏi.
  5. Mở ra 3 hướng nghiên cứu chuyên sâu mới: (i) Phân tích ổn định mố cầu trong không gian 3D; (ii) Tối ưu hóa tương tác mố - móng cọc trên sườn dốc; (iii) Ứng dụng thuật toán học máy trong dự báo độ ổn định công trình cầu vượt địa hình đồi núi dốc.
  6. Cung cấp bộ công cụ tính toán số học hoàn chỉnh, có giá trị ứng dụng thực tiễn cao cho công tác thiết kế, thẩm tra và giảng dạy kỹ thuật địa kỹ thuật - công trình giao thông tại Việt Nam và khu vực.