Tổng quan về luận án

Sự bùng nổ của các đô thị hiện đại đặt ra yêu cầu cấp thiết về phát triển mạng lưới giao thông ngầm nhằm giải quyết triệt để vấn đề ùn tắc và ô nhiễm môi trường. Tại Thành phố Hồ Chí Minh, quy hoạch mạng lưới đường sắt đô thị với 6 tuyến metro, 3 tuyến xe điện mặt đất và 2 tuyến đường sắt nhẹ trên cao—tiêu biểu là Tuyến số 1 (Bến Thành – Suối Tiên) và Tuyến số 2 (Bến Thành – Tham Lương)—đang đặt ra những thách thức kỹ thuật rất lớn khi thi công trong điều kiện địa chất phức tạp và mật độ xây dựng dày đặc. Công nghệ khiên đào cân bằng áp lực đất (Earth Pressure Balance - EPB TBM) được xem là giải pháp tối ưu cho môi trường đất yếu đô thị. Tuy nhiên, việc vận hành khiên đào tiềm ẩn nguy cơ mất ổn định gương hầm ở hai thái cực: áp lực chống đỡ không đủ dẫn đến sụp đổ chủ động (active failure), hoặc áp lực kích đẩy quá mức gây phá hoại bị động (passive failure) làm trồi bề mặt đất và phá hủy kết cấu công trình lân cận.

Khoảng trống nghiên cứu then chốt (research gap) mà luận án tiến sĩ kỹ thuật của NCS. Nguyễn Anh Tuấn (chuyên ngành Địa kỹ thuật xây dựng, Trường Đại học Bách khoa – ĐHQG TP.HCM, năm 2018) tập trung giải quyết là: Trong khi phần lớn các công trình kinh điển (Broms & Bennermark, 1967; Davis et al., 1980; Horn, 1961; Kovári & Anagnostou, 1996; Chambon & Corte, 1994) tập trung vào bài toán áp lực chủ động chống sụp đổ gương đào, các nghiên cứu chuyên sâu về áp lực phá hoại bị động, cơ chế trượt đùn và quy luật biến dạng khối đất trước gương hầm trong địa tầng hai lớp (cát chặt vừa phủ trên sét cứng) vẫn còn rất khan hiếm.

Luận án xác lập ba câu hỏi nghiên cứu trọng tâm:

  1. $RQ_1$: Cơ chế phá hoại cục bộ và hình học mặt trượt bị động trước gương hầm phát triển như thế nào trong quá trình khiên đào tiến lên phía trước?
  2. $RQ_2$: Quy luật biến thiên của áp lực phá hoại bị động chuẩn hóa $N_{\gamma m}$ theo tỷ lệ chuyển vị $S_x/D$ và tỷ lệ độ sâu tầng phủ $C/D$ diễn ra ra sao?
  3. $RQ_3$: Làm thế nào để thiết lập một mô hình toán học giải tích xác định chính xác áp lực phá hoại bị động phục vụ kiểm soát vận hành máy EPB-TBM trong điều kiện địa chất tương tự TP.HCM?

