Tổng quan luận án

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, tiêu chuẩn thiết kế cầu 22TCN 272-05, tiêu chuẩn AASHTO LRFD của Hoa Kỳ, tiêu chuẩn thiết kế công trình chịu động đất TCVN 9386:2012 hay tiêu chuẩn châu Âu EN 1998 (Eurocode 8), việc phân tích ứng xử của công trình có xét đến tương tác kết cấu - đất nền (Soil-Structure Interaction - SSI) hầu như chưa được bắt buộc hoặc chỉ dừng lại ở mức độ khuyến nghị.

Thực tế thiết kế thường chia tách thành hai bước độc lập: kỹ sư kết cấu phân tích công trình phần trên với giả thiết chân cột/móng ngàm cứng hoặc đặt trên các gối cố định để xuất nội lực chân công trình; sau đó kỹ sư địa kỹ thuật sử dụng các giá trị nội lực này để tính toán móng trên nền đất. Quy trình này bỏ qua tính mềm của đất nền và sự tương tác qua lại giữa kết cấu và môi trường địa kỹ thuật. Khi chịu tải trọng động đất, các phản ứng quán tính và chuyển động sóng trong đất gây ra các hiện tượng phi tuyến hình học (như hiện tượng tách rời đáy móng khỏi mặt đất - uplift) và phi tuyến vật liệu (trạng thái dẻo của đất, tập trung ứng suất).

Tại Việt Nam, các ghi nhận địa chấn từ Viện Vật lý địa cầu cho thấy nguy cơ động đất tại nhiều khu vực (thủy điện Sông Tranh, các tỉnh Điện Biên, Hòa Bình, Sơn La và vùng biển ngoài khơi Vũng Tàu), với khả năng xảy ra động đất cấp VIII theo thang MSK-1964 ở miền Bắc. Phần lớn các nghiên cứu trong nước trước đây chỉ áp dụng mô hình liên kết đàn hồi tuyến tính (hệ lò xo Winkler độc lập), chưa phản ánh đúng bản chất phi tuyến phức tạp khi chịu tải trọng động lực học mạnh. Phương pháp phần tử vĩ mô (macro-element) là hướng tiếp cận dung hòa giữa độ chính xác của phương pháp phần tử hữu hạn phi tuyến trực tiếp và tính tinh gọn của phương pháp tham số tập trung, nhưng chưa được hệ thống hóa và kiểm chứng thực nghiệm đầy đủ trong điều kiện nghiên cứu tại Việt Nam.

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

Luận án xác định các mục tiêu cụ thể:

  1. Nghiên cứu cơ sở lý thuyết về 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, xây dựng mô hình phần tử vĩ mô đại diện cho hệ móng - đất nền có xét đến đồng thời phi tuyến hình học và phi tuyến vật liệu.
  2. Thiết lập hệ phương trình vi phân chuyển động cho hệ móng - đất nền và hệ kết cấu phần trên - móng - đất nền dưới dạng mô hình thông số tập trung (Lumped-Parameter Model - LPM), giải bằng phương pháp tích phân số Newmark.
  3. Thiết kế và thực hiện chuỗi thí nghiệm bàn rung (shaking table test) trên mô hình thu nhỏ đối với hệ móng - đất nền và hệ kết cấu phần trên - móng - đất nền chịu kích thích động đất để thu thập dữ liệu ứng xử thực tế.
  4. Kiểm chứng mô hình lý thuyết đề xuất thông qua so sánh với kết quả thực nghiệm bàn rung và phân tích ảnh hưởng của độ cứng đất nền đến phản ứng động học của kết cấu phần trên.

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

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

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

Các hướng tiếp cận nghiên cứu trên thế giới và trong nước

