Tổng quan về luận án

Sự bùng nổ đô thị hóa và nhu cầu tối ưu hóa không gian ngầm tại các siêu đô thị đòi hỏi sự chuyển dịch từ các kết cấu hầm tròn truyền thống sang hầm có tiết diện đặc biệt. Đường hầm tiết diện chữ nhật cong (sub-rectangular tunnel) thi công bằng công nghệ khiên đào TBM (Tunnel Boring Machine) phi tròn xuất hiện như một giải pháp đột phá, giúp nâng cao hệ số sử dụng không gian mặt cắt ngang, giảm khối lượng đào đắp và kiểm soát lún bề mặt hiệu quả hơn so với hầm tròn tương đương diện tích (như trường hợp tuyến hầm ngầm cao tốc Thượng Hải). Tuy nhiên, khả năng chịu tải trọng động đất của loại kết cấu hình học phức tạp này đặt ra thách thức lớn đối với kỹ thuật công trình ngầm.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định: Trong khi các giải pháp giải tích kinh điển (Wang, 1993; Penzien, 2000) và phương pháp số hiện hành chủ yếu tập trung vào hầm tròn (circular tunnels) hoặc hầm chữ nhật (rectangular tunnels), ứng xử động lực học và nội lực gia tăng trong vỏ hầm chữ nhật cong dưới tác động của sóng địa chấn hầu như chưa được giải quyết một cách tường minh. Các nghiên cứu hiện hữu về hầm chữ nhật cong (Do et al., 2014; Huang et al., 2017; Du et al., 2021) chỉ giới hạn ở điều kiện tải trọng tĩnh (static loading). Hơn nữa, phương pháp phản lực đàn hồi (Hyperstatic Reaction Method - HRM) tuy có ưu thế vượt trội về tốc độ tính toán và tối ưu hóa kết cấu nhưng hoàn toàn thiếu một sơ đồ tải trọng tĩnh tương đương (quasi-static loading scheme) dành cho hình học chữ nhật cong chịu tải trọng động đất.

Nghiên cứu giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết (Hypotheses - H):

  • RQ1: Cơ chế biến dạng và phân bố nội lực gia tăng (mô men uốn $\Delta M$, lực dọc $\Delta N$) của vỏ hầm chữ nhật cong khác biệt như thế nào so với hầm tròn tương đương dưới tác động của sóng cắt địa chấn lan truyền vuông góc?
    • H1: Tiết diện chữ nhật cong với bán kính biến thiên ở vòm ($R_1$), góc vai ($R_2$) và tường biên ($R_3$) sẽ tạo ra sự tập trung ứng suất và biến dạng cục bộ ("sub-ovaling") tại vị trí góc vai, phá vỡ tính đối xứng góc $45^\circ$ vốn thấy ở hầm tròn.
  • RQ2: Tương tác tiếp xúc đất - vỏ hầm (điều kiện trượt hoàn toàn full-slip và không trượt no-slip) cùng các thông số địa kỹ thuật ($E_s, \nu_s$) và địa chấn ($a_H$) chi phối độ lớn nội lực cực trị ra sao?
    • H2: Điều kiện không trượt (no-slip) luôn sinh ra lực dọc gia tăng lớn nhất, trong khi mô đun đàn hồi của đất ($E_s$) và gia tốc cực đại nền ($a_H$) kiểm soát phi tuyến tỷ lệ độ mềm kết cấu (flexibility ratio $F$).
  • RQ3: Làm thế nào để thiết lập một sơ đồ tải trọng tĩnh tương đương mới trong khuôn khổ phương pháp HRM nhằm dự báo chính xác nội lực kết cấu với chi phí tính toán thấp?
    • H3: Một sơ đồ phân bố tải trọng chủ động kết hợp hệ số lò xo đàn hồi biến thiên theo hàm góc phương vị và bán kính cục bộ có thể tái hiện chính xác kết quả mô phỏng số sai phân hữu hạn (FDM).

Khung lý thuyết của luận án dựa trên nguyên lý tương tác động lực học đất - kết cấu (Soil-Structure Interaction - SSI), lý thuyết sóng đàn hồi biến dạng phẳng 2D, và phương pháp phần tử hữu hạn lò xo phi tuyến HRM. Nghiên cứu thực hiện phân tích số chi tiết mô hình 2D FDM (FLAC) đối chứng với lời giải giải tích của Wang (1993), khảo sát dải thông số rộng: gia tốc $a_H$ từ $0.1g$ đến $0.4g$, mô đun đất $E_s$ từ $20\text{ MPa}$ đến $200\text{ MPa}$, chiều dày vỏ $t$ từ $0.3\text{ m}$ đến $0.6\text{ m}$, và độ sâu chôn hầm $H$. Kết quả cung cấp công cụ tính toán giảm hơn 90% thời gian mô phỏng so với mô hình động lực học liên tục, tạo tiền đề tối ưu hóa thiết kế kháng chấn cho hạ tầng giao thông ngầm đô thị.

