Phân Tích Tương Tác Nền Móng và Kết Cấu Nhà Cao Tầng: Bước Tiến Trong Mô Hình Hóa Bằng Phương Pháp Phần Tử Hữu Hạn (FEM)

Tóm tắt nghiên cứu (Research Summary)

Nghiên cứu "Phân tích tương tác nền móng và kết cấu bên trên công trình nhà cao tầng bằng phương pháp phần tử hữu hạn" (Mã đề tài: SV2021-74, Đại học Sư phạm Kỹ thuật TP.HCM) tập trung giải quyết câu hỏi then chốt: Mô hình ngàm cứng truyền thống tại chân cột/vách có phản ánh chính xác ứng xử thực tế của nhà cao tầng, và sự biến đổi động học, nội lực diễn ra như thế nào khi kể đến sự làm việc đồng thời của hệ Đất nền – Móng – Kết cấu bên trên (Soil-Structure Interaction - SSI)?

Bằng việc áp dụng phương pháp phần tử hữu hạn (FEM) trên phần mềm Robot Structural Analysis Professional (RSAP) của Autodesk, nhóm tác giả đã xây dựng mô hình kết cấu khung - vách không gian của công trình 15 tầng (Chung cư Lý Thường Kiệt, cao 52.7m). Hệ cọc bê tông cốt thép 400x400mm cùng đất nền được mô phỏng chi tiết bằng các phần tử thanh và hệ lò xo đàn hồi có độ cứng biến thiên theo chiều sâu, dựa trên lý thuyết Joseph E. Bowles và dữ liệu địa chất SPT thực tế.

Kết quả nghiên cứu chỉ ra rằng: việc xét đến tương tác đất – kết cấu làm giảm độ cứng tổng thể của công trình, kéo dài chu kỳ dao động riêng ($T$), gia tăng đáng kể chuyển vị đỉnh và chuyển vị lệch tầng dưới tác động của tải trọng gió và động đất. Đồng thời, nội lực trong các cấu kiện dầm và cột ở các tầng dưới xảy ra hiện tượng phân phối lại rõ rệt. Nghiên cứu cung cấp cơ sở khoa học và quy trình thực hành thiết kế tin cậy, giúp tối ưu hóa an toàn và chi phí cho các công trình cao tầng hiện đại.

Bối cảnh và tầm quan trọng (Context & Significance)

Trong bối cảnh đô thị hóa nhanh chóng, các công trình nhà cao tầng mọc lên ngày càng nhiều với quy mô và chiều cao liên tục gia tăng. Để đảm bảo tính an toàn và tối ưu kinh tế, việc mô phỏng chính xác ứng xử chịu lực của công trình dưới tác động của tải trọng trọng trường cùng các tải trọng ngang nguy hiểm (gió động, động đất) là nhiệm vụ sống còn trong kỹ thuật kết cấu.

Tuy nhiên, thực tế thiết kế kết cấu tại Việt Nam hiện nay vẫn chủ yếu áp dụng giả thiết đơn giản hóa: tách rời kết cấu bên trên và hệ móng bên dưới. Trong đó, chân cột và chân vách được giả định là các liên kết ngàm cứng tuyệt đối ($Rigid\ Base$). Phản lực chân cột sau khi giải bài toán phần thân sẽ được gán làm tải trọng đầu vào độc lập để tính toán móng bằng các bảng tính giải tích hoặc phần mềm tách rời (như SAFE).

Giả thiết ngàm cứng này đã bỏ qua hoàn toàn:

  1. Độ lún đàn hồi và độ lún lệch của nền đất xung quanh cọc.
  2. Khả năng xoay và chuyển vị ngang của đài cọc.
  3. Sự tương tác động lực học liên tục giữa môi trường địa kỹ thuật và hệ kết cấu thân.

Khoảng trống nghiên cứu và thực hành này dẫn đến việc đánh giá sai lệch chu kỳ dao động tự nhiên, bỏ sót các xung lực nguy hiểm do hiệu ứng mềm hóa nền móng, hoặc phân bố thép cục bộ không sát với trạng thái làm việc thực. Do đó, việc ứng dụng mô hình tương tác đồng thời Nền – Móng – Thân bằng phương pháp phần tử hữu hạn là một đòi hỏi cấp thiết, kịp thời và mang tính đột phá nhằm tiệm cận với các tiêu chuẩn thiết kế tiên tiến trên thế giới (như Eurocode 8, ASCE 7).

Methodology và approach (Phương pháp tiếp cận & Mô hình hóa)