Luận án hệ thống hóa các hướng nghiên cứu tương tác kết cấu - đất nền qua ba nhóm phương pháp phân tích chính:

  1. Phương pháp trực tiếp (Direct Method): Toàn bộ nền đất, móng và kết cấu được chia lưới phần tử hữu hạn (Datta, 2010). Phương pháp này mô tả chi tiết hình học và vật liệu nhưng đòi hỏi khối lượng tính toán lớn, cấu hình máy tính cao và nhạy cảm với điều kiện biên hấp thụ sóng.
  2. Phương pháp kết cấu phụ (Substructure Method): Dựa trên nguyên lý cộng tác dụng do Kausel (1978) đề xuất, phân chia bài toán thành: xác định chuyển động động học móng không khối lượng, xác định độ cứng lò xo đất (Pais, 1988; Gazetas, 1991; Mylonakis, 2006) và giải phản ứng kết cấu. Nhược điểm cốt lõi là chỉ áp dụng cho hệ làm việc tuyến tính.
  3. Phương pháp lai và Phần tử vĩ mô (Macro-element):
    • Nova (1991): Đề xuất phần tử đàn dẻo đẳng hướng dưới tải trọng nén lệch tâm và nghiêng, xây dựng hàm mặt chảy $f(F)$ và quy luật chảy $g(F)$ chuẩn hóa theo tải trọng giới hạn thẳng đứng.
    • Cremer (2001, 2002): Phát triển phần tử vĩ mô cho móng băng chịu tải trọng chu kỳ có xét đến hiện tượng nhổ móng (uplift) và trạng thái dẻo của đất, áp dụng phân tích cầu Viaduc de l'Arc (Pháp).
    • Paolucci (1997, 2008): Ứng dụng phần tử vĩ mô 3 bậc tự do vào bài toán động đất kết hợp với mô hình 4 bậc tự do của kết cấu, phát triển hàm suy giảm bề rộng móng hiệu dụng $\delta(\theta^p)$ để mô phỏng phi tuyến hình học từ dữ liệu dự án PWRI.
    • Grange (2009): Mở rộng phần tử vĩ mô sang không gian 3D cho móng tròn, chữ nhật và móng băng, kiểm chứng với số liệu dự án Camus IV.
    • Chatzigogos (2009, 2011): Xây dựng phương trình liên tục cho phần tử vĩ mô xét cơ chế trượt tại mặt tiếp xúc móng - đất.
    • Figini (2012): Đề xuất mô hình suy giảm độ cứng dựa trên độ lệch tâm ban đầu $e_0$, kiểm nghiệm trên dữ liệu Camus IV, TRISEE-1 và TRISEE-2.

Về phương diện thực nghiệm, các công trình quốc tế sử dụng bàn rung và máy quay ly tâm tiêu biểu gồm: Maugeri (2000) tại Đại học Bristol; Tsukamoto (2012) khảo sát độ lún và áp lực nước lỗ rỗng; Ibsen (2015) nghiên cứu móng giếng chìm; Barari (2017) nghiên cứu móng trụ đơn xa bờ tại Đại học Aalborg; Adamidis (2018) và Zeybek (2017, 2020) nghiên cứu hóa lỏng trên mô hình ly tâm; Anastasopoulos (2012) nghiên cứu trụ đơn và hệ khung.

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) chủ yếu sử dụng mô hình lò xo tuyến tính hoặc thí nghiệm mô hình công trình ngầm trong hộp cát.

Khoảng trống nghiên cứu luận án giải quyết

Luận án là công trình đầu tiên tại Việt Nam tiếp cận bài toán SSI dưới tải trọng động đất bằng phần tử vĩ mô (đề xuất thuật ngữ tiếng Việt chính thức cho khái niệm macro-element). Luận án khắc phục các hạn chế của việc phân tích phần tử hữu hạn truyền thống (vốn phức tạp, tốn kém tài nguyên tính toán) và phương pháp giải tích tuyến tính (không phản ánh hiện tượng nhổ móng và biến dạng dẻo), đồng thời xây dựng một chuỗi thí nghiệm bàn rung độc lập để kiểm chứng trực tiếp độ tin cậy của thuật toán đề xuất.

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

Cơ sở lý thuyết phần tử vĩ mô

Phần tử vĩ mô thay thế toàn bộ hệ thống móng và nền đất xung quanh bằng một nút phần tử tập trung đặt tại cao trình đáy móng, phân tách miền đất thành hai vùng:

  • Vùng gần (near field): Khu vực đất sát đáy móng, nơi xuất hiện các biến dạng phi tuyến hình học (tách đáy móng/uplift) và phi tuyến vật liệu (chảy dẻo). Vùng này được biểu diễn bằng hệ lò xo phi tuyến.
  • Vùng xa (far field): Khu vực đất ở xa móng, ứng xử biến dạng nằm trong giới hạn đàn hồi tuyến tính và năng lượng dao động tiêu tán qua sóng lan truyền. Vùng này được mô hình bằng hệ lò xo và giảm chấn nhớt tuyến tính mắc song song.

