Giới thiệu dự án

Trong bối cảnh đô thị hóa nhanh chóng tại Việt Nam, đặc biệt là khu vực Đồng bằng sông Cửu Long và TP. Hồ Chí Minh với đặc trưng địa tầng trầm tích trẻ, lớp bùn sét yếu phân bố dày từ 15m đến hơn 40m, công tác khảo sát và thiết kế địa cơ nền móng đóng vai trò quyết định đến 40–50% độ an toàn và chi phí xây dựng công trình. Các sự cố công trình như lún lệch quá mức, nứt gãy dầm móng hay sự cố mất ổn định sức chịu tải cọc thường bắt nguồn từ việc xử lý số liệu địa chất thiếu chính xác và mô hình hóa tương tác giữa kết cấu móng với nền đất (Soil-Structure Interaction - SSI) chưa phản ánh đúng bản chất vật lý.

Đồ án "Tính toán & Thiết kế Nền Móng Công trình" do sinh viên Nguyễn Thành Trí (MSSV: 1814466, Lớp Kỹ thuật Xây dựng, Trường Đại học Bách Khoa – ĐHQG TP.HCM) thực hiện dưới sự hướng dẫn của TS. Lê Trọng Nghĩa tập trung giải quyết bài toán phức hợp: chuẩn hóa thống kê khối lượng mẫu địa chất khổng lồ từ 4 hố khoan sâu 80m – 100m, sau đó thiết kế tối ưu hai hệ kết cấu móng điển hình là Móng băng giao thoa trên nền đàn hồi và Móng cọc đài thấp chịu lực phức tạp.

graph TD
    A[Khảo sát 4 hố khoan LK1-LK4] --> B[Xử lý thống kê chuẩn TCVN / LINEST]
    B --> C[Phân tầng 5 lớp địa chất & Xác định chỉ tiêu cơ lý]
    C --> D[Thiết kế Móng Băng]
    C --> E[Thiết kế Móng Cọc BTCT]
    D --> D1[Mô hình dầm trên nền đàn hồi Winkler SAP2000]
    D --> D2[Kiểm tra ổn định, tính lún cộng lớp phân tố]
    E --> E1[Tính sức chịu tải Pvl, SPT, Chỉ tiêu cơ lý]
    E --> E2[Kiểm tra cẩu lắp, tải ngang & móng khối quy ước]

Mục tiêu kỹ thuật cụ thể của đồ án:

  1. Chuẩn hóa thống kê địa chất: Phân tích dữ liệu từ 4 hố khoan (LK1: 80m/40 mẫu, LK2: 80m/40 mẫu, LK3: 100m/50 mẫu, LK4: 100m/48 mẫu) qua 5 lớp địa tầng theo tiêu chuẩn quy phạm, loại trừ sai số thô bằng phân phối Student ($n < 25$) và phân phối chuẩn ($n \ge 25$).
  2. Xác định thông số chống cắt: Ứng dụng phương pháp hồi quy tuyến tính (hàm LINEST) để tìm cặp thông số lực dính kết ($c$) và góc ma sát trong ($\varphi$) theo hai trạng thái giới hạn (TTGH I và TTGH II).
  3. Mô hình hóa móng băng: Phân tích nội lực dầm móng băng liên tục trên nền lò xo đàn hồi Winkler bằng phần mềm SAP2000 v22, tính toán độ lún theo phương pháp cộng lớp phân tố ($S \le [S] = 8.0\text{ cm}$).
  4. Thiết kế móng cọc toàn diện: Xác định sức chịu tải cọc theo vật liệu ($P_{vl}$), xuyên tiêu chuẩn ($P_{spt}$) và chỉ tiêu cơ lý ($P_{cp}$); kiểm tra móng khối quy ước, xuyên thủng đài cọc, uốn - cắt trong quá trình cẩu lắp dựng và phản ứng chịu tải trọng ngang.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng và phương pháp luận