Đề tài triển khai phương pháp nghiên cứu kết hợp giữa cơ sở lý thuyết cơ học đất - nền móng, mô phỏng số 3D bằng phương pháp phần tử hữu hạn và phân tích so sánh định lượng trực tiếp giữa hai mô hình:

  • Mô hình 1 (Tách rời): Chân hệ chịu lực đứng được ngàm cứng tại mặt đài móng.
  • Mô hình 2 (Đồng thời SSI): Toàn bộ kết cấu phần thân, đài móng, cọc và đất nền xung quanh được mô hình hóa trong cùng một không gian tính toán.
+-------------------------------------------------------------------+
|                        QUY TRÌNH MÔ HÌNH HÓA                       |
+-------------------------------------------------------------------+
| 1. Khảo sát địa chất & Xác định sức chịu tải cọc (SPT, cường độ)  |
| 2. Tính hệ số nền Cz & Độ cứng lò xo đất theo Joseph E. Bowles    |
| 3. Xây dựng mô hình 3D trong Robot Structural Analysis (RSAP)     |
|    - Cột, dầm, cọc: Phần tử thanh (Bar elements)                  |
|    - Sàn, vách, đài móng: Phần tử tấm/vỏ (Plate/Shell elements)   |
|    - Đất nền: Lò xo đàn hồi (Elastic Springs) phi tuyến/lặp       |
| 4. Phân tích Modal, Gió động (TCVN 2737), Động đất (TCVN 9386)    |
| 5. Phân tích lặp hội tụ độ cứng lò xo và xuất kết quả so sánh    |
+-------------------------------------------------------------------+

Thiết lập thông số và công cụ tính toán

  • Công trình mẫu: Tòa nhà chung cư 15 tầng (+52.700m), kích thước mặt bằng 25.2m x 21.0m, hệ kết cấu bê tông cốt thép cấp độ bền B35/B30 kết hợp vách cứng lõi thang máy.
  • Mô hình hóa phần tử: Dầm, cột và cọc vuông 400x400mm sử dụng phần tử thanh (Bar). Sàn dày 150mm, vách dày 300mm và đài móng sử dụng phần tử tấm vỏ (Plate/Shell).
  • Xác định độ cứng lò xo đất: Sử dụng công thức bán thực nghiệm của Joseph E. Bowles: $$C_z = A_s + B_s \cdot Z^n$$ Trong đó, $A_s$ và $B_s$ phụ thuộc vào các chỉ tiêu cơ lý của từng lớp đất (dung trọng $\gamma$, lực dính $c$, góc ma sát trong $\varphi$, hệ số hình dạng móng $S_c, S_\gamma$, và hệ số sức chịu tải $N_c, N_q, N_\gamma$). Độ cứng của gối đàn hồi tại nút thứ $i$ được tính bằng: $$k_i = C_{zi} \cdot A_i$$
  • Tính toán tải trọng tác động: Tĩnh tải, hoạt tải theo TCVN 2737:1995; Thành phần động của tải trọng gió được tính toán theo phương pháp chu kỳ dao động riêng cho từng phương X, Y; Tải trọng động đất được phân tích bằng phương pháp phổ phản ứng đàn hồi đa mode (Response Spectrum Analysis) tương ứng cấp động đất VII (thang MSK-64, $a_g = 0.0774g$).

Quy trình phân tích lặp được thực hiện trên RSAP cho đến khi phản lực, độ lún và độ cứng lò xo đất đạt sự hội tụ ổn định, đảm bảo tính chuẩn xác và độ tin cậy cơ học cao.

Phát hiện chính (Key Findings)

Nghiên cứu đã đưa ra các kết quả định lượng cụ thể khi so sánh giữa Mô hình ngàm cứng (Mô hình 1) và Mô hình tương tác đồng thời (Mô hình 2):

1. Kéo dài chu kỳ dao động và làm mềm hệ kết cấu

Khi xét đến tính biến dạng đàn hồi của nền đất, độ cứng tổng thể của toàn công trình bị suy giảm:

  • Chu kỳ dao động riêng cơ bản ($T_1$) tăng lên đáng kể ở tất cả các dạng dao động (Mode 1, Mode 2, Mode 3).
  • Tần số dao động tự nhiên ($f$) giảm tương ứng, phản ánh trạng thái công trình dao động "mềm" hơn so với giả định ngàm cứng lý tưởng.