Hệ vector lực thu gọn $\mathbf{F} = [V \quad M/B \quad N]^T$ và vector chuyển vị thu gọn tương ứng $\mathbf{u} = [u \quad B\theta \quad v]^T$ được chuẩn hóa qua sức chịu tải thẳng đứng cực hạn của móng $N_{max} = B q_{max}$, trong đó: $$q_{max} = \frac{1}{2} \gamma B N_\gamma + q_0 N_q + c N_c$$

Trạng thái giới hạn dẻo được định nghĩa qua hàm mặt chảy $f(\mathbf{F}) \le 0$ và quy luật chảy phi liên kết $g(\mathbf{F})$. Khi $f(\mathbf{F}) < 0$, phần tử làm việc đàn hồi với ma trận độ cứng $\mathbf{K}^{el}$. Khi $f(\mathbf{F}) = 0$ và $df(\mathbf{F}) = 0$, ma trận độ cứng đàn dẻo $\mathbf{K}^{ep}$ được xác định thông qua đạo hàm của hàm chảy và thế dẻo theo vector lực thu gọn: $$\mathbf{K}^{ep} = \mathbf{K}^{el} - \frac{\mathbf{K}^{el} \frac{\partial g}{\partial \mathbf{F}} \left(\frac{\partial f}{\partial \mathbf{F}}\right)^T \mathbf{K}^{el}}{\left(\frac{\partial f}{\partial \mathbf{F}}\right)^T \mathbf{K}^{el} \frac{\partial g}{\partial \mathbf{F}}}$$

Sự suy giảm diện tích tiếp xúc móng - đất do phi tuyến hình học được theo dõi thông qua việc cập nhật bề rộng làm việc hiệu dụng $B' = B(1 - \delta)$, trong đó $\delta$ phụ thuộc vào góc xoay dẻo tích lũy $\theta^p$.

Phương trình chuyển động và phương pháp tích phân số

Phương trình vi phân chuyển động tổng quát của hệ nhiều bậc tự do: $$\mathbf{M}\ddot{\mathbf{u}} + \mathbf{C}\dot{\mathbf{u}} + \mathbf{K}\mathbf{u} = -\mathbf{M}_s \mathbf{I} \ddot{u}_g(t)$$

Trong đó:

  • $\mathbf{M}, \mathbf{C}, \mathbf{K}$: Lần lượt là ma trận khối lượng, ma trận cản nhớt (xác định theo mô hình cản Rayleigh) và ma trận độ cứng tiếp tuyến của toàn hệ (kết hợp độ cứng phần trên và độ cứng phần tử vĩ mô).
  • $\mathbf{M}_s$: Ma trận khối lượng liên kết với thành phần kích thích gia tốc nền $\ddot{u}_g(t)$.
  • $\mathbf{I}$: Vector chỉ số hướng tác dụng của gia tốc nền.

Hệ phương trình phi tuyến được giải từng bước thời gian $\Delta t$ bằng thuật toán Newmark kết hợp vòng lặp hiệu chỉnh lặp Newton-Raphson trên nền tảng phần mềm MATLAB.