Tiêu chí phân tích Phương pháp giải tích cổ điển (Móng cứng tuyệt đối) Phương pháp phần tử hữu hạn 3D (PLAXIS 3D) Phương pháp đề xuất trong Đồ án (Winkler Beam + SAP2000 + LINEST)
Mô phỏng tương tác đất - móng Bỏ qua biến dạng uốn dầm móng; phản lực nền giả định tuyến tính đều Mô phỏng liên tục môi trường phi tuyến (Mohr-Coulomb, Hardening Soil); rất phức tạp Mô phỏng dầm uốn liên tục trên gối tựa đàn hồi ($k_z$), phản ánh chính xác phân bố moment
Khối lượng tính toán & Thời gian Nhanh nhưng sai số cục bộ chân cột lớn ($> 35%$) Rất tốn thời gian lập lưới, chi phí bản quyền cao Tối ưu hóa giữa độ chính xác kết cấu và tốc độ thiết kế kỹ thuật
Xử lý số liệu địa chất Lấy giá trị trung bình đại số, dễ sai lệch khi biến động Cần thông số mô hình nâng cao ($E_{50}, E_{ur}$) khó thí nghiệm Áp dụng bình sai thống kê chuẩn mực, loại trừ mẫu lệch theo tiêu chuẩn TCVN
Khả năng áp dụng thực tế Hạn chế cho nhà nhiều tầng, tải trọng lệch trục Phù hợp nghiên cứu chuyên sâu, hố đào sâu phức tạp Tiêu chuẩn hóa cho đồ án kỹ sư và tư vấn thiết kế kết cấu dân dụng

Ma trận ưu tiên yêu cầu thiết kế (MoSCoW):

  • Must have (Bắt buộc): Đáp ứng TTGH I về độ bền vật liệu và sức chịu tải nền đất ($P_{max} \le 1.2 R_{tc}$, $P_{tb} \le R_{tc}$); TTGH II về độ lún ($S \le 8\text{ cm}$, $\Delta S / L \le 0.002$).
  • Should have (Nên có): Tự động hóa gán lò xo phản lực nền $k_{zi} = C_{zi} \times b \times l_i$ vào các nút phần tử dầm trong SAP2000.
  • Could have (Có thể mở rộng): Kiểm tra điều kiện cẩu lắp cọc tại 2 điểm móc ($0.207L$) và 1 điểm dựng ($0.293L$).
  • Won't have (Chưa thực hiện): Mô hình hóa hiệu ứng nhóm cọc phi tuyến theo thời gian (tính cố kết thấm 3D phụ thuộc thời gian $t$).

Thiết kế hệ thống tính toán & Công nghệ áp dụng

  • Nền tảng phần mềm: CSI SAP2000 v22.0.0 (phân tích phần tử hữu hạn khung - dầm móng), Microsoft Excel 365 (Engine xử lý ma trận LINEST, tự động tính toán địa chất), AutoCAD 2022.
  • Hệ thống Tiêu chuẩn Kỹ thuật:
    • TCVN 9362:2012: Tiêu chuẩn thiết kế nền nhà và công trình.
    • TCVN 10304:2014: Tiêu chuẩn thiết kế móng cọc.
    • TCVN 5574:2018: Tiêu chuẩn thiết kế kết cấu bê tông và bê tông cốt thép.
    • TCVN 9351:2012: Đất xây dựng – Phương pháp thí nghiệm xuyên tiêu chuẩn (SPT).
classDiagram
    class GeologicalLayer {
        +String layerName
        +float gamma
        +float Gs
        +float W
        +float WL
        +float WP
        +float c
        +float phi
        +calculateStatisticalProperties()
    }
    class StripFoundation {
        +float length_L
        +float width_B
        +float height_h
        +float subgrade_k
        +calculateSettlement()
        +exportSAP2000Model()
    }
    class PileFoundation {
        +float pileSection
        +float pileLength
        +float P_material
        +float P_soil_SPT
        +checkPunchingShear()
        +checkLiftingStress()
    }
    GeologicalLayer --> StripFoundation : Cung cấp chỉ tiêu tính lún & k
    GeologicalLayer --> PileFoundation : Cung cấp ma sát thành fs & qp

Implementation và kết quả

Quy trình phân tích dữ liệu & Thuật toán cốt lõi

1. Xử lý thống kê và hồi quy thông số chống cắt

Đối với lớp 1 (Bùn sét, xám đen, trạng thái chảy - 41 mẫu thí nghiệm) và lớp 2 (Sét, xám đen, trạng thái chảy - 16 mẫu), chỉ tiêu cơ lý được xử lý qua thuật toán lọc bỏ giá trị ngoại lai:

$$\text{Điều kiện loại bỏ sai số:} \quad |A_i - A_{tb}| > \nu \cdot \sigma$$

Hệ số lực dính ($c$) và góc ma sát trong ($\varphi$) được xác định từ phương trình định luật Coulomb $\tau = \sigma \tan\varphi + c$ thông qua ma trận bình phương tối thiểu:

=LINEST(Tau_Array, Sigma_Array, TRUE, TRUE)

Kết quả tính toán cho Lớp 1 thu được:

  • Dung trọng tiêu chuẩn: $\gamma_{tc} = 14.64\text{ kN/m}^3$ (Độ lệch quân phương $\sigma = 0.337$, hệ số biến động $v = 0.050$).
  • Tỷ trọng hạt: $G_{s} = 2.007 \rightarrow G_{stc} = 2.62\text{ g/cm}^3$.
  • Độ ẩm tự nhiên: $W_{tb} = 89.599% \rightarrow W_{tc} = 89.00%$.
  • Giới hạn Atterberg: $W_L = 81.86%$, $W_P = 43.10% \rightarrow$ Chỉ số dẻo $I_P = W_L - W_P = 38.76%$, Độ sệt $B = (W - W_P)/I_P > 1.0$ (trạng thái chảy).
  • Thông số chống cắt trạng thái giới hạn II: $\varphi^{II}{tt} = [3^\circ18' - 3^\circ42']$, $c^{II}{tt} = 7.14\text{ kPa}$.

2. Thuật toán gán độ cứng nền lò xo Winkler cho Móng Băng

Độ cứng lò xo tại mỗi nút chia nhỏ ($Joint_i$) của dầm móng băng trong mô hình SAP2000 được rời rạc hóa:

$$K_{z,i} = C_{z} \times B_m \times \Delta L_i$$

Trong đó $C_z$ là hệ số nền ($kN/m^3$), $B_m$ là bề rộng cánh móng ($m$), và $\Delta L_i$ là khoảng cách phần tử chia nhỏ ($0.2m - 0.5m$).

def calculate_winkler_springs(beam_length, b_width, c_z, element_size=0.25):
    """Tính toán phân bố độ cứng lò xo Winkler cho dầm móng băng SAP2000"""
    num_elements = int(beam_length / element_size)
    num_nodes = num_elements + 1
    spring_stiffness = []
    
    for i in range(num_nodes):
        if i == 0 or i == num_nodes - 1:
            trib_length = element_size / 2.0  # Nút biên
        else:
            trib_length = element_size        # Nút trong
        
        kz = c_z * b_width * trib_length
        spring_stiffness.append((i, i * element_size, kz))
    
    return spring_stiffness

# Thực thi mẫu với dầm dài 18m, rộng 1.6m, Cz = 12000 kN/m3
springs = calculate_winkler_springs(18.0, 1.6, 12000, 0.25)
print(f"Tổng số nút lò xo: {len(springs)}, Độ cứng nút trong: {springs[1][2]} kN/m")

Tổng số nút lò xo: 73, Độ cứng nút trong: 4800.0 kN/m

Kết quả tính toán và thẩm tra thiết kế

Mô hình lực cắt (Q) và Moment (M) dầm móng băng trên SAP2000:
    Tải cột P1=1200kN       P2=1800kN       P3=1500kN
         ↓               ↓               ↓
   ================================================= (Dầm móng 0.4x0.8m)
   WWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWWW (Lò xo Winkler Kz)
   -------------------------------------------------
   Biểu đồ M: [~~~ Gối uốn âm ~~~]      [~~~ Nhịp uốn dương ~~~]
   Moment nhịp max: M_max = 348.5 kNm -> Chọn 4D22 lớp dưới
   Moment gối max: M_gối = -412.3 kNm -> Chọn 5D25 lớp trên

Bảng tổng hợp kết quả tính toán hai phương án móng:

Hạng mục kiểm tra kết cấu Giá trị tính toán thực tế Giới hạn tiêu chuẩn quy định Kết luận kỹ thuật
Ứng suất đáy móng băng ($P_{max}$) $142.6\text{ kPa}$ $P_{tc,cp} = 168.0\text{ kPa}$ Đạt an toàn ($\approx 84.8%$ công suất)
Độ lún móng băng ($S_{mb}$) $6.42\text{ cm}$ $[S] \le 8.00\text{ cm}$ Đạt yêu cầu TTGH II
Sức chịu tải cọc đơn $350\times350$ ($Q_a$) $780.0\text{ kN}$ $P_{vl} = 1250\text{ kN}; P_{spt} = 845\text{ kN}$ Chọn $Q_{tk} = 780\text{ kN}$ theo đất nền
Ứng suất đáy đài khối quy ước $210.5\text{ kPa}$ $R_{đáy} = 265.0\text{ kPa}$ Đảm bảo ổn định tổng thể
Độ lún đài cọc móng khối quy ước $2.85\text{ cm}$ $[S] \le 8.00\text{ cm}$ Kiểm soát biến dạng tốt
Khả năng chống đâm thủng đài cọc $P_{xuyên} = 620\text{ kN}$ $P_{chống} = 1150\text{ kN}$ Đạt hệ số an toàn $F_s = 1.85$
Ứng suất cẩu lắp cọc (2 điểm móc) $M_{cẩu} = 28.4\text{ kNm}$ $[M_{nứt}] = 42.0\text{ kNm}$ Không xuất hiện vết nứt trước khi ép

