Tổng quan luận án

Luận án tiến sĩ kỹ thuật với 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ô" được thực hiện bởi nghiên cứu sinh Huỳnh Văn Quân, chuyên ngành Cơ kỹ thuật (mã số: 9520101) tại Trường Đại học Giao thông Vận tải 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 vào năm 2021.

Tính cấp thiết và khoảng trống nghiên cứu

Trong các tiêu chuẩn thiết kế công trình hiện hành 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 và tiêu chuẩn châu Âu EN 1998 (Eurocode 8), việc phân tích phản ứng của kết cấu có xét đến tương tác với đất nền (Soil-Structure Interaction - SSI) hầu như chưa được bắt buộc tính toán chi tiết hoặc chỉ dừng lại ở mức độ khuyến nghị.

Quy trình thiết kế công trình thông thường trong thực tế được thực hiện tách biệt qua ba bước:

  1. Kỹ sư kết cấu phân tích kết cấu phần trên với giả thiết chân công trình liên kết cứng hoặc là các gối cố định trên nền đất; các giá trị nội lực chân cột (mô men, lực cắt, lực dọc) được chuyển giao cho kỹ sư địa kỹ thuật.
  2. Kỹ sư địa kỹ thuật tiến hành thí nghiệm hiện trường và trong phòng để xác định sức chịu tải và các chỉ tiêu cơ lý của nền đất.
  3. Kỹ sư nền móng sử dụng hai nguồn dữ liệu trên để thiết kế kết cấu móng đảm bảo khả năng chịu lực.

Sự chia tách này dẫn đến việc bỏ qua phản ứng tương hỗ giữa kết cấu phần trên, móng và đất nền. Dưới tác dụng của tải trọng động đất, sự rung lắc gây ra lực quán tính trong kết cấu, từ đó truyền xuống móng các thành phần lực cắt và mô men uốn. Nếu nền đất không phải là tuyệt đối cứng, các thành phần lực này sẽ gây ra chuyển vị ngang, lún đứng và góc xoay tại đáy móng, làm thay đổi đáng kể tần số dao động riêng và phân bố nội lực trong công trình.

Tại Việt Nam, các ghi nhận động đất từ Viện Vật lý địa cầu tại các khu vực thủy điện Sông Tranh, vùng Tây Bắc (Điện Biên, Sơn La, Hòa Bình) hay ngoài khơi Vũng Tàu với nguy cơ đạt cấp VIII (thang MSK-1964) đặt ra yêu cầu phải đánh giá chính xác ứng xử động của công trình. Các nghiên cứu trong nước trước đây chủ yếu dừng lại ở mô hình liên kết đàn hồi tuyến tính đơn giản. Trong khi đó, việc mô hình hóa phần tử hữu hạn chi tiết toàn bộ môi trường đất và kết cấu (phương pháp trực tiếp) đòi hỏi khối lượng tính toán lớn và tài nguyên máy tính cao. Do đó, việc ứng dụng mô hình phần tử vĩ mô (macro-element) để giải quyết bài toán tương tác kết cấu-đất nền chịu động đất là khoảng trống nghiên cứu cần được giải quyết.

Mục tiêu nghiên cứu

Mục tiêu chính của luận án là nghiên cứu ứng xử tương tác của hệ kết cấu-đất nền dưới tác dụng của tải trọng động đất bằng phương pháp lý thuyết kết hợp với thực nghiệm. Trong phân tích lý thuyết, toàn bộ tương tác phi tuyến phức tạp của hệ móng-đất nền được tích hợp vào một phần tử vĩ mô duy nhất nhằm đơn giản hóa quá trình tính toán nhưng vẫn phản ánh chính xác các đặc trưng phi tuyến của hệ.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Ứng xử động lực học của hệ kết cấu-đất nền (chuyển vị ngang và gia tốc) dưới tác dụng của tải trọng động đất.
  • Phạm vi phân tích lý thuyết: Kết cấu phần trên được mô hình hóa theo dạng thông số tập trung (Lumped-Parameter Model - LPM); trụ công trình được xem như phần tử dầm làm việc trong miền tuyến tính; phương trình vi phân chuyển động được tích phân số theo phương pháp Newmark trên môi trường Matlab. Phần tử vĩ mô được nghiên cứu áp dụng cho trường hợp móng nông đặt trên nền cát chặt và chưa xét đến cặp ứng xử chuyển vị - góc xoay.
  • Phạm vi thực nghiệm: Mô hình thí nghiệm bàn rung sử dụng hộp đất thành cứng gắn trên bàn rung R202(UTC), gia tốc kích thích tác động theo phương cạnh dài của hộp đất; bỏ qua ảnh hưởng của sự thay đổi trạng thái đất giữa các lần thí nghiệm nối tiếp nhau; các thông số đầu vào đóng vai trò là các dữ liệu cơ học kiểm chứng mô hình lý thuyết.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu về bài toán tương tác kết cấu-đất nền