2. Biến động tải trọng ngang (Gió động & Động đất)

  • Tải trọng gió động: Do tần số riêng $f$ giảm, hệ số động lực $\xi$ và lực quán tính thay đổi, dẫn đến thành phần động của tải trọng gió phân bố lại theo chiều cao công trình so với mô hình ngàm cứng.
  • Tải trọng động đất: Khi chu kỳ dao động $T$ tăng, điểm làm việc trên phổ phản ứng thiết kế đàn hồi ($S_d(T)$) dịch chuyển về phía chu kỳ dài. Khối lượng tham gia dao động hữu hiệu ở các mode cơ bản thay đổi, làm biến đổi lực cắt đáy tổng cộng ($F_b$) và lực động đất phân bố tại từng tầng.

3. Gia tăng chuyển vị ngang đỉnh và chuyển vị lệch tầng

Đây là một trong những phát hiện quan trọng nhất về mặt an toàn kết cấu:

  • Chuyển vị đỉnh công trình: Dưới tác dụng của tải trọng gió và động đất, chuyển vị ngang đỉnh ở Mô hình 2 tăng vượt bậc so với Mô hình 1 (do cộng hưởng chuyển vị biến dạng thân và góc xoay/chuyển vị ngang tại đài móng).
  • Chuyển vị lệch tầng ($\Delta/h$): Tăng mạnh tại các tầng dưới (từ tầng hầm đến tầng 4), đặt ra yêu cầu kiểm tra nghiêm ngặt điều kiện ổn định P-Delta và khả năng nứt của các cấu kiện bao che, vách ngăn kiến trúc.
+---------------------------------------------------------------------+
|         BẢNG SO SÁNH TỔNG QUAN GIỮA HAI MÔ HÌNH TÍNH TOÁN           |
+------------------------------+------------------+-------------------+
| Tiêu chí so sánh             | Mô hình 1 (Ngàm) | Mô hình 2 (SSI)   |
+------------------------------+------------------+-------------------+
| Độ cứng tổng thể             | Cao hơn          | Giảm (Mềm hóa)    |
| Chu kỳ dao động riêng (T)    | Ngắn hơn         | Kéo dài hơn       |
| Chuyển vị đỉnh công trình    | Bé hơn           | Tăng rõ rệt       |
| Chuyển vị lệch tầng dưới     | Nhỏ              | Tăng đáng kể      |
| Phân phối mô-men uốn dầm     | Cố định theo độ  | Phân phối lại     |
|                              | cứng tiết diện   | do lún lệch cọc   |
| Lực dọc cột biên và lõi      | Theo diện tải    | Tái phân phối     |
+------------------------------+------------------+-------------------+

4. Tái phân phối nội lực trong dầm và cột

  • Cấu kiện Dầm: Giá trị mô-men uốn tại gối và nhịp của các dầm tầng dưới (như dầm B1, B2, B3 từ tầng 1 đến tầng 3) thay đổi rõ rệt do ảnh hưởng của chuyển vị lún không đều giữa các cụm cọc dưới cột và vách.
  • Cấu kiện Cột: Lực dọc ($N$) và mô-men uốn ($M$) tại chân cột có sự chuyển dịch tải trọng: một số cột chịu tải lớn hơn dự kiến trong khi các cột lân cận giảm tải, phản ánh đúng cơ chế tự điều tiết độ cứng của hệ kết cấu không gian.

Đóng góp khoa học (Scientific Contributions)

Lĩnh vực đóng góp Chi tiết giá trị khoa học và thực tiễn
Đóng góp lý thuyết Hệ thống hóa cơ sở cơ học tiếp xúc và tương tác Đất – Móng – Kết cấu (SSI); làm sáng tỏ mối liên hệ giữa các chỉ tiêu cơ lý đất nền (SPT, $c$, $\varphi$) với ma trận độ cứng động học của công trình cao tầng.
Cải tiến phương pháp Đề xuất quy trình tính toán độ cứng lò xo đất theo mô hình Joseph E. Bowles kết hợp phương pháp lặp hội tụ trên phần mềm Robot Structural Analysis (RSAP), mở ra hướng ứng dụng mới thay thế cho các phần mềm truyền thống.
Ứng dụng thực tiễn Cung cấp tài liệu tham khảo trực quan giúp các kỹ sư kết cấu nhận diện nguy cơ mất an toàn khi sử dụng giả thiết ngàm cứng; hỗ trợ tối ưu hóa kích thước cọc, đài móng và hàm lượng cốt thép cấu kiện.
Định hướng tiêu chuẩn Đóng góp luận cứ kỹ thuật cho việc chuẩn hóa quy trình phân tích phi tuyến/đàn hồi tương tác nền móng trong hệ thống Tiêu chuẩn Xây dựng Việt Nam (TCVN).