Các giả thuyết nghiên cứu được kiểm định bao gồm:

  • $H_1$: Mặt trượt bị động trong môi trường đất cát không phát triển theo dạng nêm phẳng cổ điển mà mở rộng thành hình phễu không gian ba chiều với góc nghiêng bề mặt phá hoại so với phương ngang xấp xỉ $(45^\circ - \varphi'/2)$.
  • $H_2$: Áp lực phá hoại bị động chuẩn hóa tăng phi tuyến theo hành trình kích đẩy và đạt trạng thái ổn định tới hạn (steady-state plateau) tại ngưỡng chuyển vị xác định, đồng thời gia tăng đáng kể khi chiều sâu tầng phủ $C/D$ tăng.

Khung lý thuyết của nghiên cứu được tích hợp từ Lý thuyết phân tích giới hạn động học (Kinematical Limit Analysis - Upper & Lower Bound Theorems) của Leca & Dormieux (1990) và Soubra (2002), kết hợp với Định lý cân bằng giới hạn (Horn, 1961) và Nguyên lý mô hình hóa tương đồng vật lý ly tâm (Centrifuge Scaling Principles). Luận án triển khai hai mô hình thí nghiệm ly tâm tỷ lệ thu nhỏ 1/100 tại Trường Đại học Khoa học và Công nghệ Hong Kong (HKUST) với hai tỷ số tầng phủ $C/D = 1,5$ và $C/D = 3,3$, mô phỏng tốc độ tiến gương hầm $0,2\text{ mm/s}$ (tương đương $15\text{ m/ngày}$ trong nguyên mẫu thực tế) với hành trình dịch chuyển tối đa $40\text{ mm}$. Nghiên cứu đã định lượng hóa sự gia tăng áp lực bị động từ $1,6%$ đến $39%$ khi tăng chiều sâu $C/D$, xác lập phạm vi ảnh hưởng trực tiếp $1,5D$ và vùng ảnh hưởng tổng thể $3D$, đóng góp cơ sở dữ liệu thực nghiệm chuẩn mực cho kỹ thuật ngầm đô thị.

Literature Review và Positioning

Lịch sử nghiên cứu ổn định gương hầm được định hình bởi nhiều trường phái học thuật. Hướng tiếp cận bán thực nghiệm bắt đầu từ Broms & Bennermark (1967) với việc đề xuất hệ số ổn định $N = (\sigma_{ob} - \sigma_T)/S_u$, trong đó xác định ngưỡng mất ổn định gương hầm trong đất sét khi $N > 6$. Atkinson & Mair (1981) và Mair (1981) mở rộng khung phân tích này sang môi trường đất cát bằng cách thiết lập các hệ số ổn định $T_c, T_s, T_\gamma$ tương đồng với hệ số sức chịu tải móng ($N_c, N_q, N_\gamma$), đồng thời đưa ra khái niệm Hệ số tải trọng (Load Factor): $$LF = \frac{\sigma_{v0} - \sigma_i}{\sigma_{v0} - \sigma_{ic}}$$ nhằm đánh giá mức độ an toàn của áp lực chống đỡ gương đào $\sigma_i$ so với giá trị tới hạn $\sigma_{ic}$.

Về mặt lý thuyết cơ học, Horn (1961) đặt nền móng cho phương pháp cân bằng giới hạn với mô hình gồm một khối nêm lăng trụ phía trước gương và khối lăng trụ thẳng đứng phía trên kéo dài lên mặt đất. Mô hình này sau đó được Kovári & Anagnostou (1996) hoàn thiện để tính toán áp lực chống đỡ tối thiểu cho khiên cân bằng áp lực vữa và EPB. Ở nhánh phân tích giới hạn (Limit Analysis), Davis và cộng sự (1980) đề xuất các cơ chế phá hoại cận trên và trường ứng suất cận dưới cho hầm tròn biến dạng phẳng. Bước đột phá về mô hình hóa ba chiều được thiết lập bởi Leca & Dormieux (1990) với cơ chế phá hoại gồm các khối hình nón cụt trượt tương đối, và tiếp tục được Soubra (2002) tối ưu hóa thành cơ chế 5 khối động học, mang lại các giá trị cận trên chặt chẽ hơn cho cả bài toán sụp đổ chủ động ($N_{\gamma c}, N_{sc}$) và bài toán đẩy trồi bị động ($N_{\gamma b}, N_{sb}$). Gần đây hơn, Mollon và cộng sự (2009) đã áp dụng kỹ thuật tối ưu hóa không gian đa khối 3D để mô tả biên phá hoại liên tục trong đất dính-ma sát.

Trong y văn tồn tại những tranh luận sâu sắc giữa các giả định lý thuyết và hiện tượng quan sát thực tế:

  1. Hình học mặt trượt: Các mô hình giải tích cổ điển (Horn, 1961; Kovári & Anagnostou, 1996) giả định mặt trượt phẳng, trong khi các quan sát thực nghiệm và mô phỏng số (Gregor Idinger et al., 2011; Mollon et al., 2009) chỉ ra rằng mặt phá hoại thực tế có dạng cong phi tuyến ba chiều.
  2. Trọng tâm nghiên cứu: Đại đa số các công trình chỉ khảo sát trạng thái chủ động nhằm ngăn chặn sụt lún mặt đất, hoàn toàn bỏ ngỏ diễn biến ứng suất và cơ chế biến dạng khi áp lực đẩy của kích khiên đào vượt ngưỡng cân bằng, gây phá hoại bị động trong nền địa tầng phân lớp.

Về mặt thực nghiệm quốc tế, luận án định vị và so sánh kết quả với các nghiên cứu tiêu biểu:

  • Nghiên cứu của Chambon & Corte (1994) trên máy ly tâm với các tỷ lệ $C/D$ từ $0,5$ đến $4,0$ trong cát Toyoura đã chứng minh mối quan hệ tuyến tính giữa đường kính hầm và áp lực phá hoại chủ động, nhưng chưa định lượng cơ chế động học khi gương hầm chủ động ép trồi vào khối đất.
  • Nghiên cứu của Kuwano (2004) tại gia tốc $80g$ cho thấy vùng phá hoại mở rộng từ đỉnh hầm lên mặt đất trong khoảng $3D - 3,5D$ khi giảm áp lực chống đỡ.
  • Nghiên cứu của Gregor Idinger và cộng sự (2011) tại $50g$ với chuyển vị piston $d_p = 5\text{ mm}$ ($5%D$) chỉ ra áp lực chống đỡ thay đổi tới $50%$ giữa các cấp $C/D = 0,5 - 1,5$, đồng thời ghi nhận xu hướng trượt xuống của khối nêm thay vì trượt ngang.
  • Thí nghiệm của Wong và cộng sự (2012) tại HKUST dưới gia tốc $100g$ cho $C/D = 2,2$ và $4,3$ trong cát bão hòa xác nhận áp lực bị động chuẩn hóa $N_{\gamma m}$ tiệm cận nghiệm cận trên của Soubra (2002) tại $S_x/D = 0,8$.

Luận án của Nguyễn Anh Tuấn đã định vị chính xác khoảng trống này: khảo sát trực tiếp quá trình ép trồi bị động liên tục trong nền hai lớp (cát chặt vừa nằm trên sét cứng) dưới gia tốc ly tâm cao, kết hợp kỹ thuật đo đạc quang học trường hạt (PIV) nhằm lượng hóa toàn diện cả trường ứng suất và trường chuyển vị.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án tạo ra bước tiến quan trọng trong việc mở rộng và kiểm chứng lý thuyết phân tích giới hạn của Leca & Dormieux (1990) và Soubra (2002) đối với trạng thái phá hoại bị động. Nghiên cứu chứng minh rằng nghiệm giải tích tĩnh cổ điển không phản ánh đầy đủ quá trình phát triển lũy tiến của khối trượt khi khiên đào chuyển vị tịnh tiến.

Một đóng góp đột phá về mặt cơ học đất là việc xác lập góc nghiêng của mặt phá hoại bị động trước gương hầm tròn. Trong lý thuyết Rankine cổ điển cho tường chắn biến dạng phẳng chịu áp lực bị động, góc tạo bởi mặt trượt và phương ngang là $(45^\circ - \varphi'/2)$ đối với mặt phẳng đứng đẩy đất. Tuy nhiên, đối với bài toán gương hầm hình tròn ba chiều bị đẩy ép vào lòng đất, luận án đã kiểm chứng bằng dữ liệu ảnh động học thực nghiệm rằng góc nghiêng thực tế của mặt trượt bị động vẫn tuân theo quy luật xấp xỉ $(45^\circ - \varphi'/2)$, trong đó $\varphi'$ là góc ma sát trong hữu hiệu của lớp đất chứa hầm.

Khung lý thuyết được cụ thể hóa thông qua ba mệnh đề định lượng:

  • Mệnh đề 1: Cơ chế phá hoại cục bộ trước gương hầm tương tự như cơ chế cắt cục bộ (local shear failure), sau đó phát triển liên tục và khi lý tưởng hóa thành đường bao trượt trơn sẽ tạo thành một hình phễu không gian 3D mở rộng hướng lên mặt đất.
  • Mệnh đề 2: Áp lực phá hoại bị động chuẩn hóa $N_{\gamma m}$ không phải là hằng số tĩnh mà là một hàm biến thiên theo tỷ số chuyển vị chuẩn hóa $S_x/D$, trải qua hai giai đoạn: gia tăng nhanh và đạt trạng thái ổn định tiệm cận tại ngưỡng tới hạn $S_x/D = 0,7$.
  • Mệnh đề 3: Chiều sâu chôn hầm $C/D$ tác động phi tuyến lên khả năng kháng bị động của khối đất; việc tăng $C/D$ từ $1,5$ lên $3,3$ làm gia tăng giá trị $N_{\gamma m}$ từ $1,6%$ ở giai đoạn đầu lên đến $39%$ khi đạt trạng thái ổn định.

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

Khung phân tích của luận án tích hợp ba trụ cột phương pháp luận:

  1. Lý thuyết biến dạng động học môi trường liên tục: Sử dụng kỹ thuật PIV để ghi nhận trường vector dịch chuyển vi mô của các hạt đất cát theo thời gian thực.
  2. Lý thuyết tương đồng vật lý địa kỹ thuật (Geotechnical Centrifuge Similitude): Tái tạo chính xác trường ứng suất bản thân nguyên mẫu tại mọi độ sâu thông qua trường trọng lực nhân tạo $N \cdot g$.
  3. Phương pháp phân tích hồi quy thực nghiệm đa biến: Xây dựng phương trình giải tích xác định hệ số áp lực bị động chuẩn hóa $N_{\gamma m}$ dựa trên tổ hợp tham số không thứ nguyên: $$N_{\gamma m} = f\left(\frac{C}{D}, \frac{S_x}{D}, \varphi'\right)$$ với các điều kiện biên xác định: áp dụng cho nền đất cát trạng thái chặt vừa phủ trên tầng sét cứng, tỷ lệ độ sâu $1,5 \le C/D \le 3,3$, và bỏ qua ảnh hưởng mất mát thể tích đào phía đuôi khiên.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), sử dụng phương pháp mô hình hóa vật lý ly tâm (Centrifuge Physical Modeling) để tái hiện hiện tượng cơ học đất nguyên mẫu với độ tin cậy cao nhất.

Thiết kế nghiên cứu đa tầng bậc (Multi-level Experimental Design) bao gồm:

  • Mức độ vi mô/trung mô (Meso-level): Quan sát trường vector chuyển vị của các cụm hạt đất cát bằng camera kỹ thuật số phân giải cao gắn phía ngoài thành hộp kính.
  • Mức độ tiếp xúc cục bộ (Interface-level): Ghi nhận trực tiếp áp lực pháp tuyến tác dụng lên bề mặt đĩa cắt của gương hầm thông qua hệ thống cảm biến lực (load cells).
  • Mức độ vĩ mô toàn thể (Macro-level): Đo đạc biến dạng lún/trồi của toàn bộ bề mặt địa hình bằng mảng biến trở đo dịch chuyển tuyến tính (LVDT).

Kế hoạch thực nghiệm gồm hai bài thí nghiệm ly tâm trọng điểm tại máy ly tâm địa kỹ thuật dạng cánh tay đòn 400g-ton của HKUST:

  • Thí nghiệm T1: Tỷ lệ chiều sâu tầng phủ $C/D = 1,5$.
  • Thí nghiệm T2: Tỷ lệ chiều sâu tầng phủ $C/D = 3,3$.

Quy trình nghiên cứu rigorous

Quy trình chuẩn bị và vận hành mô hình thí nghiệm được kiểm soát nghiêm ngặt theo các tiêu chuẩn thực nghiệm địa kỹ thuật quốc tế:

  • Chế bị mẫu đất: Mô hình được thực hiện trong thùng thép chữ nhật cứng có kích thước hình học chính xác, một mặt gắn kính acrylic trong suốt cường lực nhằm phục vụ quan trắc quang học. Lớp sét cứng phía dưới được đầm nén và cố kết đồng nhất; lớp cát hạt thô/vừa phía trên được tạo lập bằng phương pháp rải cát mưa (sand pluviation/raining method) nhằm kiểm soát chặt chẽ dung trọng khô $\gamma_d$ và đạt trạng thái chặt vừa đồng đều ($e \approx 0,65 - 0,75$).
  • Mô hình hóa đường hầm và hệ gia tải: Hầm mô hình tỷ lệ thu nhỏ $1/100$ được chế tạo từ hợp kim nhôm cứng, trang bị đầu kích dẫn động bằng động cơ servo điều khiển số, cho phép gương hầm tịnh tiến ra phía trước với vận tốc không đổi $v = 0,2\text{ mm/s}$ (tương đương tốc độ đào thực tế $15\text{ m/ngày}$ theo nguyên tắc tỷ lệ thời gian động học). Quá trình ép dừng lại khi chuyển vị đạt $40\text{ mm}$ ($S_x/D = 0,8$).
  • Hệ thống thiết bị đo đạc:
    • Mảng cảm biến LVDT (ký hiệu LA, LB, LC, LD...) bố trí dọc tim tuyến và mặt cắt ngang trên bề mặt để ghi nhận biến dạng trồi đất liên tục.
    • Kỹ thuật PIV (Particle Image Velocimetry): Camera kỹ thuật số tốc độ cao chụp ảnh liên tục bề mặt mẫu đất qua cửa sổ trong suốt. Thuật toán tương quan ảnh số (Digital Image Correlation) tính toán trường dịch chuyển và biến dạng trượt của khối đất trước gương theo từng bước gia tải $S_x/D = 0,1; 0,3; 0,5; 0,7; 0,8$.
    • Cảm biến lực (Load cells) gắn trực tiếp tại đầu gương hầm, được hiệu chuẩn tuyến tính trước thí nghiệm để ghi nhận chính xác áp lực bị động $N_{\gamma m}$.
  • Kiểm soát sai số và độ tin cậy: Giảm thiểu ma sát thành thùng bằng màng bôi trơn chuyên dụng; hiệu chỉnh biến dạng của thùng chứa dưới gia tốc hướng tâm; kiểm chứng chéo (triangulation) giữa số liệu chuyển vị đo từ LVDT trên mặt đất với tích phân trường chuyển vị từ PIV.

Data và phân tích

Dữ liệu thực nghiệm được chuẩn hóa dựa trên các định luật tỉ lệ ly tâm (Centrifuge Scaling Laws): với hệ số gia tốc $N = 100$, tỉ lệ ứng suất là $1:1$, tỉ lệ chuyển vị là $1:100$, và tỉ lệ lực là $1:100^2$.

Các thông số cơ bản được xử lý thống kê:

  • Tỷ số chuyển vị chuẩn hóa: $S_x/D$ (với $S_x$ là hành trình tiến gương, $D$ là đường kính hầm).
  • Hệ số áp lực bị động chuẩn hóa: $$N_{\gamma m} = \frac{\sigma_{pm}}{\gamma D}$$ trong đó $\sigma_{pm}$ là áp lực bị động trung bình đo tại gương, $\gamma$ là dung trọng đất.

Phương pháp phân tích hồi quy phi tuyến đa biến (Non-linear regression analysis) trên phần mềm chuyên dụng được áp dụng để xác định các hệ số tương quan $K_0, K_1, \dots, K_7$ và $A_1, A_2, \dots, A_4$, thiết lập phương trình dự báo đạt hệ số xác định $R^2 > 0,95$.

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

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

Từ kết quả phân tích ảnh PIV, số liệu cảm biến lực và đo đạc LVDT, luận án đã rút ra 5 phát hiện cốt lõi:

  1. Hình thái mặt trượt không gian dạng phễu: Dữ liệu thực nghiệm xác nhận nguyên văn: "Khi các cơ chế phá hoại quan sát được lý tưởng hóa bằng đường liền, cơ chế phá hoại đất trước gương hầm có dạng phễu." Đất ở sát mặt gương bị đẩy dịch chuyển tịnh tiến ra phía trước, trong khi đất ở xa hơn bị ép trồi xiên lên trên mặt đất tạo thành một vùng phá hoại cắt cục bộ mở rộng.
  2. Quy luật góc nghiêng mặt trượt bị động: Mặt trượt bị động tạo với mặt phẳng ngang một góc xác định: $$\theta \approx 45^\circ - \frac{\varphi'}{2}$$ với $\varphi'$ là góc ma sát trong hữu hiệu của lớp đất cát. Đây là minh chứng thực nghiệm then chốt giúp chuẩn hóa mô hình hình học cho các phân tích cân bằng giới hạn 3D.
  3. Quy mô không gian của vùng ảnh hưởng: Bề rộng vùng phá hoại trực tiếp trên mặt đất xuất hiện cách gương hầm một khoảng xấp xỉ $1,5D$. Phạm vi ảnh hưởng tổng thể do quá trình kích ép khiên đào tác động đến các công trình lân cận lan rộng tới $3,0D$.
  4. Động học phát triển áp lực bị động $N_{\gamma m}$:
    • Đối với hầm đặt nông ($C/D = 1,5$): $N_{\gamma m}$ tăng đều khi $S_x/D < 0,3$, sau đó tốc độ tăng giảm dần và đạt tới trạng thái ổn định (steady state) khi $S_x/D = 0,7$.
    • Đối với hầm đặt sâu ($C/D = 3,3$): $N_{\gamma m}$ tăng rất nhanh và mạnh trong suốt giai đoạn $S_x/D < 0,6$, và cũng đạt trạng thái ổn định tại ngưỡng $S_x/D = 0,7$.
  5. Mức độ gia tăng áp lực theo độ sâu: Khi chiều sâu đặt hầm tăng từ $C/D = 1,5$ lên $3,3$, giá trị áp lực phá hoại bị động tại gương hầm gia tăng tương ứng từ $1,6%$ (ở giai đoạn chuyển vị nhỏ ban đầu) lên tới $39%$ (ở giai đoạn phá hoại toàn diện).

Phương trình toán học thực nghiệm tính toán áp lực phá hoại bị động cho đất cát chặt vừa được luận án xác lập có dạng: $$N_{\gamma m} = f\left(\frac{C}{D}, \frac{S_x}{D}, \varphi'\right) = \left[ -0{,}012\left(\frac{C}{D}\right)^3 + 0{,}105\left(\frac{C}{D}\right)^2 - 0{,}28\left(\frac{C}{D}\right) + 1{,}15 \right] \cdot \dots$$

Implications đa chiều

  • Về mặt học thuật: Luận án cung cấp bộ dữ liệu thực nghiệm ly tâm chuẩn mực đầu tiên về phá hoại bị động trong nền phân lớp, hoàn thiện phổ lý thuyết ổn định gương hầm từ trạng thái sụp đổ chủ động (active collapse) đến trạng thái trồi đùn bị động (passive blow-out).
  • Về phương pháp luận: Khẳng định tính ưu việt của việc kết hợp mô hình ly tâm gia tốc cao ($100g$) với công nghệ phân tích ảnh quang học PIV trong việc bóc tách cơ chế biến dạng nội tại của khối đất mà các thí nghiệm 1g truyền thống không thể mô phỏng chính xác.
  • Về thực tiễn kỹ thuật thi công: Cung cấp công cụ tính toán giới hạn áp lực chống đỡ trên (Upper support pressure limit) cho kỹ sư vận hành máy TBM. Việc khống chế áp lực buồng đào không vượt quá ngưỡng $N_{\gamma m}$ giúp triệt tiêu nguy cơ nứt vỡ mặt đường, trồi đất và làm hư hỏng móng cọc công trình lân cận.
  • Về chính sách và tiêu chuẩn: Làm cơ sở định lượng để xây dựng quy trình kiểm soát an toàn thi công hầm metro trong đô thị cho các cơ quan quản lý (Ban Quản lý Đường sắt Đô thị MAUR, Sở Giao thông Vận tải).

Limitations và Future Research

Luận án thừa nhận một số giới hạn nghiên cứu xuất phát từ điều kiện thực nghiệm:

  1. Mô hình địa tầng lý tưởng hóa: Nghiên cứu chỉ tập trung vào mô hình nền hai lớp (cát chặt vừa nằm trên sét cứng), chưa khảo sát các địa tầng phức tạp hơn như cát lẫn sét, thấu kính bùn hữu cơ hoặc đất sét dẻo chảy.
  2. Điều kiện thủy lực tĩnh: Thí nghiệm tiến hành trong điều kiện đất cát khô trên nền sét, chưa xét đến tác động của áp lực nước lỗ rỗng động và gradient thủy lực xuất hiện khi đào hầm dưới mực nước ngầm cao có dòng thấm.
  3. Giới hạn cơ chế hình học khiên: Thí nghiệm mô phỏng chuyển vị tịnh tiến liên tục của gương hầm tròn phẳng, chưa xét đến ma sát dọc thân khiên (skin friction), độ côn của vỏ khiên (shield taper), mất mát thể tích đào phía đuôi khiên (tail void gap) và áp lực ép vữa chèn lấp (grouting pressure).
  4. Chế độ kích tiến cố định: Vận tốc đẩy được duy trì không đổi ($0,2\text{ mm/s}$), chưa xét đến chu kỳ ngắt quãng thực tế khi dừng đào để lắp ghép các phân đoạn vỏ hầm (segmental lining).

Chương trình nghiên cứu tương lai (Future Research Agenda) cần tập trung vào 4 hướng:

  • Hướng 1: Mở rộng thí nghiệm ly tâm cho mô hình ba chiều có xét đến dòng thấm thủy lực và áp lực nước ngầm dâng cao.
  • Hướng 2: Tích hợp mô phỏng số 3D kết hợp giữa Phương pháp Phần tử Hữu hạn và Phần tử Rời rạc (FEM - DEM) để phân tích chi tiết tương tác tiếp xúc giữa các răng cắt của đầu khiên và hạt đất.
  • Hướng 3: Khảo sát ảnh hưởng của tải trọng động từ công trình cao tầng lân cận và cọc móng hiện hữu nằm trong phạm vi vùng ảnh hưởng $1,5D - 3,0D$.
  • Hướng 4: Nghiên cứu ảnh hưởng của quá trình bơm vữa áp lực cao sau đuôi khiên kết hợp với áp lực gương đào đến trường biến dạng tổng thể của khối đất.

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

Nghiên cứu tạo ra tác động sâu rộng trên nhiều phương diện:

  • Tác động học thuật: Luận án là tài liệu tham khảo nền tảng cho các nghiên cứu tiếp theo về công trình ngầm đô thị tại Việt Nam, mở ra tiềm năng trích dẫn cao trong các tạp chí chuyên ngành quốc tế uy tín (như Tunnelling and Underground Space Technology, Géotechnique, Acta Geotechnica).
  • Chuyển giao công nghệ và công nghiệp: Cung cấp giải pháp tính toán trực tiếp cho các tổng thầu thi công TBM (như Shimizu, Maeda, Fecon, v.v.) trong việc lập biện pháp thi công và cài đặt thông số áp lực vận hành buồng cân bằng áp lực đất EPB cho Tuyến Metro số 1 và Tuyến số 2 TP.HCM.
  • Hiệu quả kinh tế - xã hội: Việc kiểm soát chính xác áp lực bị động giúp ngăn ngừa triệt đẻ các sự cố trồi sụt mặt đường, bảo vệ an toàn tuyệt đối cho các công trình di sản kiến trúc (Nhà hát Lớn, Trụ sở UBND Thành phố) và mạng lưới hạ tầng kỹ thuật ngầm, tiết kiệm hàng triệu USD chi phí xử lý sự cố và đền bù thiệt hại.

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

  • Nghiên cứu sinh và Giảng viên ngành Địa kỹ thuật: Thụ hưởng bộ dữ liệu thực nghiệm ly tâm chất lượng cao, quy trình ứng dụng công nghệ PIV và công thức giải tích $N_{\gamma m}$ phục vụ công tác giảng dạy và phát triển các đề tài nghiên cứu chuyên sâu.
  • Kỹ sư thiết kế và Tư vấn giám sát: Sử dụng công thức đề xuất để xác định dải áp lực vận hành an toàn $[\sigma_{active}, \sigma_{passive}]$ của buồng đào khiên TBM, thiết lập ngưỡng cảnh báo sớm trong hệ thống quan trắc công trình.
  • Ban Quản lý dự án đường sắt đô thị (MAUR): Có cơ sở khoa học độc lập, chuẩn xác để thẩm định biện pháp thi công, đánh giá rủi ro địa kỹ thuật và ban hành các chỉ dẫn kỹ thuật thi công hầm qua các khu vực đô thị nhạy cảm.
  • Cộng đồng xã hội: Hưởng lợi từ sự an toàn của hệ thống hạ tầng đô thị, giảm thiểu tối đa tình trạng gián đoạn giao thông và rủi ro sụt lún công trình trong quá trình xây dựng ngầ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ì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng và kiểm chứng thành công Lý thuyết phân tích giới hạn cận trên (Upper-Bound Limit Analysis) của Leca & Dormieux (1990) và Soubra (2002) cho bài toán phá hoại bị động khiên đào trong nền hai lớp. Nghiên cứu xác lập rằng mặt trượt bị động 3D có hình phễu với góc nghiêng bề mặt phá hoại so với phương ngang xấp xỉ $\theta \approx 45^\circ - \varphi'/2$, đồng thời chứng minh áp lực bị động $N_{\gamma m}$ phát triển phi tuyến theo chuyển vị và đạt trạng thái ổn định tiệm cận tại $S_x/D = 0,7$.

2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với các công trình quốc tế?

So với các nghiên cứu của Chambon & Corte (1994) (chỉ tập trung vào phá hoại chủ động trong cát đồng nhất) và Kuwano (2004) (nghiên cứu mất áp lực chống đỡ ở 80g), luận án đã tạo bước đột phá khi: (1) Mô phỏng trực tiếp quá trình kích đẩy bị động liên tục ($v = 0,2\text{ mm/s}$) trong nền phân tầng (cát trên sét); (2) Ứng dụng đồng thời kỹ thuật PIV để ghi nhận trường chuyển vị liên tục và hệ cảm biến lực đa điểm tại mặt gương để đo trực tiếp áp lực $N_{\gamma m}$; (3) Thực hiện ở gia tốc trọng trường cao $100g$, tái hiện hoàn hảo ứng suất nguyên mẫu.

3. Phát hiện bất ngờ nhất (counter-intuitive) có sự hỗ trợ của dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là sự tồn tại của một "ngưỡng chuyển vị ổn định" (steady-state displacement threshold) đồng nhất tại $S_x/D = 0,7$ cho cả hai trường hợp chiều sâu hầm khác nhau ($C/D = 1,5$ và $C/D = 3,3$). Mặc dù tỷ lệ $C/D = 3,3$ làm tăng áp lực phá hoại bị động lên tới $39%$ so với $C/D = 1,5$, khối đất đều đạt đến trạng thái phá hoại chảy dẻo hoàn toàn tại cùng một tỷ lệ chuyển vị chuẩn hóa $S_x/D = 0,7$.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (replication protocol) đầy đủ không?

Có. Luận án trình bày chi tiết toàn bộ thông số kỹ thuật: thiết bị ly tâm 400g-ton tại HKUST, kích thước thùng chứa, quy trình rải cát mưa kiểm soát dung trọng, thông số lớp sét, kích thước mô hình hầm tỷ lệ $1/100$, vận tốc kích đẩy servo $0,2\text{ mm/s}$, sơ đồ lắp đặt mạng lưới LVDT, thông số camera và điểm kiểm soát quang học phục vụ thuật toán PIV.

5. Khung chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Chương trình nghiên cứu dài hạn tập trung vào: (1) Phát triển mô hình giải tích 3D ngẫu nhiên xét đến tính dị hướng của địa tầng; (2) Tích hợp thuật toán trí tuệ nhân tạo (AI/Machine Learning) dự báo áp lực gương hầm thời gian thực dựa trên dữ liệu quan trắc PIV và LVDT; (3) Mở rộng nghiên cứu tương tác động học giữa khiên TBM và hệ thống cọc móng công trình ngầm lân cận trong đô thị thông minh.

Kết luận

  1. Luận án đã làm sáng tỏ toàn diện cơ chế phá hoại bị động của khối đất trước gương hầm trong nền hai lớp (cát chặt vừa trên sét cứng) dưới tác động của khiên đào EPB-TBM thông qua mô hình ly tâm tiên tiến tỷ lệ 1/100 tại HKUST.
  2. Xác định chính xác hình học vùng phá hoại đất có dạng phễu không gian ba chiều, với góc nghiêng bề mặt phá hoại so với phương ngang xấp xỉ $(45^\circ - \varphi'/2)$.
  3. Định lượng hóa phạm vi phá hoại trực tiếp trên mặt đất là $1,5D$ và phạm vi ảnh hưởng tổng thể đến các công trình lân cận là $3,0D$.
  4. Phát hiện quy luật biến thiên của hệ số áp lực bị động chuẩn hóa $N_{\gamma m}$, đạt trạng thái ổn định tới hạn tại $S_x/D = 0,7$, và tăng từ $1,6%$ đến $39%$ khi tăng chiều sâu chôn hầm từ $C/D = 1,5$ lên $3,3$.
  5. Thiết lập thành công công thức giải tích thực nghiệm xác định áp lực phá hoại bị động $N_{\gamma m} = f(C/D, S_x/D, \varphi')$, mang lại công cụ tính toán đắc lực cho kỹ sư công trình ngầm.
  6. Mở ra 3 nhánh nghiên cứu mới: mô hình hóa tương tác đất - khiên có xét đến dòng thấm áp lực cao, mô phỏng vi cơ học hạt FEM-DEM cho quá trình cắt phá đất, và phát triển hệ thống điều khiển thông minh tự động điều chỉnh áp lực buồng đào TBM theo thời gian thực.