Phương pháp nghiên cứu thực nghiệm

  • Nguyên lý tương tự hình học và động học: Thiết kế mô hình dựa trên định luật tỷ lệ mô hình hóa trường trọng lực đơn của Meymand (1998) với hệ số tỷ lệ hình học $\lambda$.
  • Cấu hình thí nghiệm:
    • Bàn rung: Bàn rung một phương R202(UTC) đặt tại phòng thí nghiệm.
    • Hộp đất: Hộp chứa đất dạng vách cứng, kích thước đảm bảo giảm thiểu phản xạ sóng từ thành hộp tới mô hình móng.
    • Địa chất nền: Cát khô đầm chặt, xác định thành phần cấp phối hạt và góc ma sát trong thông qua thí nghiệm trong phòng.
    • Tín hiệu động đất: Kích thích sóng địa chấn lịch sử (như sóng Tolmezzo, El-Centro 1940) được nạp vào bộ điều khiển bàn rung.
    • Hệ thống đo: Cảm biến đo gia tốc áp điện và cảm biến đo chuyển vị điện cảm LVDT kết nối qua bộ thu nhận và xử lý số liệu động.

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 tổng hợp bức tranh toàn cảnh về cơ sở lý thuyết và các nghiên cứu thực nghiệm trong lĩnh vực SSI:

  • Phân tích sự đối lập trong quan điểm thiết kế truyền thống giữa kỹ sư kết cấu (coi móng là ngàm cứng) và kỹ sư địa kỹ thuật (coi kết cấu phần trên là tải trọng tĩnh đặt lên gối mềm).
  • Trình bày điều kiện cần xét đến bài toán SSI dựa trên vận tốc sóng cắt $V_s$, chiều cao công trình $h$, và tần số dao động cơ bản $f$: $V_s / (f \cdot h) \le 200$, tương đương với $V_s \le 600\text{ m/s}$ đối với nhà nhiều tầng thông thường hoặc $V_s \le 57580/\lambda_s$ đối với kết cấu thép dạng công-xon.
  • Phân biệt bản chất của tương tác động học (kinematic interaction - sự tán xạ, khúc xạ sóng tới qua móng) và tương tác quán tính (inertial interaction - lực quán tính từ khối lượng phần trên truyền ngược xuống nền).
  • Làm rõ cơ chế hình thành hai dạng phi tuyến tại vùng tiếp xúc: phi tuyến hình học (nhổ móng/uplift khi mô men lật vượt giới hạn) và phi tuyến vật liệu (chảy dẻo nền đất, tập trung ứng suất).
  • Đánh giá ưu nhược điểm của các mô hình phần tử vĩ mô kinh điển (Nova, Cremer, Paolucci, Grange, Chatzigogos, Figini), xác định cơ sở để luận án kế thừa và phát triển mô hình tính toán.

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ô

Chương 2 trình bày chi tiết việc thiết lập mô hình toán học và giải thuật phần tử vĩ mô:

  • Xây dựng vector lực thu gọn và vector chuyển vị thu gọn không thứ nguyên cho móng phẳng tiếp xúc nền cát.
  • Xác định biểu thức tính sức chịu nén cực hạn $q_{max}$ theo lý thuyết Terzaghi có hiệu chỉnh theo hình dạng và chiều sâu chôn móng $D_f$.
  • Thiết lập hàm dẻo $f(\mathbf{F})$ định nghĩa mặt giới hạn chịu tải trong không gian lực $(N, V, M)$ và thế dẻo $g(\mathbf{F})$ điều khiển quy luật biến dạng dẻo trượt và xoay.
  • Xây dựng thuật toán xác định ma trận độ cứng tiếp tuyến của phần tử vĩ mô xét đến việc suy giảm diện tích tiếp xúc khi móng bị nhổ lên khỏi mặt đất trong từng bước gia tải động đất.
  • Thiết lập hệ phương trình vi phân chuyển động cho 2 trường hợp:
    1. Hệ móng - đất nền đơn lẻ: Mô hình 2 bậc tự do (chuyển vị ngang $u$ và góc xoay $\theta$).
    2. Hệ kết cấu phần trên - móng - đất nền: Mô hình thông số tập trung kết hợp khối lượng kết cấu phần trên, trụ cầu/cột đàn hồi và phần tử vĩ mô tại móng.
  • Lập sơ đồ khối thuật toán và chương trình tính toán số tích phân Newmark trong môi trường MATLAB, kiểm tra tính ổn định thông qua ví dụ kiểm chứng với phần mềm Cyclic1D/CyclicTP.

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 bàn rung:

  • Thiết kế mô hình thu nhỏ: Sử dụng các hệ số tỷ lệ hình học $\lambda$ theo lý thuyết của Meymand (1998), tham khảo kích thước công trình cầu nguyên mẫu Sanguinneti (Sonora, California, Hoa Kỳ) từ các nghiên cứu của Drosos (2012) và Anastasopoulos (2012, 2014).
  • Chế tạo thiết bị và chuẩn bị mẫu:
    • Hộp đất vách cứng được gia công và liên kết chắc chắn vào mặt bàn rung R202(UTC).
    • Đất nền là cát đồng nhất, đầm nén theo từng lớp để kiểm soát độ chặt và dung trọng khối lượng tích thể $\gamma$.
    • Móng mô hình bằng vật liệu kim loại có độ cứng cao; trụ kết cấu chế tạo từ thép hình đảm bảo làm việc trong miền đàn hồi; khối lượng kết cấu phần trên được gắn tập trung tại đỉnh trụ.
  • Bố trí cảm biến đo: Cảm biến gia tốc đặt tại sàn bàn rung, đáy hộp cát, đỉnh móng và đỉnh khối lượng kết cấu phần trên; các cảm biến đo chuyển vị LVDT bố trí theo phương ngang để ghi nhận dao động tương đối.
  • Kịch bản thí nghiệm:
    • Nhóm chuỗi thí nghiệm hệ móng - đất nền (T12, T13, T14, T15): Khảo sát gia tốc và chuyển vị đỉnh móng dưới các mức biên độ gia tốc kích thích khác nhau.
    • Nhóm chuỗi thí nghiệm hệ kết cấu phần trên - móng - đất nền (T21, T22, T23, T24, T25, T26): Khảo sát phản ứng động của toàn hệ với các điều kiện độ sâu chôn móng khác nhau (đặt trên mặt đất $D_f = 0$, chôn sâu $50%$ chiều cao móng, và chôn sâu $100%$ chiều cao móng) dưới tác động của phổ gia tốc động đất Tolmezzo.
