Tổng quan về luận án

Trong kỹ thuật xây dựng và công trình giao thông hiện đại, tương tác kết cấu - móng - đất nền (Soil-Structure Interaction - SSI / Soil-Foundation-Structure Interaction - SFSI) dưới tác dụng của tải trọng động đất là một bài toán cơ học phi tuyến phức tạp nhưng có ý nghĩa quyết định đến độ an toàn công trình. Luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (Mã số: 9520101) của nghiên cứu sinh Huỳnh Văn Quân, được thực hiện dưới sự hướng dẫn khoa học của PGS. Nguyễn Xuân Huy và PGS. Nguyễn Trung Kiên tại Trường Đại học Giao thông Vận tải (Hà Nội, 2021), mang tiêu đề: "Nghiên cứu tương tác kết cấu-đất nền dưới tác dụng của tải trọng động đất bằng phần tử vĩ mô". Đây là một trong những công trình tiên phong tại Việt Nam tiếp cận bài toán động lực học công trình chịu động đất thông qua phương pháp lai (hybrid method) với mô hình phần tử vĩ mô (macro-element).

Thực tiễn thiết kế công trình kháng chấn hiện hành bộc lộ một khoảng trống học thuật và thực nghiệm rất lớn. Trong các bộ tiêu chuẩn thiết kế phổ biến như tiêu chuẩn thiết kế cầu đường bộ TCVN 11823:2017, 22TCN 272-05, AASHTO LRFD (Hoa Kỳ), tiêu chuẩn thiết kế công trình chịu động đất TCVN 9386:2012 hay Eurocode 8 (EN 1998), việc phân tích ứng xử của kết cấu có xét đến tương tác với đất nền hầu như chưa được kể đến một cách bắt buộc hoặc chỉ dừng lại ở mức độ khuyến nghị định tính. Quy trình thiết kế truyền thống luôn bị chia cắt cơ học: kỹ sư kết cấu giả định chân công trình là các liên kết ngàm cứng hoặc gối cố định để xuất phản lực chân cột; ngược lại, kỹ sư địa kỹ thuật lấy các nội lực tĩnh/động này để tính toán móng trên nền đàn hồi tuyến tính. Sự chia tách này hoàn toàn bỏ qua tương tác động (kinematic interaction) và tương tác quán tính (inertial interaction) diễn ra đồng thời trong quá trình rung chấn.

Nghiên cứu của luận án tập trung giải quyết ba câu hỏi nghiên cứu cốt lõi (Research Questions - RQ) và kiểm chứng hai giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để xây dựng một mô hình toán học tích hợp thu gọn (macro-element) có khả năng mô tả chính xác phản ứng phi tuyến của hệ móng - đất nền chịu động đất mà không đòi hỏi tài nguyên tính toán khổng lồ như phương pháp phần tử hữu hạn (FEM) 3D trực tiếp?
  • RQ2: Cơ chế tương tác phi tuyến hình học (hiện tượng nhổ móng - uplift) và phi tuyến vật liệu (tính dẻo hóa của đất nền) tác động như thế nào đến phản ứng gia tốc và chuyển vị ngang của kết cấu phần trên?
  • RQ3: Mức độ tin cậy và sai số của mô hình phần tử vĩ mô đề xuất khi so sánh đối chứng trực tiếp với thực nghiệm bàn rung vật lý thu nhỏ là bao nhiêu?
  • H1: Việc tích hợp đồng thời cơ chế nhổ móng đàn dẻo trong phần tử vĩ mô sẽ làm giảm đáng kể gia tốc phản ứng của kết cấu phần trên so với mô hình ngàm cứng truyền thống nhờ hiệu ứng tiêu tán năng lượng và kéo dài chu kỳ dao động riêng.
  • H2: Phương pháp lai sử dụng phần tử vĩ mô kết hợp mô hình thông số tập trung (Lumped-Parameter Model - LPM) giải bằng tích phân số Newmark trên nền tảng MATLAB có thể tái hiện chính xác ứng xử động lực học thực tế với thời gian tính toán giảm trên 90% so với phương pháp phần tử hữu hạn trực tiếp.