Đổi mới và đóng góp

  1. Chuẩn hóa quy trình tiền xử lý số liệu địa chất bằng hàm ma trận: Thay vì nhập liệu thủ công từng giá trị trung bình, tác giả đã xây dựng bảng tính tự động hóa bằng thuật toán LINEST kết hợp phân phối Student, giảm thiểu $90%$ nguy cơ sai sót do thao tác người dùng khi xử lý hơn 170 mẫu thí nghiệm cắt phẳng ($\sigma - \tau$).
  2. Khắc phục nhược điểm của phương pháp móng cứng: Ứng dụng mô hình dầm trên nền lò xo đàn hồi (Winkler Foundation) trên SAP2000 giúp phản ánh sự phân phối lại nội lực thực tế; giảm $14.2%$ lượng bê tông dầm móng và $11.8%$ lượng cốt thép dọc so với tính toán theo sơ đồ phân phối phản lực hình thang truyền thống.
  3. Phân tích toàn diện trạng thái thi công cọc: Bổ sung bước kiểm tra nội lực khi cẩu kéo tại công trường ($M = 0.0214 \cdot q \cdot L^2$) và cẩu dựng tại giá búa ($M = 0.043 \cdot q \cdot L^2$), ngăn ngừa nguy cơ nứt gãy thân cọc trước khi ép vào lòng đất.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng kỹ thuật

Quy trình và mô hình tính toán của đồ án có khả năng áp dụng trực tiếp cho các dự án xây dựng dân dụng từ 5 đến 15 tầng tại khu vực đất yếu:

  • Khu vực Nhà Bè, Quận 7, TP.HCM: Chiều dày bùn sét lớn $\rightarrow$ Ứng dụng mô hình móng cọc ép/khoan cắm vào tầng cát pha sỏi sạn (Lớp 3 và Lớp 4).
  • Khu vực tải trọng trung bình, móng công trình phụ trợ: Áp dụng giải pháp móng băng bản rộng kết hợp xử lý đệm cát hoặc cừ tràm giảm lún.
Lộ trình triển khai ứng dụng thiết kế:
[Tuần 1: Tiếp nhận hố khoan] -> [Tuần 2: Xử lý LINEST] -> [Tuần 3: Chạy mô hình SAP2000] -> [Tuần 4: Bản vẽ & Thi công]

Phân tích Hiệu quả Kinh tế (Cost-Benefit Analysis)

  • Chi phí thiết kế: Tự động hóa bảng tính giúp rút ngắn thời gian tính toán từ 45 giờ xuống còn 12 giờ làm việc của kỹ sư.
  • Chi phí vật tư: Tối ưu hóa chiều sâu chôn đài và tiết diện cọc ($350\times 350\text{ mm}$ thay vì $400\times 400\text{ mm}$ dự toán sơ bộ) giúp tiết kiệm xấp xỉ $120.000.000\text{ VNĐ}$ cho một khối công trình $500\text{ m}^2$ sàn.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Mô hình lò xo độc lập: Giả thiết nền Winkler bỏ qua sự liên kết trượt giữa các lò xo liền kề (hiệu ứng Pasternak), dẫn đến phản lực nền ngoài mép móng bị triệt tiêu đột ngột.
  • Tính toán lún đàn hồi: Phương pháp cộng lớp phân tố chưa xét đến yếu tố lún cố kết thứ cấp (lún dão từ biến) theo thời gian của đất sét yếu.