Nhóm thí nghiệm Ký hiệu mẫu Đối tượng khảo sát Kích thích / Tình trạng chôn móng Thông số ghi nhận chính
Móng - Đất nền T12, T13, T14, T15 Móng đặt trên nền cát Tăng dần mức biên độ gia tốc kích thích Gia tốc cực đại đỉnh móng, độ lún
Kết cấu - Móng - Đất nền T21, T22, T23 Hệ toàn thể (LPM + móng) Móng đặt trên mặt đất ($D_f = 0$) Chuyển vị ngang, gia tốc đỉnh kết cấu
Kết cấu - Móng - Đất nền T24, T25, T26 Hệ toàn thể (LPM + móng) Chôn móng $0%$, $50%$ và $100%$ chiều cao Chuyển vị, gia tốc đỉnh, hiện tượng tách móng

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 việc đối chiếu kết quả mô phỏng số với số liệu đo thực nghiệm và phân tích độ nhạy tham số:

  • Xác định các thông số tương đương của mô hình thực nghiệm (độ cứng đàn hồi, khối lượng, hệ số cản, góc ma sát trong, dung trọng đất) đưa vào chương trình giải tích.
  • So sánh kết quả phân tích lý thuyết với thí nghiệm:
    • Hệ móng - đất nền: So sánh lịch sử thời gian và giá trị gia tốc cực đại tại đỉnh móng cho các trường hợp T12-00, T13-00, T14-00, T15-00; đánh giá sai số đỉnh cực đại và độ trễ pha dao động.
    • Hệ kết cấu phần trên - móng - đất nền: So sánh dạng sóng chuyển vị và gia tốc đỉnh kết cấu giữa mô hình đề xuất của luận án, mô hình của Paolucci (2008) và dữ liệu thực nghiệm bàn rung (các trường hợp T21-00 đến T26-00).
  • Phân tích ảnh hưởng của độ cứng đất nền $K_S$:
    • Khảo sát sự biến thiên của phản ứng động học kết cấu khi độ cứng nền thay đổi ở các mức: $1% k_1$, $5% k_1$, $10% k_1$, và $50% k_1$ (với $k_1$ là độ cứng đàn hồi tham chiếu).
    • Kết quả cho thấy khi độ cứng nền $K_S$ giảm (nền mềm hơn), chu kỳ dao động riêng của toàn hệ kéo dài, làm giảm lực quán tính và gia tốc đỉnh truyền lên kết cấu phần trên, tuy nhiên chuyển vị ngang tổng cộng tại đỉnh kết cấu lại tăng lên đáng kể do chuyển vị xoay của móng đóng góp vào.

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

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

  1. Luận án là công trình nghiên cứu hệ thống đầu tiên tại Việt Nam tiếp cận và chuẩn hóa thuật ngữ khoa học "phần tử vĩ mô" (macro-element) trong bài toán tương tác động lực học kết cấu - đất nền (SSI).
  2. Xây dựng hoàn chỉnh mô hình tính toán phần tử vĩ mô 2D tích hợp ma trận độ cứng tiếp tuyến có khả năng mô phỏng đồng thời hai hiệu ứng phi tuyến quan trọng: phi tuyến hình học (hiện tượng móng tách rời - nhổ lên khỏi mặt đất) và phi tuyến vật liệu của đất nền dưới đáy móng.
  3. Thiết lập hệ phương trình vi phân chuyển động phi tuyến kết hợp giữa mô hình tham số tập trung (LPM) của kết cấu phần trên và phần tử vĩ mô của móng - đất, giải bằng thuật toán tích phân Newmark trên MATLAB với thời gian tính toán nhanh vượt trội so với mô hình phần tử hữu hạn chi tiết nhưng vẫn đảm bảo độ tin cậy.

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

  1. Cung cấp bộ cơ sở dữ liệu thực nghiệm bàn rung phong phú (với các chuỗi thí nghiệm T12–T15 và T21–T26) về gia tốc và chuyển vị của móng nông và hệ kết cấu - móng đặt trên nền cát chịu kích thích động đất thực tế.
  2. Kiểm chứng và khẳng định tính khả thi của mô hình phần tử vĩ mô: sai số về giá trị cực đại và dạng sóng dao động giữa mô phỏng giải tích và thực nghiệm bàn rung nằm trong phạm vi chấp nhận được của cơ học công trình.
  3. Làm sáng tỏ định lượng vai trò của độ cứng tương tác $K_S$ và mức độ chôn sâu của móng ($D_f$) đối với việc phân bổ gia tốc và biến dạng ngang của kết cấu phần trên, cung cấp cơ sở khoa học để các kỹ sư thiết kế công trình giao thông và dân dụng cân nhắc hiệu ứng SSI thay vì mặc định giả thiết ngàm cứng.

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 dừng lại ở bài toán phẳng 2D, áp dụng cho móng nông đặt trên nền cát chặt đồng nhất, chưa xét đến đất dính (sét) hoặc nền đất phân tầng phức tạp.
  • Chưa xét đến cặp tương tác chéo giữa chuyển vị ngang và góc xoay trong ma trận độ cứng của phần tử vĩ mô.
  • Thí nghiệm bàn rung sử dụng hộp đất vách cứng chịu kích thích một phương ngang đơn thuần; chưa xét đến ảnh hưởng của sóng động đất đa chiều hoặc sự biến đổi dần đặc tính cơ lý của đất (hóa lỏng, suy giảm độ chặt) sau nhiều chu kỳ rung liên tiếp.