Literature Review và Positioning

Các nghiên cứu địa chấn công trình ngầm bắt nguồn từ nhận định của Owen và Scholl (1981): "The behavior of an underground structure during a seismic event can be approximated to the one of an elastic beam subjected to deformations imposed by the surrounding ground". Dưới tác động của sóng địa chấn, ba dạng biến dạng chính xuất hiện: nén/dãn dọc trục (axial compression/extension), uốn dọc trục (longitudinal bending), và biến dạng méo mặt cắt ngang (ovalling/racking). Trong đó, Penzien (2000) và Hashash et al. (2001) khẳng định biến dạng méo mặt cắt ngang do sóng cắt ngang (S-wave) gây ra là nguyên nhân nguy hiểm nhất dẫn đến phá hủy kết cấu vỏ hầm.

Tổng quan các trường phái nghiên cứu cho thấy những tranh luận học thuật sâu sắc:

  1. Trường phái Lời giải giải tích (Closed-form solutions): Wang (1993) phát triển lời giải bán giải tích cho hầm tròn dựa trên hệ số độ mềm $F$ và hệ số chịu nén $C$, chỉ ra rằng "relative flexibility as a key parameter for understanding the seismic-induced distortions of underground structures interacting with the surrounding ground". Tuy nhiên, có sự bất đồng lớn giữa Wang (1993) và Penzien (2000) đối với điều kiện tiếp xúc không trượt (no-slip). Hashash et al. (2005) và Kouretzis et al. (2006) thông qua các mô hình số đã chứng minh rằng giải pháp của Penzien đánh giá thấp lực dọc khi xem xét ma sát tiếp xúc, đồng thời khuyến nghị áp dụng phương pháp của Wang cho điều kiện no-slip. Đối với hầm chữ nhật, Wang (1993) thiết lập hệ số méo (racking coefficient $R_r$), chứng minh $R_r$ phụ thuộc chủ yếu vào độ cứng tương đối giữa kết cấu và môi trường hơn là đặc trưng chuyển động đất.
  2. Trường phái Mô hình vật lý thực nghiệm (Physical testing): Các thí nghiệm mô hình ly tâm (centrifuge tests) của Lanzano et al. (2012), Bilotta et al. (2014), Tsinidis et al. (2015, 2016) trên hầm tròn và hộp chữ nhật trong cát khô ghi nhận sự tồn tại của nội lực dư tích lũy (residual internal forces) sau động đất, phụ thuộc chặt chẽ vào gia tốc đỉnh $a_H$. Các thí nghiệm bàn rung 1g (Guobo et al., 2012; Xin et al., 2018) nghiên cứu hiện tượng trượt nứt, khớp nối phân đoạn (segmental joints) và hiệu ứng tương tác công trình lân cận, song gặp hạn chế nghiêm trọng về quy luật đồng dạng vật liệu khi vỏ hầm bước vào giai đoạn nứt dẻo.
  3. Trường phái Phương pháp số và Tĩnh tương đương (Numerical & Quasi-static methods): Sederat et al. (2009) dùng phần mềm ADINA phân tích hầm tròn, kết luận rằng điều kiện tiếp xúc quyết định trực tiếp tới mô men và lực dọc cực trị. Do et al. (2015) mở rộng mô hình FDM FLAC đánh giá độ cứng quay của khớp lắp ghép. Đối với hầm chữ nhật cong, Jianbin Li (2017) và Huang et al. (2017) thực hiện thí nghiệm gia tải quy mô thực (full-scale loading test) để đánh giá trọng lượng bản thân và lắp ghép phân đoạn, trong khi Do et al. (2014) và Du et al. (2021) áp dụng HRM cho điều kiện tải trọng tĩnh.
Phương pháp phân tích Hầm tròn (Circular) Hầm chữ nhật (Rectangular) Hầm chữ nhật cong (Sub-rectangular)
Giải tích (Analytical) Wang (1993); Penzien (2000); Bobet (2003) Wang (1993); Penzien (2000); Huo et al. (2006) Chưa có lời giải giải tích tường minh
Mô hình vật lý (Physical) Lanzano et al. (2012); Tsinidis et al. (2016) Amorosi et al. (2014); Ulgen et al. (2015) Huang et al. (2017) (chỉ xét tải trọng tĩnh)
Mô hình số / HRM (Numerical) Sederat et al. (2009); Do et al. (2015) Hashash et al. (2001); Chen et al. (2020) Do et al. (2014); Luận án Phạm Văn Vĩ (2022)