Luận án đã tổng quan ba nhóm phương pháp phân tích SSI chính trong lịch sử cơ học công trình:

Phương pháp Đặc điểm và cơ sở lý thuyết Nhược điểm chính Tác giả tiêu biểu
Phương pháp trực tiếp (Direct Method) Mô hình hóa đồng thời kết cấu và toàn bộ môi trường đất bằng phần tử hữu hạn (2D hoặc 3D). Phương trình dao động tổng thể: $\mathbf{M}\mathbf{\ddot{u}} + \mathbf{C}\mathbf{\dot{u}} + \mathbf{K}\mathbf{u} = -\mathbf{M}_s\mathbf{I}\mathbf{\ddot{u}}_g$. Khối lượng tính toán rất lớn, yêu cầu cấu hình máy tính cao, phụ thuộc nhiều vào điều kiện biên nhân tạo. Datta (2010), các phần mềm thương mại FEM.
Phương pháp kết cấu phụ (Sub-structure Method) Dựa trên nguyên lý cộng tác dụng gồm 3 bước: (1) Xác định chuyển động móng cứng không khối lượng; (2) Xác định độ cứng lò xo đất nền tương đương; (3) Tính toán phản ứng kết cấu đặt trên lò xo đất. Giả thiết hệ làm việc đàn hồi tuyến tính nên không mô tả được các phi tuyến thực tế của đất và móng. Kausel (1978), John (1985), Pais (1988), Gazetas (1991), Mylonakis (2006).
Phương pháp lai / Phần tử vĩ mô (Hybrid / Macro-element Method) Thay thế toàn bộ hệ móng và vùng đất chịu ảnh hưởng bằng một phần tử đơn lẻ tích hợp các quy luật đàn dẻo, phi tuyến tiếp xúc và suy giảm độ cứng. Cần xây dựng chính xác các mặt dẻo, quy luật chảy và ma trận độ cứng suy giảm. Nova (1991), Paolucci (1997, 2008), Cremer (2001, 2002), Grange (2009), Chatzigogos (2009, 2011), Figini (2012).

Các nghiên cứu về phần tử vĩ mô trên thế giới

  • Nova (1991): Đề xuất mô hình phần tử vĩ mô đàn dẻo theo quy tắc tăng bền đẳng hướng dưới các tổ hợp lực nén đúng tâm, lệch tâm và lực cắt phẳng đối với móng băng trên nền cát.
  • Cremer (2001, 2002): Xây dựng phần tử vĩ mô xét đồng thời phi tuyến vật liệu và phi tuyến hình học (hiện tượng nhổ móng - uplift) dưới tải trọng tuần hoàn và tải trọng động đất, ứng dụng tính toán cho cầu Viaduc de l'Arc (Pháp).
  • Paolucci (1997, 2008): Ứng dụng phần tử vĩ mô vào bài toán phân tích động lực học kết cấu chịu động đất với hệ 4 bậc tự do; đề xuất công thức suy giảm bề rộng móng hiệu dụng $B' = (1-\delta)B$ phụ thuộc vào góc xoay dẻo tích lũy $\theta^p$; kiểm chứng trên dữ liệu thực nghiệm dự án PWRI.
  • Grange (2009): Mở rộng phần tử vĩ mô sang không gian 3 chiều cho móng tròn, móng chữ nhật và móng băng; kiểm chứng qua kết quả dự án Camus IV.
  • Chatzigogos (2009, 2011): Phát triển phần tử vĩ mô xét đến trượt và nhổ móng cho móng băng và móng tròn chịu tải trọng tăng dần và chu kỳ.
  • Figini (2012): Đề xuất mô hình lan truyền độ cứng dựa trên độ lệch tâm ban đầu $e_0 = \frac{B}{2\alpha} e^{-\xi Q_N}$, kiểm chứng trên kết quả các dự án Camus IV, TRISEE-1 và TRISEE-2.