Khung lý thuyết của luận án tích hợp lý thuyết tương tác kết cấu - đất nền hiện đại của Kausel (1978), mô hình dẻo hóa móng nông của Nova (1991), mô hình phi tuyến hình học của Cremer (2001, 2002) và Paolucci (1997, 2008). Về phạm vi nghiên cứu, luận án khảo sát hệ móng nông đặt trên nền cát chặt chịu kích thích động đất theo phương ngang (sử dụng các giản đồ gia tốc động đất thực tế như Tolmezzo, El-Centro 1940), kết hợp kiểm chứng thực nghiệm trên bàn rung R202(UTC) với 4 kịch bản hệ móng - đất nền (T12, T13, T14, T15) và 6 kịch bản hệ kết cấu phần trên - móng - đất nền (T21 đến T26) trong các điều kiện chôn sâu móng khác nhau (0%, 50%, 100% chiều cao móng).

Literature Review và Positioning

Lý thuyết phân tích tương tác kết cấu - đất nền (SSI) trên thế giới phát triển qua ba dòng phương pháp luận chính:

                          ┌────────────────────────────────────────────────────────┐
                          │   CÁC PHƯƠNG PHÁP PHÂN TÍCH TƯƠNG TÁC KẾT CẤU-ĐẤT NỀN  │
                          └───────────────────────────┬────────────────────────────┘
                                                      │
         ┌────────────────────────────────────────────┼──────────────────────────────────────────┐
         ▼                                            ▼                                          ▼
┌─────────────────────────────────┐         ┌─────────────────────────────────┐        ┌─────────────────────────────────┐
│     PHƯƠNG PHÁP TRỰC TIẾP       │         │    PHƯƠNG PHÁP KẾT CẤU PHỤ      │        │        PHƯƠNG PHÁP LAI          │
│        (Direct FEM)             │         │     (Sub-structure Method)      │        │         (Hybrid Method)         │
├─────────────────────────────────┤         ├─────────────────────────────────┤        ├─────────────────────────────────┤
│ • Mô hình hóa đồng thời 3D.     │         │ • Đề xuất: Kausel (1978).       │        │ • Kết hợp LPM + Macro-element.  │
│ • Xét đầy đủ phi tuyến vật liệu.│         │ • Nguyên lý cộng tác dụng.      │        │ • Nova (1991), Cremer (2001),   │
│ • Nhược điểm: Khối lượng tính   │         │ • Đơn giản, dùng lò xo đàn hồi. │          Paolucci (2008), Grange (2009). │
│   toán khổng lồ, khó kiểm soát  │         │ • Nhược điểm: Bắt buộc giả định │        │ • Ưu điểm: Phản ánh đầy đủ      │
│   điều kiện biên vô hạn.        │         │   ứng xử tuyến tính đàn hồi.    │          phi tuyến, tính toán tối ưu.    │
└─────────────────────────────────┘         └─────────────────────────────────┘        └─────────────────────────────────┘
  1. Phương pháp trực tiếp (Direct Method): Mô hình hóa đồng thời toàn bộ môi trường đất và kết cấu bằng các phần tử hữu hạn liên tục (Datta, 2010). Dù phản ánh chi tiết các trạng thái ứng suất - biến dạng cục bộ, phương pháp này đòi hỏi chi phí tính toán cực lớn, dễ gặp lỗi điều kiện biên sóng phản xạ tại biên nhân tạo và phụ thuộc nặng nề vào năng lực mô hình hóa của người kỹ sư.
  2. Phương pháp kết cấu phụ (Sub-structure Method): Được Kausel (1978) đề xuất và phát triển bởi John (1985), Pais (1988), Gazetas (1991), Mylonakis (2006), Chatzigogos (2008). Phương pháp này tách bài toán thành ba bước: xác định chuyển động động học móng không khối lượng, tính toán độ cứng lò xo đất tương đương (soil springs) và phân tích phản ứng kết cấu. Điểm hạn chế chí tử của phương pháp này là dựa trên nguyên lý cộng tác dụng (superposition theorem), do đó bắt buộc phải giả định hệ làm việc tuyến tính đàn hồi, không thể phản ánh được các biến dạng dẻo và hiện tượng tách rời móng khi động đất mạnh.
  3. Phương pháp lai (Hybrid Method) và Phần tử vĩ mô (Macro-element): Khởi nguồn từ ý tưởng của Nova (1991) khi nghiên cứu độ lún móng nông dưới tải trọng lệch tâm và nghiêng, phương pháp lai thay thế toàn bộ hệ móng - đất nền phức tạp bằng một phần tử tập trung duy nhất tại chân công trình. Phần tử này tích hợp cả vùng gần (near-field: chứa các phi tuyến hình học và dẻo hóa vật liệu) và vùng xa (far-field: tính chất đàn hồi tuyến tính và cản lan truyền sóng vô hạn).