Luận án định vị chính xác tại điểm giao thoa chưa từng được khám phá: Mở rộng cơ sở lý thuyết kháng chấn từ hầm tròn/chữ nhật sang kết cấu hình học chữ nhật cong, đồng thời nâng cấp phương pháp HRM từ trạng thái tĩnh sang trạng thái động đất tĩnh tương đương.

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

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

Luận án đóng góp vào hệ thống lý thuyết cơ học công trình ngầm thông qua ba trụ cột:

  1. Phát triển lý thuyết biến dạng méo phi đối xứng (Sub-ovalling Kinematics): Nghiên cứu chỉ ra rằng dưới tác dụng của ứng suất cắt trường xa $\tau = \gamma_{\max} G$, cơ chế biến dạng của hầm chữ nhật cong không tuân theo quy luật oval hóa đối xứng trục $45^\circ$ như công thức của Wang (1993) trên hầm tròn. Do sự thay đổi đột ngột bán kính cong giữa vòm đỉnh ($R_1$), vòm vai ($R_2$) và vách đứng ($R_3$), dạng hình học bị biến dạng uốn cục bộ tại góc vai ($R_2 \ll R_1, R_3$). Hiện tượng này tạo ra mô hình phân bố nội lực gia tăng mới với cực trị mô men uốn tập trung cục bộ tại các góc chuyển tiếp.
  2. Xác lập mô hình toán học tương tác tiếp xúc phi tuyến: Nghiên cứu mở rộng hệ thức tính toán ứng suất tiếp xúc của Kouretzis et al. (2006) sang biên hình học nhiều bán kính cong, cung cấp bằng chứng định lượng về sự thay đổi tỷ lệ giữa lực dọc nén và kéo khi xét hệ số ma sát tiếp xúc từ full-slip đến no-slip.
  3. Chuyển dịch mô hình tính toán kháng chấn (Paradigm Shift): Thay vì phải giải bài toán động lực học phi tuyến theo lịch sử thời gian (time-history dynamic analysis) tốn kém tài nguyên tính toán, luận án chứng minh rằng trạng thái ứng suất - biến dạng giới hạn của vỏ hầm chữ nhật cong hoàn toàn có thể được mô phỏng chính xác thông qua sơ đồ tựa tĩnh tương đương được hiệu chỉnh theo trường biến dạng cắt cực đại $\gamma_{\max} = V_{\max}/C_m$.
[Trường biến dạng cắt tự do gamma_max] 
[Ứng suất cắt trường xa: tau = G * gamma_max]
[Mô hình số liên tục FDM (FLAC3D)]      [Khung phân tích mới HRM Tĩnh tương đương]
  - Trạng thái biến dạng phẳng 2D          - Phần tử thanh dầm đàn hồi (EA, EI)
  - Điều kiện Full-slip & No-slip         - Hệ lò xo phản lực nền phi tuyến (kn, ks)
  - Biến dạng méo Sub-ovalling            - Sơ đồ tải trọng chủ động hiệu chỉnh (a, b, beta)
      [Kiểm chứng tương thích & Tối ưu hóa kết cấu vỏ hầm]

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

Khung phân tích tích hợp ba lý thuyết nền tảng: Lý thuyết lan truyền sóng trong môi trường đàn hồi, Lý thuyết tương tác kết cấu - môi trường Winkler mở rộng, và Phương pháp phần tử hữu hạn dầm - lò xo HRM (Hyperstatic Reaction Method).

Trong mô hình HRM cải tiến:

  • Kết cấu vỏ hầm được rời rạc hóa thành các phần tử thanh dầm 2 nút với 3 bậc tự do mỗi nút ($u, v, \theta$), đặc trưng bởi độ cứng chống uốn $EI$ và độ cứng chống nén $EA$.
  • Tương tác với môi trường đất xung quanh được mô hình hóa qua hệ thống lò xo pháp tuyến ($k_n$) và lò xo tiếp tuyến ($k_s$). Áp lực phản lực thụ động của đất $p$ liên hệ phi tuyến với chuyển vị pháp tuyến $\delta$ qua độ cứng ban đầu $\eta_{n,0}$ và áp lực giới hạn $p_{\lim}$:

$$p = \eta_{n,0} \cdot \delta \left(1 - \frac{\eta_{n,0} \cdot \delta}{2 p_{\lim}}\right)$$

  • Điểm đột phá nằm ở việc xây dựng trường tải trọng chủ động động đất tương đương tác dụng lên các phân đoạn biên của hầm chữ nhật cong thông qua bộ tham số hiệu chỉnh vô thứ nguyên ($a, b, \beta$). Các tham số $\beta_1, \beta_2, \beta_3, \beta_4$ được xây dựng tương ứng cho từng cung bán kính ($R_1, R_2, R_3$), giúp sơ đồ tải trọng tự động thích ứng với sự biến thiên hình học của mặt cắt ngang.