Đối tượng quan tâm (Target Audience)

  • Kỹ sư Kết cấu & Đơn vị Tư vấn Thiết kế: Cung cấp phương pháp luận chuẩn xác để lập mô hình kết cấu cho các dự án nhà cao tầng, công trình phức bách thảo, các khu vực có địa chất yếu hoặc vùng có rủi ro động đất.
  • Nhà nghiên cứu, Giảng viên & Học viên Cao học/Sinh viên ngành Kỹ thuật Xây dựng: Nguồn tài liệu học thuật giá trị về ứng dụng phương pháp phần tử hữu hạn nâng cao, tính toán động lực học công trình và cơ học đất ứng dụng.
  • Chủ đầu tư & Các Ban Quản lý Dự án: Hiểu rõ tầm quan trọng của việc đầu tư cho công tác khảo sát địa chất và phân tích kết cấu chuyên sâu, giúp giảm thiểu rủi ro lún nứt công trình và tối ưu hóa chi phí vật liệu xây dựng.

Các câu hỏi thường gặp (FAQ)

1. Phát hiện quan trọng nhất của nghiên cứu này là gì?

Phát hiện quan trọng nhất là: Mô hình ngàm cứng truyền thống đánh giá thiếu an toàn về mặt chuyển vị (chuyển vị đỉnh và chuyển vị lệch tầng bị đánh giá thấp hơn thực tế) và sai lệch về mặt phân phối nội lực ở các tầng dưới cùng của nhà cao tầng do bỏ qua độ lún đàn hồi và tương tác cọc - đất nền.

2. Phương pháp xác định độ cứng lò xo theo Joseph E. Bowles có gì ưu việt hơn TCVN 10304:2014?

Công thức của TCVN 10304:2014 phụ thuộc nhiều vào hệ số tỷ lệ $k$ quy ước với dải tra bảng biến động rất rộng (dễ gây sai số lớn). Trong khi đó, phương pháp của Joseph E. Bowles tích hợp trực tiếp các chỉ tiêu cơ lý đất nền cụ thể ($c$, $\varphi$, $\gamma$, ứng suất hữu hiệu $q'$) và hình học cọc, mang lại độ tin cậy và phản ánh sát thực tế địa chất hơn.

3. Kết quả của nghiên cứu có thể khái quát hóa cho mọi loại nhà cao tầng không?

Quy luật chung về sự "mềm hóa" kết cấu (tăng chu kỳ $T$, tăng chuyển vị) có thể khái quát hóa cho hầu hết các công trình nhà cao tầng móng cọc. Tuy nhiên, mức độ biến thiên nội lực cụ thể sẽ phụ thuộc trực tiếp vào địa tầng đất nền, sơ đồ mặt bằng kiến trúc, tỷ lệ độ cứng giữa khung và vách lõi của từng công trình.

4. Hướng phát triển tiếp theo của đề tài là gì?

Nghiên cứu có thể mở rộng sang việc khảo sát ứng xử phi tuyến vật liệu của đất nền (sử dụng mô hình Mohr-Coulomb hoặc Hardening Soil trong Plaxis 3D liên kết với RSAP), phân tích cọc chịu tải trọng ngang lớn kết hợp tải trọng động đất theo miền thời gian (Time-history analysis).

5. Có thể áp dụng ngay quy trình này vào thiết kế thực tế hay không?

Hoàn toàn có thể. Quy trình tính toán lò xo đất từ dữ liệu SPT và địa chất hiện trường kết hợp mô hình hóa trên các phần mềm FEM hiện đại (RSAP, ETABS) đã được chuẩn hóa từng bước trong đề tài, giúp kỹ sư dễ dàng tích hợp vào quy trình thiết kế hàng ngày.

Kết luận (Conclusion)

Nghiên cứu của nhóm tác giả đã khẳng định vai trò quyết định của việc xét đến tương tác Nền móng – Đất nền – Kết cấu bên trên trong phân tích ứng xử công trình cao tầng. Việc loại bỏ giả thiết ngàm cứng lý tưởng hóa để chuyển sang mô hình tương tác đồng thời bằng phương pháp phần tử hữu hạn là bước đi tất yếu, giúp nâng cao độ tin cậy, bảo đảm an toàn công trình trước tải trọng gió động và động đất, đồng thời tối ưu hóa chi phí đầu tư.

Các kỹ sư và nhà thiết kế được khuyến nghị ứng dụng mô hình SSI vào các dự án thực tế để đảm bảo chất lượng công trình tiệm cận với các tiêu chuẩn xây dựng quốc tế hiện đại.