Hướng nghiên cứu tiếp

  • Mở rộng xây dựng phần tử vĩ mô không gian 3D chịu tải trọng động đất tác dụng đồng thời theo ba phương.
  • Nghiên cứu mô hình phần tử vĩ mô cho móng cọc, móng giếng chìm và móng bè trên nền đất yếu có xét đến hiện tượng tăng áp lực nước lỗ rỗng và hóa lỏng đất.
  • Hoàn thiện tích hợp thuật toán phần tử vĩ mô thành các mô-đun phần tử biên/phần tử đặc biệt trong các phần mềm kết cấu chuyên dụng thương mại hoặc mã nguồn mở.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và học viên sau đại học: Tài liệu cung cấp tổng quan tài liệu chuyên sâu về cơ học tính toán động đất, phương pháp mô hình hóa phần tử vĩ mô, và kỹ thuật lập trình tích phân Newmark giải phương trình vi phân chuyển động phi tuyến.
  • Đối với giảng viên và cơ sở nghiên cứu: Luận án là tài liệu tham khảo giá trị cho các môn học chuyên đề: Động lực học công trình, Tương tác kết cấu - môi trường, Kỹ thuật động đất tại các trường đại học khối kỹ thuật xây dựng và giao thông; đồng thời cung cấp phương pháp luận thiết kế thí nghiệm bàn rung mô hình thu nhỏ theo định luật tỷ lệ hình học.
  • Đối với kỹ sư thiết kế kết cấu và địa kỹ thuật: Cung cấp phương pháp đánh giá nhanh ảnh hưởng của nền đất mềm và hiệu ứng tách móng đến chuyển vị đỉnh công trình, giúp tối ưu hóa kích thước móng và bố trí chi tiết kháng chấn hợp lý hơn so với giả thiết ngàm cứng truyền thống.

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

1. Phần tử vĩ mô (macro-element) khác gì so với phương pháp phần tử hữu hạn trực tiếp (FEM) trong phân tích SSI?