Điều kiện biên lý thuyết xác định rõ: Môi trường đất đồng nhất, đẳng hướng, ứng xử đàn hồi tuyến tính; vỏ hầm làm việc trong giai đoạn đàn hồi chưa nứt; chuyển vị và biến dạng giả định ở trạng thái biến dạng phẳng (plane strain).

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (positivism paradigm) với quy trình diễn dịch kết hợp giữa mô phỏng số học liên tục chuyên sâu và phát triển mô hình giải thuật phần tử hữu hạn rút gọn.

Thiết kế nghiên cứu đa tầng (multi-level design) bao gồm 3 cấp độ:

  • Cấp độ 1 (Kiểm chuẩn nền tảng): Xây dựng mô hình số sai phân hữu hạn 2D FDM (FLAC) cho hầm tròn chịu tải trọng tĩnh tương đương, so sánh trực tiếp với nghiệm giải tích giải tích của Wang (1993) cho hai điều kiện biên tiếp xúc (full-slip và no-slip) nhằm chuẩn hóa lưới phần tử và điều kiện biên.
  • Cấp độ 2 (Phân tích tham số số học chuyên sâu): Thiết lập mô hình số FDM hầm chữ nhật cong (tham chiếu kích thước từ hầm metro Thượng Hải: chiều rộng $w$, chiều cao $h$), mô phỏng trường ứng suất cắt tự do thông qua các lực cắt tác dụng tại biên mô hình $\tau$. Thực hiện khảo sát tham số đa biến diện rộng.
  • Cấp độ 3 (Thiết lập và Chuẩn hóa giải thuật HRM mới): Tích hợp sơ đồ tải trọng tĩnh tương đương vào thuật toán HRM, tối ưu hóa các hàm hệ số tải trọng ($\beta, a, b$) thông qua phương pháp khớp hàm phi tuyến (nonlinear curve fitting) đối chiếu với dữ liệu FDM.
                    QUY TRÌNH NGHIÊN CỨU RIGOROUS

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu số học được thực hiện nghiêm ngặt:

  1. Thiết lập miền tính toán: Mô hình FDM có kích thước đủ lớn ($100\text{ m} \times 100\text{ m}$) để loại bỏ ảnh hưởng phản xạ biên. Kích thước phần tử nhỏ nhất $\Delta z_{\min}$ tại lân cận vỏ hầm được kiểm soát chặt chẽ nhằm đảm bảo độ chính xác của gradient ứng suất tiếp xúc.
  2. Mô phỏng trường cắt địa chấn: Trường biến dạng cắt tự do cực đại $\gamma_{\max}$ được tính toán từ vận tốc sóng cắt đỉnh $V_{\max}$ và mô đun trượt của đất $G_{\max} = E_s / [2(1 + \nu_s)]$. Ứng suất cắt biên tương đương $\tau = G_{\max} \cdot \gamma_{\max}$ được áp đặt lên đỉnh và hai bên biên mô hình dưới điều kiện biến dạng phẳng.
  3. Mô hình hóa tiếp xúc (Interface Modeling): Giao diện giữa đất và vỏ hầm được thiết lập bằng các phần tử tiếp xúc đặc biệt (interface elements) với hai trạng thái cực hạn:
    • Full-slip: Độ cứng tiếp tuyến $k_s = 0$, ma sát tiếp xúc bằng 0.
    • No-slip: Đất và vỏ hầm liên kết chặt chẽ, không có chuyển vị trượt tương đối.
  4. Kiểm tra độ tin cậy và tính hội tụ: Quá trình lặp số FDM sử dụng tiêu chí hội tụ tỷ số mất cân bằng lực cực đại nhỏ hơn $10^{-5}$. Giải thuật HRM được lập trình tự động hóa, kiểm tra tính độc lập của bước chia phần tử thanh dầm (tối thiểu 60-120 phần tử trên chu vi hầm).

Data và phân tích

Dữ liệu tham số đầu vào được trích xuất từ các dự án thực tế và tiêu chuẩn quốc tế:

  • Thông số đất nền: Dung trọng $\gamma = 18\text{ - }20\text{ kN/m}^3$; Mô đun đàn hồi $E_s = 20, 50, 100, 200\text{ MPa}$; Hệ số Poisson $\nu_s = 0.3$; Hệ số áp lực hông $K_0 = 0.5\text{ - }0.7$.
  • Thông số kết cấu vỏ hầm: Bê tông cốt thép C30/C40 có mô đun đàn hồi $E = 30\text{ - }35\text{ GPa}$, hệ số Poisson $\nu = 0.15\text{ - }0.20$, chiều dày thay đổi $t = 0.30, 0.45, 0.60\text{ m}$.
  • Thông số địa chấn: Gia tốc đỉnh bề mặt $a_H = 0.1g, 0.2g, 0.3g, 0.4g$, suy giảm theo độ sâu theo bảng chuẩn hóa của Power et al. (1996) và Sadigh & Egan (1998).