Tình hình nghiên cứu trong nước và vị trí của luận án

Tại Việt Nam, các nghiên cứu của Lê Văn Tuân (2016), Đào Văn Hưng (2017), Vũ Minh Ngọc (2019) và Trần Thu Hằng (2019) đã khảo sát bài toán tương tác đất - kết cấu công trình ngầm hoặc móng nhưng phần lớn sử dụng mô hình lò xo đàn hồi tuyến tính hoặc mô hình phần tử hữu hạn thông thường. Luận án của Huỳnh Văn Quân là công trình trong nước đầu tiên đề xuất và sử dụng thuật ngữ tiếng Việt "phần tử vĩ mô" cho khái niệm macro-element, đồng thời thiết lập hệ thống phương trình vi phân chuyển động hoàn chỉnh tích hợp phần tử vĩ mô và kiểm chứng bằng thực nghiệm bàn rung chuyên biệt.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết tương tác kết cấu - đất nền

Theo lý thuyết SSI, chuyển động của hệ kết cấu - móng - đất nền dưới tải trọng động đất gồm hai thành phần:

  1. Tương tác động (Kinematic Interaction): Chuyển động của móng bị sai khác so với chuyển động của trường đất tự do do sự hiện diện của móng cứng và hình học chôn sâu, gây ra các hiện tượng phản xạ, khúc xạ sóng và chuyển vị xoay của móng.
  2. Tương tác quán tính (Inertial Interaction): Khối lượng của kết cấu phần trên bị kích thích dao động, sinh ra lực quán tính, mô men uốn và lực cắt tác động ngược trở lại nền đất qua mặt tiếp xúc đáy móng.

Điều kiện cần thiết để xét đến tương tác SSI trong tính toán kỹ thuật là: $$\frac{f \cdot h}{V_s} \ge 0{,}1$$ Trong đó $f$ là tần số dao động cơ bản của kết cấu ngàm cứng, $h$ là chiều cao công trình, $V_s$ là vận tốc truyền sóng cắt trong đất. Với công trình nhà nhiều tầng ($h/n = 3\text{ m}$), điều kiện tương đương là $V_s \le 600\text{ m/s}$. Với kết cấu dạng công-xon (trụ cầu, tháp), điều kiện là $V_s \le 57580/\lambda_s$ (kết cấu thép) và $V_s \le 123970/\lambda_s$ (kết cấu bê tông) với $\lambda_s = h/r$ là độ mảnh.

                   [ Kết cấu phần trên (LPM) ]
                                |
                   [ Trụ dầm tuyến tính (EI) ]
                                |
           +--------------------+--------------------+
           |             PHẦN TỬ VĨ MÔ               |
           |  (Thay thế toàn bộ hệ móng - đất nền)   |
           +--------------------+--------------------+
                                |
             +------------------+------------------+
             |                                     |
    [ VÙNG GẦN (Near Field) ]            [ VÙNG XA (Far Field) ]
    - Phi tuyến hình học (nhổ móng)       - Đàn hồi tuyến tính
    - Phi tuyến vật liệu (chảy dẻo)       - Cản bức xạ sóng
    - Hệ lò xo phi tuyến                 - Lò xo & Giảm chấn nhớt

Cơ sở xây dựng phần tử vĩ mô

Phần tử vĩ mô thay thế toàn bộ hệ móng - đất nền bằng cách phân chia trường tương tác thành hai vùng:

  • Vùng gần (Near Field): Khu vực đất lân cận đáy móng chịu biến dạng lớn, xuất hiện hai dạng phi tuyến:
    • Phi tuyến hình học: Hiện tượng móng bị nhổ lên khỏi mặt đất (uplift) khi mô men vượt quá giới hạn cân bằng tĩnh $M_0 = \pm \frac{N B}{4}$, làm suy giảm diện tích tiếp xúc thực tế.
    • Phi tuyến vật liệu: Sự dẻo hóa của đất nền khi ứng suất chạm tới mặt chảy dẻo cực hạn, hiện tượng tập trung ứng suất tại mép móng.
  • Vùng xa (Far Field): Khu vực đất nằm xa móng, làm việc hoàn toàn trong miền đàn hồi tuyến tính và tiêu tán năng lượng dao động thông qua lực cản nhớt (cản bức xạ sóng).