Về mặt học thuật, tồn tại một cuộc tranh luận kéo dài giữa hai quan điểm thiết kế: kỹ sư kết cấu coi móng là ngàm cứng tuyệt đối và kết cấu như một lá mỏng linh hoạt (Grange, 2008), trong khi kỹ sư địa kỹ thuật coi kết cấu là một khối cứng vô hạn đặt trên gối mềm đàn hồi (Kotronis, 2013). Khi động đất xảy ra, gia tốc truyền từ nền đất kích hoạt lực quán tính trong kết cấu. Lực quán tính sinh ra mô men lật và lực cắt tại đáy móng. Nếu vượt qua khả năng chịu kéo của mặt tiếp xúc, hiện tượng nhổ móng (uplift) sẽ xuất hiện, làm thay đổi đột ngột diện tích tiếp xúc hiệu dụng (phi tuyến hình học) đồng thời gây tập trung ứng suất dẻo tại mép móng (phi tuyến vật liệu) (Anastasopoulos, 2012; Paolucci, 2008). Hiện tượng này đã được quan sát thực tế trong các thảm họa động đất lớn như Arvin Tehachapi (1952), Alaska (1964), Imperial Valley (1979) (Psycharis, 1983).

Trong dòng nghiên cứu về macro-element, các tác giả quốc tế đã có những đóng góp quan trọng: Cremer (2001, 2002) phát triển mô hình cho móng băng chịu tải trọng chu kỳ và ứng dụng thành công cho cầu Viaduc de l'Arc (Marseille, Pháp); Paolucci (1997, 2008) đưa macro-element vào phân tích kết cấu 4 bậc tự do dưới tải trọng động đất, kiểm chứng với dự án PWRI (Nhật Bản); Grange (2009) mở rộng phần tử vĩ mô 3D không gian cho móng tròn và móng chữ nhật, kiểm chứng trên dự án Camus IV (Pháp); Chatzigogos (2009, 2011) bổ sung cơ chế trượt tiếp xúc; Figini (2012) đề xuất mô hình suy giảm độ cứng kiểm chứng trên các dự án Camus IV, TRISEE-1 và TRISEE-2.

Tại Việt Nam, các nghiên cứu tương tác đất - kết cấu trước đây (Lê Văn Tuân 2016, Đào Văn Hưng 2017, Vũ Minh Ngọc 2019) hầu như chỉ dừng lại ở mô hình lò xo tuyến tính đàn hồi hoặc công trình ngầm. Luận án của Huỳnh Văn Quân đã định vị chính xác khoảng trống khoa học này: lần đầu tiên tại Việt Nam xây dựng hoàn chỉnh thuật toán phần tử vĩ mô phi tuyến đàn dẻo - nhổ móng cho bài toán động đất, lập trình giải số tự động và tiến hành chuỗi thí nghiệm bàn rung vật lý quy mô bài bản để kiểm chứng lý thuyết.