Kỹ thuật phân tích dữ liệu ứng dụng thuật toán hồi quy phi tuyến đa biến để xác định quy luật biến thiên của 4 hệ số hình học $\beta_1, \beta_2, \beta_3, \beta_4$ định hình cho hàm tải trọng $\beta(\theta)$, cùng các hệ số cường độ $a$ và $b$. Độ sai lệch tương đối giữa HRM và FDM được đánh giá qua sai số toàn phương trung bình (Root Mean Square Error - RMSE) và hệ số xác định ($R^2 > 0.96$).

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

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

  1. Sự tập trung dị thường của mô men uốn tại góc vai vòm ($R_2$): Khác với hầm tròn có phân bố mô men cực trị tại các góc $\theta = 45^\circ, 135^\circ, 225^\circ, 315^\circ$, hầm chữ nhật cong ghi nhận mô men uốn cực đại $M_{\max}$ dịch chuyển mạnh về phía cung vai vòm ($R_2$). Tại vị trí này, độ cong thay đổi đột ngột làm tăng mô men gia tăng lên từ $28%$ đến $45%$ so với mặt cắt tròn có cùng diện tích sử dụng không gian ngầm dưới cùng cấp gia tốc $a_H = 0.2g$.
  2. Quy luật chi phối tuyệt đối của điều kiện tiếp xúc đối với lực dọc ($N_{\max}$): Dữ liệu phân tích chứng minh rằng điều kiện no-slip tạo ra lực dọc nén gia tăng $N_{\max}$ lớn hơn từ $2.5$ đến $3.8$ lần so với điều kiện full-slip. Điều này khẳng định ma sát tiếp xúc đất - vỏ hầm đóng vai trò truyền ứng suất tiếp tuyến rất lớn trong suốt quá trình truyền sóng cắt.
  3. Mối quan hệ phi tuyến giữa độ cứng đất $E_s$ và nội lực động đất: Khi mô đun đàn hồi của đất tăng từ $20\text{ MPa}$ lên $200\text{ MPa}$ (đất chuyển từ mềm sang cứng), hệ số độ mềm $F$ giảm mạnh, dẫn đến lực dọc gia tăng $N_{\max}$ tăng gần như tuyến tính nhưng mô men uốn gia tăng $M_{\max}$ có xu hướng bão hòa và giảm biên độ biến dạng tương đối do đất cứng cản trở chuyển vị méo của hầm.
  4. Hiệu ứng phi đối xứng của chiều dày vỏ hầm $t$: Tăng chiều dày vỏ hầm từ $0.3\text{ m}$ lên $0.6\text{ m}$ làm tăng mô men quán tính $I$ gấp 8 lần, khiến độ cứng chống uốn $EI$ tăng vọt. Kết quả là mô men uốn động đất gia tăng $M_{\max}$ tăng hơn $180%$, trong khi lực dọc $N_{\max}$ chỉ tăng nhẹ ($15\text{ - }25%$). Điều này đưa ra cảnh báo quan trọng: Việc tăng chiều dày vỏ hầm để gia cường chịu lực tĩnh có thể phản tác dụng, làm tăng mức độ tổn thương do uốn dưới tải trọng động đất.
  5. Độ chính xác vượt trội của sơ đồ tải trọng HRM mới: So sánh kết quả tính toán giữa HRM và FDM trên 12 kịch bản hình học và địa chất khác nhau cho thấy sai số cực trị mô men uốn và lực dọc không vượt quá $6.8%$, trong khi thời gian thực thi thuật toán giảm từ hàng chục phút (trên FDM liên tục) xuống dưới 2 giây (trên mô hình HRM 1D).
So sánh phân bố Mô men uốn gia tăng (Delta M) và Lực dọc (Delta N):
- Hầm tròn: Biến dạng Ovaling đối xứng hoàn hảo 4 đỉnh cực trị tại góc 45°, 135°, 225°, 315°.
- Hầm chữ nhật cong: Biến dạng Sub-ovaling bất đối xứng, đỉnh cực trị dồn cục bộ tại vòm vai (R2) và chân vách (R3).