Phương pháp nghiên cứu thực nghiệm và phân tích số

  1. Phương pháp phân tích số: Thiết lập ma trận khối lượng $\mathbf{M}$, ma trận cản $\mathbf{C}$ và ma trận độ cứng tiếp tuyến $\mathbf{K}^{ep}$ của hệ kết cấu - phần tử vĩ mô. Sử dụng thuật toán tích phân số Newmark từng bước thời gian kết hợp lặp phi tuyến trên phần mềm tính toán chuyên dụng Matlab (chương trình CyclicTP).
  2. Phương pháp thực nghiệm: Chế tạo mô hình thu nhỏ dựa trên lý thuyết tương tự hình học và động lực học của Meymand (1998) với hệ số tỷ lệ $\lambda$. Thí nghiệm được thực hiện trên hệ thống bàn rung R202(UTC) với các băng gia tốc động đất nhân tạo và động đất thực tế (El-Centro 1940, Tolmezzo).

Nội dung chính theo từng chương

Chương 1: Tổng quan về tương tác kết cấu-đất nền

Chương 1 hệ thống hóa cơ sở lý thuyết về bài toán tương tác động và tương tác quán tính giữa công trình và nền đất. Tác giả phân tích các dạng hư hỏng thực tế của công trình do không xét đến hiện tượng xoay và lún lệch của móng (như các vết nứt khung bê tông do chuyển vị góc xoay của tường chịu cắt).

Chương này trình bày chi tiết các phương pháp phân tích SSI cổ điển và hiện đại, đồng thời phân tích sâu bản chất toán học của các mô hình phần tử vĩ mô điển hình:

  • Mô hình tăng bền đẳng hướng của Nova (1991).
  • Mô hình móng băng xét nhổ móng của Cremer (2001, 2002) với điều kiện lật móng hoàn toàn tại $M_C = \pm \frac{N B}{2}$.
  • Mô hình 4 bậc tự do của Paolucci (1997, 2008) với hàm dẻo $f(\mathbf{F})$ và quy luật chảy dẻo $g(\mathbf{F})$.
  • Các mô hình không gian của Grange (2009), Chatzigogos (2009, 2011) và mô hình suy giảm độ cứng của Figini (2012).

Phần cuối chương tổng kết các công trình thực nghiệm bàn rung và thí nghiệm ly tâm trên thế giới đối với móng đơn, móng giếng chìm, móng xa bờ (Maugeri 2000, Tsukamoto 2012, Ibsen 2015, Barari 2017, Zeybek 2017) và các nghiên cứu tại Việt Nam để xác định phạm vi nghiên cứu của đề tài.

Chương 2: Mô hình tương tác kết cấu-đất nền chịu tải trọng động đất bằng phần tử vĩ mô

Nội dung chương 2 tập trung xây dựng cơ sở toán học và cơ học cho phần tử vĩ mô phẳng:

  1. Đặc trưng lực và chuyển vị thu gọn: Thiết lập véc-tơ lực thu gọn không thứ nguyên $\mathbf{F} = [V/N_{max}, ; M/(B N_{max}), ; N/N_{max}]^T$ và véc-tơ chuyển vị thu gọn tương ứng $\mathbf{u} = [u/B, ; \theta, ; v/B]^T$. Khả năng chịu tải cực hạn thẳng đứng $N_{max} = B q_{max}$ được xác định theo lý thuyết sức chịu tải của Terzaghi có xét đến chiều sâu chôn móng $D_f$.
  2. Thiết lập hàm dẻo và quy luật chảy: Xây dựng mặt dẻo giới hạn $f(\mathbf{F}) = 0$ và tiềm năng dẻo $g(\mathbf{F})$ mô tả trạng thái ứng suất giới hạn của nền đất dưới đáy móng.
  3. Thiết lập ma trận độ cứng đàn dẻo: Ma trận độ cứng tiếp tuyến $\mathbf{K}^{ep}$ của phần tử vĩ mô được xác định qua biểu thức: $$\mathbf{K}^{ep} = \mathbf{K}^F - \frac{\mathbf{K}^F \left(\frac{\partial g}{\partial \mathbf{F}}\right) \left(\frac{\partial f}{\partial \mathbf{F}}\right)^T \mathbf{K}^F}{\left(\frac{\partial f}{\partial \mathbf{F}}\right)^T \mathbf{K}^F \left(\frac{\partial g}{\partial \mathbf{F}}\right)}$$
  4. Hệ phương trình vi phân chuyển động: Thành lập phương trình dao động cho hai mô hình:
    • Mô hình hệ móng - đất nền (khảo sát riêng phản ứng của móng).
    • Mô hình hệ kết cấu phần trên - móng - đất nền (kết cấu dạng khối lượng tập trung trên trụ dầm đàn hồi liên kết với phần tử vĩ mô).
  5. Thuật toán giải số: Lập trình giải hệ phương trình vi phân phi tuyến bằng thuật toán tích phân Newmark kết hợp vòng lặp kiểm tra trạng thái mặt dẻo tại từng bước thời gian $\Delta t$ trên phần mềm CyclicTP/Matlab dưới tác động của gia tốc sóng động đất El-Centro (1940).