Hướng phát triển đề xuất

  1. Nâng cấp mô hình nền lò xo sang mô hình nền hai tham số (Pasternak/Vlasov) hoặc tích hợp mô hình phần tử hữu hạn tương tác đất - cọc 3D qua OpenSees / PLAXIS 3D.
  2. Viết Macro/Python Plugin tích hợp trực tiếp vào SAP2000 API để tự động trích xuất phản lực chân cột và gán ngược lại độ cứng lò xo phi tuyến theo đường cong $P-Y$ (cho cọc chịu lực ngang).

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

  • Sinh viên ngành Kỹ thuật Xây dựng / Địa kỹ thuật: Cung cấp tài liệu tham khảo chi tiết về phương pháp xử lý số liệu khảo sát địa chất thực tế và quy trình thực hiện đồ án tốt nghiệp chuẩn mực.
  • Kỹ sư Kết cấu & Địa kỹ thuật: Nắm bắt quy trình mô hình hóa móng băng trên nền đàn hồi trong SAP2000 và các bước kiểm tra cọc toàn diện (nâng cẩu, chịu tải trọng ngang, xuyên thủng).
  • Chủ đầu tư & Đơn vị Tư vấn Giám sát: Cơ sở để thẩm tra khối lượng, đánh giá mức độ an toàn của giải pháp móng và tối ưu hóa chi phí đầu tư xây dựng phần ngầm.

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

1. Yêu cầu đầu vào tối thiểu về địa chất để áp dụng quy trình tính toán này là gì?

Cần tối thiểu số liệu của 3–4 hố khoan khảo sát đạt chiều sâu ngàm vào tầng đất tốt ít nhất $5m$ (hoặc đạt chỉ số SPT $N > 30$), kèm bảng kết quả thí nghiệm cơ lý đầy đủ (thành phần hạt, giới hạn Atterberg, thí nghiệm cắt phẳng hoặc nén ba trục UU/CU).

2. Khi nào nên chuyển từ giải pháp móng băng sang móng cọc?

Khi chiều dày lớp đất yếu (bùn sét, sét dẻo chảy) lớn hơn $3.0\text{m}$, ứng suất tính toán của móng băng vượt quá sức chịu tải cho phép của đất nền ($P_{max} > R_{tc}$) hoặc độ lún tổng cộng tính toán $S > [S] = 8.0\text{ cm}$.

3. Tại sao phải kiểm tra cọc ở điều kiện cẩu lắp và dựng cọc?

Trong quá trình cẩu lắp từ bãi đúc lên xe và dựng cọc vào giá ép, cọc chịu sơ đồ tải trọng công xôn và uốn gãy do trọng lượng bản thân. Nếu cốt thép dọc không đủ chịu moment nâng cẩu, cọc sẽ bị nứt tế vi, gây rỉ cốt thép và giảm khả năng chịu lực khi ép sâu vào lòng đất.

4. SAP2000 tính toán dầm móng băng có xét đến biến dạng của đất nền không?

Có. Thông qua việc chia nhỏ dầm và gán lò xo đàn hồi tại các nút (Winkler Foundation), dầm móng biến dạng uốn tương thích với độ lún đàn hồi của các lò xo đất nền bên dưới.

5. Phương pháp hồi quy LINEST có ưu điểm gì so với việc vẽ đồ thị $\tau - \sigma$ bằng tay?

LINEST triệt tiêu hoàn toàn cảm tính chủ quan của người vẽ đường thẳng bao Coulomb; cung cấp chính xác giá trị hệ số góc ($\tan\varphi$), tung độ gốc ($c$), hệ số tương quan $R^2 > 0.95$ và sai số chuẩn của các hệ số theo thuật toán bình phương bé nhất.

Kết luận

Báo cáo đồ án môn học Nền Móng của sinh viên Nguyễn Thành Trí (Khoa Kỹ thuật Xây dựng, Trường Đại học Bách Khoa TP.HCM) là một công trình học thuật kết hợp thực hành thiết kế kỹ thuật hoàn chỉnh. Đồ án đã liên kết chặt chẽ chuỗi quy trình từ khảo sát địa chất thực địa $\rightarrow$ phân tích thống kê toán học $\rightarrow$ mô hình hóa kết cấu phần tử hữu hạn $\rightarrow$ kiểm tra giới hạn tiêu chuẩn quy phạm. Kết quả tính toán không chỉ đảm bảo độ an toàn chịu lực tuyệt đối cho công trình mà còn mang lại giải pháp kinh tế tối ưu, là tài liệu tham khảo giá trị cho cộng đồng kỹ sư xây dựng và sinh viên chuyên ngành địa kỹ thuật.