Implications đa chiều

  • Về mặt lý thuyết: Nghiên cứu mở rộng hệ thống biên lý thuyết tương tác kết cấu - môi trường từ hình học đối xứng tâm sang hình học đa bán kính cong, cung cấp bằng chứng thực nghiệm số cho các lý thuyết cơ học môi trường liên tục chịu tải trọng biến dạng cưỡng bức.
  • Về mặt phương pháp luận: Thiết lập thành công quy trình chuẩn hóa các tham số vô thứ nguyên trong phương pháp HRM, chứng minh phương pháp này hoàn toàn có khả năng mở rộng sang các dạng kết cấu ngầm phức tạp khác như hầm vòm móng ngựa (horseshoe), hầm đa cung tròn (multi-circular face - MF), hoặc hầm đôi giao nhau (DOT).
  • Về ứng dụng thực tiễn công trình: Cung cấp biểu đồ thiết kế nhanh cho kỹ sư kết cấu ngầm. Khuyến cáo trực tiếp các nhà thiết kế không được áp dụng nguyên mẫu công thức tính hầm tròn cho hầm chữ nhật cong; khu vực vòm vai ($R_2$) bắt buộc phải được bố trí cốt thép chịu uốn tăng cường hoặc sử dụng bê tông sợi tính năng cao.
  • Về mặt chính sách và quy chuẩn: Cung cấp luận cứ khoa học để bổ sung vào các tiêu chuẩn thiết kế công trình ngầm chịu động đất (như TCVN, Eurocode 7 & 8, FHWA), đưa ra các hướng dẫn định lượng về việc lựa chọn điều kiện trượt tiếp xúc trong tính toán an toàn hạ tầng ngầm đô thị.

Limitations và Future Research

Nghiên cứu công nhận một số giới hạn học thuật cần được xem xét cẩn trọng:

  1. Mô hình vật liệu tuyến tính: Cả môi trường đất và kết cấu bê tông vỏ hầm đều giả định làm việc trong miền đàn hồi tuyến tính (linearly elastic), chưa xét đến ứng xử dẻo phi tuyến cục bộ của đất nền hoặc sự suy giảm độ cứng do nứt nẻ bê tông sau nhiều chu kỳ tải trọng động đất lặp lại.
  2. Giả định trạng thái tĩnh tương đương 2D: Mô hình phân tích dựa trên trạng thái biến dạng phẳng 2D dưới tác động của sóng cắt đơn hướng, chưa mô phỏng được hiệu ứng không gian 3D của sóng địa chấn lan truyền xiên góc hoặc các thành phần sóng mặt (Rayleigh, Love) gây uốn và nén dọc trục đường hầm.
  3. Bỏ qua khớp nối phân đoạn (Segmental Joints): Mô hình xem xét vỏ hầm như một kết cấu đúc liền khối liên tục (continuous lining), chưa xét đến sự suy giảm độ cứng uốn cục bộ tại các vị trí bu-lông nối phân đoạn trong vỏ hầm lắp ghép TBM thực tế.
  4. Điều kiện địa chất đồng nhất: Đất nền được giả định đồng nhất, chưa xét đến các địa tầng phân lớp phức tạp, hiện tượng mực nước ngầm dâng cao hoặc nguy cơ hóa lỏng cát (soil liquefaction).

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 5 hướng tiếp cận:

  • Mở rộng mô hình HRM cho vỏ hầm phân đoạn lắp ghép: Tích hợp các phần tử khớp lò xo quay phi tuyến (rotational spring joints) để mô phỏng chính xác ứng xử của bu-lông và gioăng cao su truyền lực.
  • Phân tích động lực học 3D phi tuyến đầy đủ: Thực hiện mô phỏng lịch sử thời gian (time-history dynamic FDM/FEM) với trường kích động địa chấn 3 pha không gian thực.
  • Nghiên cứu tương tác trong môi trường đất hóa lỏng: Đánh giá nguy cơ nổi hầm (flotation) và biến dạng cắt lớn khi hầm chữ nhật cong nằm trong tầng cát bão hòa nước chịu chấn động mạnh.
  • Thực nghiệm kiểm chứng bàn rung hoặc mô hình ly tâm đa phương: Chế tạo mô hình vật lý hầm chữ nhật cong thu nhỏ để kiểm chứng trực tiếp cơ chế "sub-ovaling" và ứng suất tập trung tại góc vai.
  • Tích hợp trí tuệ nhân tạo (AI/Machine Learning) trong tối ưu hóa đa mục tiêu: Ứng dụng mạng nơ-ron sâu kết hợp HRM để tự động tìm kiếm hình học mặt cắt chữ nhật cong tối ưu, tối thiểu hóa đồng thời chi phí đào đất và nội lực động đất.

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