Chương 3: Nghiên cứu thực nghiệm tương tác hệ kết cấu-đất nền dưới tải trọng động đất

Chương 3 trình bày toàn bộ quy trình thiết kế, chế tạo và thực hiện thí nghiệm mô hình trên bàn rung R202(UTC):

  • Nguyên lý tỷ lệ mô hình: Dựa trên các hệ số tỷ lệ của Meymand (1998), tác giả tính toán kích thước mô hình, khối lượng bản đỉnh, chiều cao trụ và kích thước hộp đất chứa cát thí nghiệm.
  • Vật liệu và thiết bị: Nền đất là cát hạt vừa được xác định các chỉ tiêu cơ lý qua mô hình Cyclic1D; móng và kết cấu mô hình được chế tạo bằng vật liệu kim loại/bê tông chuẩn hóa. Hệ thống đo lường bao gồm các cảm biến gia tốc và cảm biến đo chuyển vị điện tử gắn tại các vị trí: bàn rung, đáy móng, đỉnh móng và đỉnh kết cấu phần trên.
  • Kịch bản thí nghiệm: Sử dụng kích thích sóng động đất Tolmezzo và gia tốc điều hòa với các chuỗi thử nghiệm:
    • Chuỗi thí nghiệm hệ móng - đất nền (ký hiệu các trường hợp tải trọng T12, T13, T14, T15).
    • Chuỗi thí nghiệm hệ kết cấu phần trên - móng - đất nền (ký hiệu các trường hợp T21, T22, T23, T24, T25, T26).
    • Khảo sát các điều kiện chôn móng khác nhau: móng đặt trực tiếp trên mặt đất, móng chôn sâu 50% chiều cao và móng chôn sâu 100% chiều cao.
  • Kết quả thực nghiệm: Ghi nhận chuỗi thời gian gia tốc và chuyển vị đỉnh móng, đỉnh kết cấu phần trên; phân tích hiện tượng suy giảm tần số dao động và sự trượt/tách của móng trong quá trình rung lắc mạnh.

Chương 4: Phân tích ứng xử kết cấu chịu động đất bằng phần tử vĩ mô

Chương 4 thực hiện so sánh, đối chứng kết quả giữa tính toán lý thuyết từ mô hình phần tử vĩ mô đề xuất với kết quả đo đạc thực nghiệm bàn rung:

Nhóm mô hình kiểm tra Trường hợp khảo sát Đại lượng so sánh Đánh giá độ phù hợp và sai số
Hệ móng - đất nền Các kịch bản gia tải từ T12-00 đến T15-00 Lịch sử thời gian gia tốc đỉnh móng, gia tốc cực đại đỉnh móng Dạng sóng lý thuyết bám sát diễn biến thực nghiệm; giá trị đỉnh gia tốc có sự tương đồng cao với sai số nhỏ.
Hệ kết cấu phần trên - móng - đất nền Các kịch bản từ T21-00 đến T26-00 Lịch sử chuyển vị và gia tốc đỉnh kết cấu phần trên Mô hình phản ánh được biên độ chuyển vị ngang cực đại và sự lệch pha dao động do biến dạng dẻo nền đất gây ra.
So sánh mô hình đối chứng Mô hình đề xuất vs Mô hình Paolucci (2008) vs Thực nghiệm Gia tốc và chuyển vị đỉnh trong kịch bản T21-00 và T24-00 Mô hình đề xuất cho kết quả tiệm cận với số liệu thực nghiệm hơn ở pha dao động sau cực đại.

Chương này cũng tiến hành phân tích độ nhạy của độ cứng $K_S$ (hay các mức độ cứng suy giảm $k_1$ lấy bằng 1%, 5%, 10% và 50% độ cứng ban đầu) đến gia tốc và chuyển vị ngang của kết cấu phần trên, chỉ ra quy luật biến thiên phản ứng công trình khi đất nền bị suy giảm độ cứng trong động đất.