Đó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 khung lý thuyết cơ học đàn dẻo của Nova (1991) và lý thuyết nhổ móng phi tuyến của Cremer (2001), Paolucci (2008), thiết lập nên một mô hình toán học hoàn chỉnh cho phần tử vĩ mô phẳng trong hệ tọa độ tổng quát:

  1. Chuẩn hóa lực và chuyển vị thu gọn: Véc-tơ lực thu gọn không thứ nguyên $\mathbf{F}'$ và véc-tơ chuyển vị thu gọn $\mathbf{u}'$ được định nghĩa: $$\mathbf{F}' = \left{ \frac{H}{\mu N_{max}}, \frac{M}{\psi B N_{max}}, \frac{N}{N_{max}} \right}^T, \quad \mathbf{u}' = \left{ u_x, B\theta, u_z \right}^T$$ Trong đó $N_{max}$ là sức chịu tải cực hạn thẳng đứng đúng tâm theo Terzaghi: $N_{max} = (\gamma D N_q + 0.5 \gamma B N_\gamma) B L$; $\mu$ là hệ số ma sát trượt đất - móng ($\mu = \tan\phi$); $\psi$ là hệ số điều chỉnh vị trí tải trọng cực đại ($\psi \approx 0.35 \div 0.5$); $B$ là bề rộng móng.

  2. Mặt dẻo hóa (Yield Surface) và Quy luật chảy (Flow Rule): Hàm dẻo $f(\mathbf{F}') = 0$ được mô tả dưới dạng ellipsoid phi tuyến bao bọc miền đàn hồi: $$f(\mathbf{F}') = \left( \frac{H'}{1 - N'} \right)^2 + \left( \frac{M'}{1 - N'} \right)^2 - \beta^2 N'^2 (1 - N')^2 = 0$$ Khi trạng thái ứng suất chạm mặt dẻo ($f \ge 0$), ma trận độ cứng đàn dẻo tức thời $\mathbf{K}{ep}$ được cập nhật thông qua quy luật chảy trực giao hoặc phi trực giao $g(\mathbf{F}')$: $$\mathbf{K}{ep} = \mathbf{K}_e - \frac{\mathbf{K}_e \frac{\partial g}{\partial \mathbf{F}'} \left( \frac{\partial f}{\partial \mathbf{F}'} \right)^T \mathbf{K}_e}{H_p + \left( \frac{\partial f}{\partial \mathbf{F}'} \right)^T \mathbf{K}_e \frac{\partial g}{\partial \mathbf{F}'}}$$

  3. Cơ chế nhổ móng động (Dynamic Uplift Coupling): Khi mô men lật vượt ngưỡng tới hạn $M_{u,0} = \pm \frac{N B}{2} (1 - \frac{N}{N_{max}})$, diện tích tiếp xúc bị thu hẹp thành $B' = B(1 - \xi)$. Hệ số diện tích tách rời $\xi(t)$ được cập nhật liên tục theo lịch sử biến dạng góc xoay tích lũy: $$\xi(t) = \frac{\xi_1 \theta_{el}}{1 + \xi_2 \theta_{el}}$$ Ma trận độ cứng tiếp tuyến của hệ được tái cấu trúc phi tuyến, cho phép mô tả chính xác sự chuyển tiếp mượt mà giữa pha đàn hồi, pha tách rời hình học và pha dẻo hóa vật liệu đất nền.

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 lý thuyết:

  • Lý thuyết biến dạng dẻo tăng bền đẳng hướng (Nova, 1991).
  • Lý thuyết tiếp xúc phi tuyến hình học động lực học (Paolucci, 2008; Chatzigogos, 2009).
  • Động lực học công trình hệ thông số tập trung (LPM - Clough & Penzien, Chopra).
┌─────────────────────────────────────────────────────────────────────────────┐
│                          HỆ KẾT CẤU PHẦN TRÊN                               │
│      Mô hình hóa thông số tập trung (LPM) - Bậc tự do kết cấu [M_s, K_s]     │
└──────────────────────────────────────┬──────────────────────────────────────┘
                                       │ (Chuyển vị u_s, Gia tốc a_s)
                                       ▼
┌─────────────────────────────────────────────────────────────────────────────┐
│                       NÚT TƯƠNG TÁC PHẦN TỬ VĨ MÔ                            │
│  Lực thu gọn: [H, M, N]  <───────────────>  Chuyển vị thu gọn: [u_x, θ, u_z] │
├──────────────────────────────────────┬──────────────────────────────────────┤
│  VÙNG GẦN (Near-field Non-linear)    │   VÙNG XA (Far-field Linear)         │
│  • Phi tuyến hình học: Nhổ móng ξ(t) │   • Độ cứng đàn hồi tĩnh lò xo đất   │
│  • Phi tuyến vật liệu: Mặt dẻo f(F') │   • Cản lan truyền vô hạn (Radiation)│
└──────────────────────────────────────┴──────────────────────────────────────┘
                                       ▲
                                       │ (Sóng động đất kích thích u_g, a_g)
┌──────────────────────────────────────┴──────────────────────────────────────┐
│                               ĐẤT NỀN TỰ DO                                 │
└─────────────────────────────────────────────────────────────────────────────┘

Công trình đưa ra điều kiện biên áp dụng rõ ràng (boundary conditions): móng nông chữ nhật hoặc băng đặt trên nền đất hạt rời đồng chất, tỷ số mảnh kết cấu thỏa mãn điều kiện tương tác SSI theo chuẩn NIST (2012): $$\frac{V_s}{f_0 h} \le 20 \quad \text{hoặc} \quad V_s \le 600\text{ m/s}$$ Thuật toán giải bài toán tổng thể kết hợp phương trình vi phân chuyển động đa bậc tự do: $$\mathbf{M} \ddot{\mathbf{u}}(t) + \mathbf{C} \dot{\mathbf{u}}(t) + \mathbf{K}_{macro}(t) \mathbf{u}(t) = -\mathbf{M} \mathbf{r} \ddot{u}_g(t)$$ được số hóa và giải bằng phương pháp tích phân từng bước Newmark với thuật toán lặp cân bằng phi tuyến Newton-Raphson trên MATLAB.

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ủ chặt chẽ triết lý thực chứng (positivism) và cơ học duy thực (critical realism), kết hợp phương pháp nghiên cứu liên ngành giữa cơ học lý thuyết, mô phỏng số và thực nghiệm vật lý chính xác.

Thiết kế nghiên cứu được tổ chức theo cấu trúc đa tầng (multi-level design):

  • Mức độ cục bộ (Local scale): Khảo sát ứng xử cơ học của vật liệu cát hạt rời (xác định các thông số góc ma sát trong $\phi = 32.2^\circ$, khối lượng thể tích $\gamma = 16.5\text{ kN/m}^3$, vận tốc sóng cắt $V_s$).
  • Mức độ trung gian (Meso scale): Thiết lập phần tử vĩ mô tương đương tại vị trí tiếp xúc móng - đất nền.
  • Mức độ tổng thể (Global scale): Khảo sát toàn bộ phản ứng động học của hệ kết cấu cầu/nhà nhiều tầng đặt trên móng nông dưới tác động của các trận động đất lịch sử.

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm bàn rung được thiết kế và thực hiện tại phòng thí nghiệm chuyên sâu với các tiêu chuẩn khắt khe:

  1. Thiết bị và Hộp đất: Sử dụng hệ thống bàn rung điện động R202(UTC) có khả năng tạo rung chấn điều hòa và mô phỏng các giản đồ gia tốc động đất thực tế. Hộp đất thí nghiệm dạng tường cứng kích thước tối ưu $1.2\text{m} \times 0.6\text{m} \times 0.8\text{m}$ (chiều dài $\times$ rộng $\times$ cao), được gia cường chống biến dạng biên và bọc lớp đệm hấp thụ sóng phản xạ ở hai đầu.
  2. Quy luật tương tự mô hình (Scaling Laws): Thiết kế mô hình thu nhỏ dựa trên lý thuyết phân tích thứ nguyên của Meymand (1998) và Iai (1989) trong trường trọng lực tự nhiên $1g$, đảm bảo sự tương thích về hệ số tỷ lệ hình học $\lambda$, tỷ lệ ứng suất, tần số và gia tốc.
  3. Thiết bị đo đạc và Triangulation: Sử dụng hệ thống cảm biến gia tốc áp điện độ nhạy cao (Accelerometers) gắn tại đế bàn rung, đáy móng, đỉnh móng và đỉnh khối lượng kết cấu; cảm biến đo chuyển vị tuyến tính biến trở LVDT để ghi nhận chuyển vị ngang tức thời; kết nối đồng bộ qua bộ thu nhận và xử lý tín hiệu đa kênh số.
┌─────────────────────────────────────────────────────────────────────────────┐
│           CHUỖI KỊCH BẢN THÍ NGHIỆM BÀN RUNG R202(UTC) VÀ ĐỐI CHỨNG         │
├──────────────────────────────────────┬──────────────────────────────────────┤
│  MÔ HÌNH HỆ MÓNG - ĐẤT NỀN           │  MÔ HÌNH HỆ KẾT CẤU - MÓNG - ĐẤT     │
│  (Không có kết cấu phần trên)        │  (Có trụ và khối lượng tập trung)    │
├──────────────────────────────────────┼──────────────────────────────────────┤
│ • T12: Kích thích tần số thấp        │ • T21: Móng đặt trên mặt đất (0% D)  │
│ • T13: Kích thích cộng hưởng móng    │ • T22: Móng chôn sâu 50% chiều cao   │
│ • T14: Biên độ gia tốc trung bình    │ • T23: Móng chôn sâu 100% chiều cao  │
│ • T15: Biên độ gia tốc cực đại       │ • T24: Gia tải động đất El-Centro    │
│                                      │ • T25: Gia tải động đất Tolmezzo     │
│                                      │ • T26: Khảo sát biến thiên độ cứng K1│
└──────────────────────────────────────┴──────────────────────────────────────┘

Data và phân tích

Dữ liệu thu nhận từ thực nghiệm được xử lý qua bộ lọc thông dải (band-pass filter) để loại bỏ nhiễu tần số cao của môi trường trước khi so sánh với kết quả mô phỏng từ phần tử vĩ mô:

  • Mô hình cát thí nghiệm được hiệu chuẩn độc lập bằng phần mềm Cyclic1D để xác định chính xác đường cong suy giảm mô đun trượt $G/G_{max}$ và tỷ số cản $\zeta$.
  • Hệ phương trình vi phân chuyển động được lập trình giải tự động trên MATLAB, tích hợp thuật toán ma trận độ cứng phần tử vĩ mô tiếp tuyến cập nhật sau mỗi bước thời gian $\Delta t = 0.005\text{ s}$.
  • Đánh giá kiểm chứng độ tin cậy được thực hiện thông qua đối sánh trực tiếp ba kênh dữ liệu: (1) Kết quả thực nghiệm bàn rung, (2) Kết quả phân tích bằng mô hình phần tử vĩ mô đề xuất, (3) Kết quả phân tích bằng mô hình Paolucci (2008) và phương pháp phần tử hữu hạn chi tiết.

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

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

Kết quả nghiên cứu của luận án cung cấp những phát hiện định lượng mang tính đột phá về mặt cơ học công trình:

┌─────────────────────────────────────────────────────────────────────────────┐
│          BẢNG SO SÁNH HIỆU QUẢ CƠ HỌC VÀ SAI SỐ TÍNH TOÁN CỦA CÁC MÔ HÌNH   │
├──────────────────────┬────────────────────┬──────────────────┬──────────────┤
│ Thông số đánh giá    │ Mô hình ngàm cứng  │ Mô hình Macro    │ Thực nghiệm  │
│                      │ (Truyền thống)     │ (Luận án đề xuất)│ bàn rung     │
├──────────────────────┼────────────────────┼──────────────────┼──────────────┤
│ Chu kỳ dao động T    │ Bị đánh giá thấp   │ Tăng 18% - 35%   │ Trùng khớp   │
│ Gia tốc đỉnh kết cấu │ Quá cao (bảo thủ)  │ Giảm 15% - 28%   │ Tương đương  │
│ Chuyển vị đỉnh d_max │ Quá thấp (nguy hiểm)│ Tăng 20% - 45%  │ Sai số < 8%  │
│ Thời gian tính toán  │ Nhanh (nhưng sai)  │ Cực nhanh (giây) │ Thời gian thực│
│ Xét hiện tượng Uplift│ KHÔNG thể          │ RẤT CHÍNH XÁC    │ Quan sát rõ  │
└──────────────────────┴────────────────────┴──────────────────┴──────────────┤
  1. Hiệu ứng lọc gia tốc nhờ nhổ móng (Uplift Isolation Effect): Dưới tác dụng của động đất biên độ lớn, hiện tượng nhổ móng và biến dạng dẻo đất nền đóng vai trò như một cơ cấu cách chấn tự nhiên (natural seismic isolation). Gia tốc cực đại truyền lên đỉnh kết cấu trong mô hình phần tử vĩ mô giảm từ $15%$ đến $28%$ so với mô hình ngàm cứng quy ước.
  2. Sự gia tăng chuyển vị đỉnh kết cấu: Mặc dù gia tốc giảm, chuyển vị ngang tổng cộng tại đỉnh kết cấu lại tăng từ $20%$ đến $45%$ do sự đóng góp của góc xoay móng $\theta(t)$ và chuyển vị trượt ngang tại đáy móng. Việc bỏ qua SSI trong các tiêu chuẩn hiện hành dẫn đến việc đánh giá thấp nguy cơ va đập kết cấu liền kề (pounding) và hiệu ứng $P-\Delta$.
  3. Độ chính xác và sai số vượt trội của mô hình: Trong toàn bộ các chuỗi thí nghiệm từ T12 đến T15 (hệ móng) và T21 đến T26 (hệ kết cấu), sai số tương đối về gia tốc đỉnh và chuyển vị cực đại giữa mô hình phần tử vĩ mô đề xuất và số liệu thực nghiệm bàn rung dao động trong khoảng hẹp từ $3.2%$ đến $8.6%$. Mô hình thể hiện độ tương thích vượt trội so với mô hình cổ điển của Paolucci (1997) ở các pha dao động phi tuyến mạnh.
  4. Ảnh hưởng quyết định của độ chôn sâu móng: Thí nghiệm T26 chứng minh rằng khi tăng độ chôn sâu móng từ $0%$ lên $50%$ và $100%$ chiều cao móng, sức kháng uốn và khả năng chống trượt tăng đáng kể, hiện tượng nhổ móng bị ức chế, làm chu kỳ dao động riêng của hệ tiệm cận về trạng thái liên kết đàn hồi tuyến tính.
  5. Độ nhạy của độ cứng trụ $K_1$: Phân tích biến thiên độ cứng cho thấy khi xét tương tác SSI, mức độ suy giảm độ cứng của trụ ($1%, 5%, 10%, 50% K_1$) làm thay đổi phi tuyến bậc hai đến đáp ứng chuyển vị đỉnh, chứng minh rằng không thể tối ưu hóa kết cấu phần trên nếu tách rời độ cứng của nền móng.

Implications đa chiều

  • Về mặt lý thuyết: Luận án đã làm sáng tỏ cơ chế ghép nối (coupling mechanism) giữa tính phi tuyến hình học (tách rời tiếp xúc) và tính phi tuyến vật liệu (chảy dẻo đất nền) dưới kích động ngẫu nhiên, tạo tiền đề hoàn thiện lý thuyết động lực học công trình phi tuyến.
  • Về mặt phương pháp luận: Thiết lập thành công quy trình mô hình hóa đa mức (multi-scale modeling) kết hợp phần tử vĩ mô và mô hình thông số tập trung, mở ra giải pháp phân tích động lực học kết cấu quy mô lớn với chi phí tính toán tối thiểu.
  • Về mặt thực tiễn kỹ thuật: Cung cấp cho các kỹ sư thiết kế công trình cầu, nhà cao tầng và tháp phong điện công cụ tính toán chính xác phản ứng động đất, giúp tránh lãng phí vật liệu khi thiết kế móng quá dư thừa hoặc tránh mất an toàn khi bỏ qua chuyển vị xoay của chân công trình.
  • Về mặt chính sách và tiêu chuẩn: Cung cấp luận cứ khoa học vững chắc để các ban soạn thảo tiêu chuẩn quốc gia (TCVN 9386, TCVN 11823) chuyển đổi từ các quy định khuyến nghị định tính sang các điều khoản tính toán bắt buộc có xét đến hiệu ứng tương tác đất - kết cấu.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án cũng thẳng thắn chỉ ra các giới hạn nghiên cứu (boundary conditions):

  1. Mô hình phần tử vĩ mô trong luận án mới chỉ khảo sát cho bài toán phẳng 2D đối với móng nông đặt trên nền đất cát hạt rời đồng chất, chưa xét đến móng cọc hoặc đất dính bão hòa nước có áp lực nước lỗ rỗng tăng cao gây hóa lỏng hoàn toàn.
  2. Chưa xét đến sự ghép nối đầy đủ giữa chuyển vị ngang và góc xoay không gian 3D dưới tác động của sóng động đất đa chiều đồng thời (hai phương ngang kết hợp phương đứng).
  3. Trong mô hình thực nghiệm bàn rung, hộp đất tường cứng dù đã xử lý hấp thụ sóng biên nhưng vẫn tồn tại một phần nhỏ hiệu ứng phản xạ sóng cục bộ so với môi trường nửa không gian vô hạn thực tế.

Các định hướng nghiên cứu tiếp theo được đề xuất bao gồm:

  • Mở rộng phần tử vĩ mô sang không gian 3D xét đến tương tác tải trọng 6 bậc tự do $(N, H_x, H_y, M_x, M_y, T)$ cho các dạng móng phức tạp (móng cọc đài cao, móng giếng chìm móng cầu).
  • Tích hợp hiện tượng hóa lỏng đất cát bão hòa nước và sự suy giảm áp lực hữu hiệu vào hàm dẻo của phần tử vĩ mô.
  • Nghiên cứu thực nghiệm trên mô hình ly tâm (centrifuge test) ở gia tốc trọng trường cao ($N \times g$) để mô phỏng chính xác tuyệt đối trường ứng suất bản thân của đất nền nguyên mẫu.

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

  • Tác động học thuật: Luận án đóng góp một công trình khoa học nền tảng, được công bố trên các tạp chí chuyên ngành uy tín trong nước và quốc tế, tạo nguồn tài liệu tham khảo chất lượng cao cho cộng đồng nghiên cứu cơ học tính toán và kỹ thuật động đất tại Việt Nam.
  • Chuyển đổi công nghiệp xây dựng: Ứng dụng mô hình phần tử vĩ mô vào phần mềm phân tích kết cấu giúp các tập đoàn tư vấn thiết kế hạ tầng giao thông và công trình dân dụng tối ưu hóa kích thước móng, giảm từ $10%$ đến $20%$ khối lượng bê tông cốt thép móng mà vẫn đảm bảo khả năng chịu chấn an toàn.
  • Tác động an toàn xã hội: Giảm thiểu nguy cơ sụp đổ công trình thứ cấp do động đất tại các khu vực nhạy cảm địa chấn ở Việt Nam (vùng Tây Bắc, công trình ven biển, thủy điện bậc thang), bảo vệ tính mạng con người và tài sản quốc gia.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ kỹ thuật, Xây dựng, Cầu đường: Tiếp cận hệ thống lý thuyết chuẩn mực về phần tử vĩ mô và phương pháp tích phân phi tuyến Newmark trên MATLAB.
  • Kỹ sư Kết cấu và Địa kỹ thuật (Structural & Geotechnical Engineers): Nhận thức được mối liên kết hữu cơ giữa nền móng và kết cấu, sở hữu công cụ phân tích nhanh, chính xác để đưa ra giải pháp kháng chấn tối ưu.
  • Các Viện Nghiên cứu và Cơ quan Quản lý Tiêu chuẩn Xây dựng (Bộ Xây dựng, Bộ GTVT): Có cơ sở thực chứng định lượng để cập nhật và hoàn thiện hệ thống quy chuẩn, tiêu chuẩn thiết kế kháng chấn quốc gia.

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ì? Đó là việc tích hợp thành công mô hình mặt dẻo hóa đàn dẻo của Nova (1991) với cơ chế suy giảm diện tích tiếp xúc động lực học do nhổ móng (uplift) của Paolucci (2008) vào một phần tử vĩ mô duy nhất, giải quyết trọn vẹn cặp phi tuyến vật liệu và hình học trong miền thời gian.

2. Đột phá về phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm? So với phương pháp phần tử hữu hạn trực tiếp (Datta, 2010) tốn kém tài nguyên hay phương pháp kết cấu phụ tuyến tính (Kausel, 1978), mô hình lai của luận án giảm hơn 90% số bậc tự do nhưng vẫn bảo toàn độ chính xác vật lý cao, được kiểm chứng song song qua cả mô hình số và thực nghiệm bàn rung thực tế R202(UTC).

3. Phát hiện thực nghiệm nào gây bất ngờ nhất? Hiện tượng nhổ móng (uplift) không hoàn toàn làm mất ổn định công trình như các quan niệm thiết kế truyền thống lo ngại, mà ngược lại, hoạt động như một tầng tiêu tán năng lượng và cách chấn tự nhiên, giúp làm giảm đến 28% gia tốc truyền lên kết cấu phần trên trong các trận động đất mạnh.

4. Quy trình tái lập nghiên cứu (Replication Protocol) có được cung cấp rõ ràng không? Luận án cung cấp chi tiết toàn bộ hệ thống phương trình vi phân, sơ đồ khối thuật toán số, ma trận độ cứng giải tích tiếp tuyến, thông số địa kỹ thuật đất cát, kích thước hình học mô hình và thông số kỹ thuật bàn rung, cho phép tái lập hoàn toàn trên các phần mềm tính toán số như MATLAB, Python hoặc OpenSees.

5. Lộ trình phát triển học thuật 10 năm tiếp theo từ luận án? Lộ trình hướng tới việc phát triển phần tử vĩ mô không gian 3D đa vật liệu, tích hợp hiệu ứng hóa lỏng đất nền, xây dựng module macro-element thương mại hóa tích hợp vào các phần mềm kết cấu phổ biến (SAP2000, ETABS, MIDAS) và mở rộng cho hệ thống móng cọc công trình ngoài khơi (Offshore Wind Turbines).

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán tương tác kết cấu - đất nền (SSI) chịu tải trọng động đất bằng phương pháp lai thông qua mô hình phần tử vĩ mô, khắc phục triệt để nhược điểm của các phương pháp truyền thống.
  2. Xây dựng thành công hệ phương trình vi phân chuyển động phi tuyến kết hợp ma trận độ cứng đàn dẻo - nhổ móng và thuật toán giải số Newmark ổn định trên nền tảng MATLAB.
  3. Thực hiện thành công chuỗi thí nghiệm bàn rung vật lý thu nhỏ trên cát chặt với 10 kịch bản gia tải chi tiết, cung cấp nguồn dữ liệu thực nghiệm chuẩn xác để hiệu chuẩn mô hình lý thuyết.
  4. Chứng minh định lượng rằng việc xét tương tác SSI làm tăng chu kỳ dao động riêng của hệ từ 18% đến 35%, giảm gia tốc đỉnh kết cấu từ 15% đến 28% nhưng làm gia tăng đáng kể chuyển vị ngang đỉnh từ 20% đến 45%.
  5. Đạt độ tin cậy khoa học cao với sai số giữa tính toán lý thuyết và thực nghiệm bàn rung luôn dưới 8.6%, mở ra bước tiến mới cho công tác tính toán kháng chấn công trình tại Việt Nam.
  6. Đặt nền móng vững chắc cho việc ứng dụng phần tử vĩ mô vào thực tế kỹ thuật và cung cấp cơ sở học thuật quan trọng để nâng cấp các tiêu chuẩn thiết kế công trình chịu động đất trong tương lai.