Nghiên cứu mang lại những ảnh hưởng sâu rộng trên nhiều phương diện:

  • Tác động học thuật (Academic Impact): Mở ra một nhánh nghiên cứu chuyên sâu về cơ học kháng chấn cho kết cấu ngầm phi tròn. Ước tính công trình tạo nguồn trích dẫn quan trọng cho các nghiên cứu thuộc nhóm Q1/Q2 về Kỹ thuật Xây dựng ngầm (Tunnelling and Underground Space Technology, Computers and Geotechnics, Soil Dynamics and Earthquake Engineering).
  • Chuyển đổi công nghệ công nghiệp (Industry Transformation): Thúc đẩy ứng dụng công nghệ TBM phi tròn tại các dự án metro đô thị. Bằng việc cung cấp công cụ tính toán HRM chính xác và siêu nhanh, các đơn vị tư vấn thiết kế có thể rút ngắn quy trình thiết kế sơ bộ từ nhiều tuần xuống vài ngày, tối ưu hóa vật liệu vỏ hầm và tiết kiệm từ $10\text{ - }15%$ chi phí xây dựng kết cấu.
  • Tác động chính sách và hạ tầng quốc gia (Policy & Resilience): Đóng góp trực tiếp vào mục tiêu phát triển hạ tầng bền vững, nâng cao khả năng chống chịu động đất của mạng lưới giao thông ngầm tại các quốc gia nằm trên vành đai địa chấn phức tạp (như Việt Nam, Nhật Bản, Đài Loan).
  • Lợi ích xã hội và môi trường (Societal Benefits): Việc tối ưu hóa hình học hầm chữ nhật cong giúp giảm khối lượng đào thải đất đá đô thị, hạn chế tối đa ảnh hưởng lún nứt công trình lân cận trên mặt đất trong quá trình thi công và vận hành an toàn qua các trận động đất lịch sử.

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

  • Nghiên cứu sinh và Giới học thuật (Doctoral & Senior Researchers): Tiếp cận khung lý thuyết hoàn chỉnh về cơ chế "sub-ovaling", bộ dữ liệu kiểm chuẩn tin cậy giữa FDM và HRM, cùng các công thức giải tích mở rộng để tiếp tục phát triển các mô hình cơ học tiên tiến.
  • Kỹ sư tư vấn thiết kế công trình ngầm (Design Engineers): Nhận được công cụ tính toán tĩnh tương đương HRM đã được đóng gói và kiểm chứng, cho phép đánh giá nhanh chóng nội lực $\Delta M, \Delta N$ mà không cần trang bị các phần mềm động lực học 3D đắt tiền và phức tạp.
  • Các nhà thầu thi công TBM và Chế tạo vỏ hầm (TBM Manufacturers & Contractors): Nắm bắt được phân bố ứng suất nguy hiểm tại vòm vai để tối ưu hóa thiết kế cốt thép, vị trí bố trí khớp nối phân đoạn và quy trình bơm vữa chèn lưng vỏ hầm (grouting).
  • Cơ quan quản lý quy hoạch đô thị và Giao thông (Policy Makers & Urban Authorities): Có cơ sở khoa học vững chắc để phê duyệt các dự án hầm giao thông tiết diện đặc biệt, đảm bảo tính an toàn dài hạn của mạng lưới không gian ngầm đô thị trước rủi ro thiên tai.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của nghiên cứu là gì và lý thuyết nào đã được mở rộng?

Đóng góp lý thuyết độc đáo nhất là việc phát hiện và hình thức hóa toán học cơ chế biến dạng "sub-ovaling" cùng sự tập trung nội lực cục bộ tại vòm vai ($R_2$) của hầm chữ nhật cong dưới tác dụng của sóng cắt địa chấn. Nghiên cứu đã mở rộng lý thuyết tương tác đàn hồi kinh điển của Wang (1993) và Kouretzis et al. (2006) – vốn chỉ áp dụng cho hầm tròn đối xứng – sang môi trường hình học đa bán kính cong phức tạp, đồng thời mở rộng phạm vi áp dụng của phương pháp phản lực đàn hồi (HRM) từ tĩnh học sang động lực học tĩnh tương đương.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu trước đây là gì?

So với nghiên cứu của Do et al. (2014) và Du et al. (2021) (vốn chỉ áp dụng HRM cho bài toán đào hầm chịu tải trọng tĩnh) và nghiên cứu của Sederat et al. (2009) (chỉ dùng FEM phân tích hầm tròn), luận án đã:

  • Tích hợp thành công sơ đồ phân bố tải trọng động đất chủ động và hệ số phản lực nền lò xo biến thiên theo hàm góc phương vị $\theta$ vào giải thuật HRM.
  • Xây dựng quy trình chuẩn hóa các tham số hình học vô thứ nguyên ($\beta, a, b$) cho phép dự báo nội lực với sai số dưới $6.8%$ so với mô hình liên tục FDM (FLAC), nhưng giảm thời gian tính toán hàng trăm lần.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu chứng minh?