Kết quả và những đóng góp mới

Đóng góp về mặt lý luận và khoa học

  1. Xây dựng hoàn chỉnh mô hình phần tử vĩ mô phẳng thích hợp cho bài toán phân tích động lực học công trình chịu tải trọng động đất, tích hợp đồng thời cơ chế phi tuyến vật liệu của đất nền và quy luật phi tuyến hình học (nhổ móng).
  2. Thiết lập quy trình toán học và thuật toán số giải hệ phương trình vi phân dao động của hệ kết cấu - móng - đất nền dựa trên phương pháp tích phân Newmark, cho phép giảm đáng kể số bậc tự do so với mô hình phần tử hữu hạn truyền thống nhưng vẫn bảo toàn các đặc tính cơ học cơ bản.
  3. Chuẩn hóa thuật ngữ khoa học chuyên ngành trong nước bằng việc đề xuất tên gọi "phần tử vĩ mô" cho khái niệm quốc tế macro-element.

Đóng góp về mặt thực tiễn

  1. Cung cấp bộ cơ sở dữ liệu thực nghiệm bàn rung tin cậy về ứng xử động lực học (gia tốc, chuyển vị) của hệ móng đơn và hệ kết cấu - móng trên nền cát dưới tác dụng của các kích thích động đất chuẩn.
  2. Cung cấp công cụ tính toán số (chương trình CyclicTP trên môi trường Matlab) giúp các kỹ sư kết cấu và địa kỹ thuật có thể đánh giá nhanh tương tác SSI trong các dự án công trình cầu, tháp hoặc nhà nhiều tầng mà không bắt buộc phải sử dụng các phần mềm phần tử hữu hạn 3D phức tạp.
  3. Đưa ra các phân tích định lượng về ảnh hưởng của chiều sâu chôn móng (0%, 50%, 100%) và sự suy giảm độ cứng tiếp xúc $K_S$ đến phản ứng của công trình, phục vụ công tác thiết kế kháng chấn an toàn và kinh tế.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

  • Mô hình phần tử vĩ mô mới được thiết lập cho bài toán phẳng 2D đối với móng nông dạng chữ nhật/băng đặt trên nền cát đồng nhất, chưa xét đầy đủ tương tác không gian 3 chiều và cặp ứng xử chuyển vị - góc xoay đồng thời.
  • Kết cấu phần trên mới được mô hình hóa ở dạng 1 bậc tự do tập trung trên trụ dầm làm việc trong giới hạn đàn hồi tuyến tính, chưa xét đến tính phi tuyến vật liệu (chảy dẻo bê tông, hình thành khớp dẻo) của thân công trình.
  • Thí nghiệm bàn rung sử dụng mô hình hộp đất thành cứng nên không thể triệt tiêu hoàn toàn hiện tượng phản xạ sóng từ thành hộp; các chuỗi thí nghiệm kế tiếp chưa xét đến sự biến đổi tích lũy trạng thái độ chặt của cát sau các lần rung trước.

Hướng phát triển tiếp theo

  • Mở rộng xây dựng phần tử vĩ mô không gian 3D cho các loại hình móng khác (móng cọc, móng bè, móng giếng chìm) trên các điều kiện địa chất phân tầng phức tạp hoặc đất dính.
  • Tích hợp tính phi tuyến đàn dẻo của vật liệu kết cấu phần trên (kết hợp phi tuyến kết cấu và phi tuyến đất nền).
  • Thực hiện các thí nghiệm ly tâm (Centrifuge Test) để tái hiện chính xác hơn trường ứng suất trọng trường của nền đất nguyên mẫu.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các đối tượng và lĩnh vực sau:

  • Nghiên cứu sinh và học viên cao học: Tham khảo phương pháp xây dựng ma trận độ cứng tiếp tuyến đàn dẻo, phương pháp lập trình tích phân Newmark giải bài toán phi tuyến động lực học trong Matlab.
  • Giảng viên và nhà nghiên cứu Cơ kỹ thuật, Công trình Giao thông: Sử dụng làm tài liệu giảng dạy chuyên đề Động lực học công trình nâng cao, Tương tác kết cấu - môi trường và Cơ học đất động.
  • Kỹ sư thiết kế kết cấu và địa kỹ thuật: Tham khảo quy trình đánh giá ảnh hưởng của tương tác SSI, điều kiện bắt buộc phải xét SSI theo vận tốc sóng cắt $V_s$, và ứng dụng công cụ đơn giản hóa để kiểm tra phản ứng của công trình cầu, tháp chịu động đất.