Phương pháp phần tử hữu hạn trực tiếp đòi hỏi chia lưới toàn bộ thể tích đất và kết cấu thành hàng ngàn đến hàng triệu phần tử vi mô với các liên kết nút phức tạp, đòi hỏi thời gian mô hình hóa và tài nguyên tính toán rất lớn. Ngược lại, phần tử vĩ mô thu gọn toàn bộ ứng xử phức tạp của móng và môi trường đất xung quanh về một nút phần tử duy nhất đặt tại đáy móng, nhưng vẫn tích hợp đầy đủ các quy luật phi tuyến hình học (uplift) và vật liệu thông qua hàm thế dẻo và mặt tải giới hạn, giúp giảm đáng kể thời gian tính toán.

2. Tiêu chí nào xác định một công trình bắt buộc phải xét đến hiệu ứng tương tác kết cấu - đất nền (SSI)?

Theo công thức kinh điển được luận án trích dẫn, bài toán SSI cần được xem xét khi tỷ số: $$\frac{V_s}{f \cdot h} \le 200$$ Trong đó $V_s$ là vận tốc sóng cắt của đất, $f$ là tần số dao động cơ bản của kết cấu khi ngàm cứng, và $h$ là chiều cao công trình. Đối với công trình nhà dân dụng nhiều tầng, điều này tương ứng với điều kiện nền đất có vận tốc sóng cắt $V_s \le 600\text{ m/s}$. Đối với kết cấu dạng công-xon (trụ cầu, tháp), tiêu chuẩn áp dụng 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à tỷ số mảnh.

3. Phi tuyến hình học và phi tuyến vật liệu của hệ móng - đất nền trong luận án được mô tả như thế nào?

  • Phi tuyến hình học: Xuất hiện khi mô men lật lớn làm móng xoay, khiến một phần đáy móng bị nhổ lên khỏi mặt tiếp xúc với đất (uplift). Diện tích làm việc tiếp xúc bị thu hẹp từ bề rộng $B$ xuống $B' = B(1 - \delta)$, làm thay đổi độ cứng của hệ.
  • Phi tuyến vật liệu: Thể hiện qua trạng thái chảy dẻo của đất dưới đáy móng khi ứng suất đạt tới sức chịu tải cực hạn $q_{max}$, được mô tả bằng hàm dẻo $f(\mathbf{F})$ và quy luật chảy dẻo $g(\mathbf{F})$.

4. Kết quả nghiên cứu chỉ ra độ cứng đất nền $K_S$ ảnh hưởng như thế nào đến phản ứng động học của kết cấu?

Khi độ cứng đất nền $K_S$ giảm (đất nền mềm hơn hoặc độ sâu chôn móng nông):

  • Tần số dao động riêng của toàn hệ giảm (chu kỳ kéo dài), dẫn đến gia tốc và lực quán tính tác dụng lên kết cấu phần trên giảm đi.
  • Tuy nhiên, góc xoay tại móng tăng lên, làm cho chuyển vị ngang tổng cộng tại đỉnh kết cấu phần trên tăng lên đáng kể so với mô hình liên kết cứng.

5. Thí nghiệm bàn rung trong luận án sử dụng thiết bị và quy chuẩn tỷ lệ nào?

Thí nghiệm được thực hiện trên bàn rung một phương R202(UTC) với hộp đất vách cứng. Tỷ lệ mô hình thu nhỏ được thiết kế dựa trên lý thuyết tỷ lệ mô hình hóa trường trọng lực đơn của Meymand (1998) với tham số tỷ lệ hình học $\lambda$, tham chiếu theo kích thước hình học của công trình cầu nguyên mẫu Sanguinneti (Sonora, California, Hoa Kỳ).

Kết luận

Luận án tiến sĩ kỹ thuật của tác giả Huỳnh Văn Quân đã giải quyết thành công bài toán phân tích tương tác động lực học kết cấu - đất nền chịu động đất thông qua việc đề xuất và ứng dụng mô hình phần tử vĩ mô 2D tích hợp đồng thời phi tuyến hình học và vật liệu. Sự kết hợp chặt chẽ giữa mô hình giải tích thông số tập trung giải bằng thuật toán Newmark và chuỗi thí nghiệm kiểm chứng trên bàn rung R202(UTC) đã khẳng định tính chính xác và hiệu quả vượt trội về mặt khối lượng tính toán của phương pháp đề xuất. Công trình đóng góp nền tảng lý luận và thực nghiệm quan trọng, mở ra hướng ứng dụng tính toán kháng chấn thực tế cho các công trình xây dựng và cầu đường tại Việt Nam.