Phát hiện bất ngờ nhất là xu hướng gia tăng nghịch lý của mô men uốn khi tăng chiều dày vỏ hầm $t$. Khi tăng $t$ từ $0.3\text{ m}$ lên $0.6\text{ m}$, mô men uốn gia tăng $M_{\max}$ tăng vọt hơn $180%$, trong khi lực dọc $N_{\max}$ chỉ biến thiên từ $15\text{ - }25%$. Nguyên nhân là việc tăng chiều dày làm tăng độ cứng chống uốn $EI$ theo hàm bậc 3 của $t$, khiến kết cấu hút nhiều năng lượng biến dạng uốn cưỡng bức hơn từ môi trường đất xung quanh trong quá trình sóng cắt truyền qua.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) rõ ràng không?

Có. Luận án cung cấp đầy đủ thông số đầu vào hình học ($R_1, R_2, R_3, w, h, t$), đặc tính cơ lý đất nền ($E_s, \nu_s, \gamma, K_0$), đặc tính kết cấu ($E, \nu$), quy trình thiết lập điều kiện biên FDM trong FLAC, các công thức giải tích xác định ứng suất cắt biên $\tau = G_{\max} \cdot \gamma_{\max}$, và toàn bộ lưu đồ thuật toán hiệu chuẩn tham số cho phương pháp HRM.

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

Lộ trình 10 năm tập trung vào 4 giai đoạn chiến lược:

  • Giai đoạn 1 (1-2 năm): Phát triển mô hình HRM phi tuyến cho vỏ hầm phân đoạn lắp ghép có khớp nối bu-lông chịu động đất.
  • Giai đoạn 2 (3-5 năm): Thực hiện thí nghiệm mô hình ly tâm động lực học để kiểm chứng các hệ số tải trọng tĩnh tương đương dưới sóng kích động ngẫu nhiên.
  • Giai đoạn 3 (6-8 năm): Xây dựng tiêu chuẩn thiết kế kháng chấn chuyên biệt cho hầm TBM phi tròn và tích hợp vào quy chuẩn xây dựng đô thị ngầm.
  • Giai đoạn 4 (9-10 năm): Phát triển phần mềm thương mại hóa tối ưu hóa tự động thiết kế kháng chấn kết cấu ngầm đô thị ứng dụng trí tuệ nhân tạo.

Kết luận

  1. Lần đầu tiên làm sáng tỏ toàn diện ứng xử kháng chấn của hầm chữ nhật cong: Nghiên cứu đã mô tả chi tiết cơ chế biến dạng "sub-ovaling" và sự phân bố nội lực gia tăng ($\Delta M, \Delta N$) dưới tác động của sóng cắt địa chấn lan truyền vuông góc, khắc phục hoàn toàn khoảng trống lý thuyết của các nghiên cứu trước đây.
  2. Khẳng định vai trò quyết định của điều kiện tiếp xúc và hình học cục bộ: Chứng minh bằng số liệu định lượng rằng điều kiện no-slip tạo ra lực dọc gia tăng lớn nhất (gấp $2.5\text{ - }3.8$ lần so với full-slip), và góc vai vòm ($R_2$) là vị trí xung yếu nhất chịu mô men uốn tập trung cao hơn từ $28\text{ - }45%$ so với hầm tròn cùng diện tích.
  3. Phát triển thành công sơ đồ tải trọng tĩnh tương đương mới cho phương pháp HRM: Thiết lập hệ thống công thức giải tích hiệu chỉnh thông qua bộ tham số vô thứ nguyên ($\beta, a, b$), cho phép phương pháp HRM tính toán nội lực động đất với độ chính xác cao (sai số $< 6.8%$ so với FDM) và tốc độ tính toán tức thời.
  4. Cung cấp chỉ dẫn thực tiễn đột phá cho thiết kế kết cấu ngầm: Đưa ra cảnh báo khoa học về việc gia tăng chiều dày vỏ hầm sẽ làm tăng mạnh mô men uốn động đất, đồng thời cung cấp các giải pháp bố trí cốt thép và cấu tạo khớp nối hợp lý cho khu vực góc vai hầm.
  5. Mở ra các nhánh nghiên cứu học thuật mới trên trường quốc tế: Tạo tiền đề vững chắc cho việc phát triển các mô hình tính toán kháng chấn cho hầm lắp ghép phân đoạn, hầm trong đất hóa lỏng và tối ưu hóa không gian ngầm đô thị bền vững trong tương lai.