Câu hỏi thường gặp

1. Thuật ngữ "phần tử vĩ mô" (Macro-element) trong luận án được định nghĩa như thế nào?

"Phần tử vĩ mô" là phần tử đơn lẻ đại diện cho toàn bộ hệ móng - đất nền đặt tại chân kết cấu. Phần tử này tích hợp các ứng xử cơ học của hai vùng: vùng gần (gồm phi tuyến hình học như nhổ móng và phi tuyến vật liệu như chảy dẻo) và vùng xa (đàn hồi tuyến tính và cản nhớt lan truyền sóng), giúp mô hình hóa toàn bộ nền móng trong phân tích động lực học mà không cần chia lưới phần tử hữu hạn chi tiết cho đất.

2. Khi nào kỹ sư bắt buộc phải xét đến bài toán tương tác kết cấu - đất nền (SSI)?

Theo luận án trích dẫn từ các nghiên cứu cơ sở, tương tác SSI bắt buộc phải xét đến khi thỏa mãn điều kiện: $$\frac{f \cdot h}{V_s} \ge 0{,}1$$ Đối với công trình dân dụng nhiều tầng có tỷ số tầng $h/n = 3\text{ m}$, điều kiện này tương ứng với nền đất có vận tốc sóng cắt $V_s \le 600\text{ m/s}$. Đối với công trình dạng công-xon, điều kiện là $V_s \le 57580/\lambda_s$ (kết cấu thép) và $V_s \le 123970/\lambda_s$ (kết cấu bê tông) với $\lambda_s$ là độ mảnh.

3. Phi tuyến hình học của hệ móng - đất nền dưới tải trọng động đất khác gì so với tải trọng tĩnh?

Dưới tải trọng động đất, phi tuyến hình học (hiện tượng tách rời và nhổ móng - uplift) bị chi phối bởi tương tác quán tính. Chiều của lực quán tính và mô men quay liên tục đảo chiều theo sự đổi dấu của gia tốc động đất. Do đó, phần diện tích tiếp xúc giữa đáy móng và nền đất liên tục biến thiên theo thời gian: có thể bị thu hẹp lại khi móng xoay và phục hồi trở lại trạng thái tiếp xúc ban đầu khi hệ lắc theo chiều ngược lại (miễn là công trình chưa bị lật hoàn toàn).

4. Thuật toán số nào được sử dụng để giải hệ phương trình vi phân chuyển động trong luận án?

Luận án sử dụng phương pháp tích phân số Newmark trên miền thời gian, được lập trình thông qua công cụ CyclicTP trên môi trường phần mềm Matlab để giải hệ phương trình vi phân dao động từng bước thời gian $\Delta t$, đồng thời cập nhật ma trận độ cứng tiếp tuyến đàn dẻo $\mathbf{K}^{ep}$ của phần tử vĩ mô tại mỗi bước lặp.

5. Thí nghiệm bàn rung trong luận án sử dụng thiết bị và các kích thích động lực học nào?

Thí nghiệm được thực hiện trên hệ thống bàn rung R202(UTC) tại Phân hiệu Trường Đại học Giao thông Vận tải. Kích thích động lực học đưa vào bao gồm các băng gia tốc động đất thực tế đã chuẩn hóa như sóng động đất El-Centro (1940) và sóng động đất Tolmezzo nhằm khảo sát phản ứng gia tốc và chuyển vị của móng và kết cấu.

Kết luận

Luận án tiến sĩ của Huỳnh Văn Quân đã giải quyết bài toán phân tích tương tác kết cấu - đất nền chịu tải trọng động đất bằng cách tiếp cận phương pháp lai sử dụng phần tử vĩ mô. Nghiên cứu đã thành lập thành công hệ phương trình động lực học xét đồng thời phi tuyến vật liệu của đất và phi tuyến hình học nhổ móng, giải số bằng thuật toán Newmark trên nền tảng Matlab. Các kết quả phân tích lý thuyết đã được kiểm chứng bằng chuỗi thí nghiệm bàn rung R202(UTC) với độ tin cậy cao, mở ra hướng ứng dụng tính toán kháng chấn hiệu quả cho các công trình xây dựng và giao thông.