Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình giao thông, kết cấu móng mố trụ cầu chiếm khoảng 25% đến 35% tổng mức đầu tư xây lắp của một dự án cầu vượt sông quy mô lớn. Bệ cọc đóng vai trò là bộ phận trung gian chịu lực đặc biệt phức tạp, có nhiệm vụ tiếp nhận toàn bộ tải trọng thẳng đứng, lực cắt và mômen uốn lên tới hàng chục nghìn kilonewton từ thân trụ để phân phối đều xuống hệ thống cọc ma sát và nền đất. Hiện nay, thực tiễn thiết kế công trình cầu tại Việt Nam vẫn chủ yếu dựa trên hai phương pháp cốt lõi: phương pháp dầm truyền thống căn cứ theo giả thiết Bernoulli và phương pháp mô hình chống giằng (Strut-and-Tie Model - STM). Mặc dù mô hình dàn ảo phản ánh tốt cơ chế truyền lực ở các vùng không liên tục, quy trình tính toán thử sai phức tạp khiến phương pháp này khó ứng dụng nhanh trong giai đoạn thiết kế sơ bộ. Ngược lại, phương pháp dầm giản đơn có ưu điểm tính toán nhanh nhưng mức độ tin cậy và sai số ứng suất chưa được kiểm chứng định lượng cụ thể.

Mục tiêu trọng tâm của nghiên cứu là khảo sát, đánh giá một cách toàn diện mức độ phù hợp của các sơ đồ tính đơn giản dựa trên lý thuyết dầm khi phân tích bệ cọc, lấy kết quả phân tích phần tử hữu hạn 3D làm hệ quy chiếu chuẩn xác. Phạm vi nghiên cứu tập trung vào kết cấu bệ cọc móng trụ T8 thuộc công trình cầu Sông Lũy thi công bằng phương pháp đúc hẫng cân bằng, thực hiện khảo sát biến thiên thông số trên 9 trường hợp hình học với chiều cao bệ từ 700 mm đến 1400 mm và khoảng cách cọc từ 2.5D đến 3.0D. Kết quả nghiên cứu có ý nghĩa khoa học và thực tiễn cao, giúp rút ngắn khoảng 40% thời gian tính toán sơ bộ cho kỹ sư thiết kế, đồng thời cung cấp hệ số hiệu chỉnh ứng suất chính xác nhằm kiểm soát an toàn kết cấu móng cầu đường.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng kết hợp các lý thuyết cơ học kết cấu kinh điển và tiêu chuẩn thiết kế hiện đại:

  • Lý thuyết dầm Euler-Bernoulli và nguyên lý Saint-Venant: Giả thiết tiết diện phẳng được áp dụng để xác định nội lực uốn và cắt trong dầm thông thường. Tuy nhiên, nguyên lý Saint-Venant chỉ ra rằng trong phạm vi từ 1.0 đến 2.0 lần kích thước hình học tính từ điểm đặt tải trọng hoặc gối tựa, trường ứng suất bị xáo trộn mạnh mẽ tạo thành vùng không liên tục (vùng D), nơi lý thuyết dầm thông thường không còn duy trì độ chính xác cao.
  • Mô hình dàn ảo (Strut-and-Tie Model) theo tiêu chuẩn ACI 318-08 và 22 TCN 272-05: Cấu trúc bệ cọc được lý tưởng hóa thành một hệ dàn không gian gồm các thanh chống bê tông chịu nén (Strut), thanh giằng cốt thép chịu kéo (Tie) và các vùng nút liên kết (Nodal zones dạng C-C-C, C-C-T). Cường độ chịu nén hữu hiệu của thanh chống được kiểm soát qua hệ số chiết giảm vật liệu bê tông.
  • Thuyết bền thế năng biến đổi hình dáng cực đại Huber - von Mises: Được sử dụng để đánh giá trạng thái ứng suất tương đương phức hợp trong không gian 3 chiều của khối bê tông bệ cọc, so sánh trực tiếp với cường độ chịu nén tính toán của vật liệu.
  • Lý thuyết tương tác nhóm cọc và nền đất: Áp dụng công thức Converse - Labarre và tiêu chuẩn BOCA 1993 quy định khoảng cách tối thiểu giữa các tim cọc ma sát $s \ge 2.5D$ (tối thiểu 760 mm) nhằm triệt tiêu hiện tượng chồng chéo vùng nén lún và suy giảm sức chịu tải của nhóm cọc.

Phương pháp nghiên cứu

  • Cỡ mẫu và trường hợp khảo sát: Nghiên cứu thiết lập hệ thống gồm 9 mô hình tính toán biến thiên toàn diện các thông số hình học: thay đổi 3 cấp chiều cao bệ cọc (H = 2.0 m, 2.5 m, 3.0 m) và 3 cấp khoảng cách tim cọc (s = 2.5D, 2.75D, 3.0D với đường kính cọc D = 1.0 m).
  • Phương pháp chọn mẫu: Sử dụng phương pháp chọn mẫu điển hình theo chủ đích (purposive sampling), lấy nguyên mẫu trụ T8 cầu Sông Lũy – một kết cấu cầu dầm hộp bê tông cốt thép dự ứng lực đúc hẫng cân bằng thực tế có cấu tạo địa chất phức tạp gồm 4 lớp đất nền đặc trưng (từ bùn sét yếu đến cát thô lẫn sỏi chặt vừa).
  • Phương pháp phân tích và lý do lựa chọn: Ứng dụng phương pháp số phần tử hữu hạn (FEM) trên phần mềm ANSYS kết hợp phần mềm chuyên dụng CDC PRO 1.0. Tác giả lựa chọn loại phần tử khối bậc hai SOLID186 (20 nút, 3 bậc tự do mỗi nút) và SOLID187 (10 nút dạng tứ diện) cho mô hình bệ cọc và đoạn ngàm cọc. Lý do lựa chọn là vì phần tử khối 3D bậc cao cho phép mô phỏng chính xác hiện tượng tập trung ứng suất cục bộ, tính phi tuyến hình học và loại bỏ hiện tượng khóa trượt cắt (shear locking) mà các phần tử thanh 1D hay tấm 2D không thể thực hiện được. Chiều dài tự do chịu uốn của cọc trong nền đất được xác định chuẩn xác qua công thức ngàm tương đương $L_M = L_0 + 7d$. Quá trình tính toán, mô phỏng và phân tích đối chứng được thực hiện liên tục trong thời gian 6 tháng (từ tháng 01/2015 đến tháng 06/2015).

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

  • Phát hiện 1 (Mô hình dầm console có độ tương thích cao nhất): Trong các sơ đồ dầm giản đơn được khảo sát, sơ đồ tính toán bệ cọc theo dạng dầm công-xôn (consol) ngàm tại mép chân thân trụ cho thấy sự tương quan phân bố ứng suất sát nhất với kết quả mô hình khối 3D. Độ lệch chuẩn về vị trí xuất hiện mômen uốn cực đại giữa hai phương pháp chỉ dao động dưới 15% tại các vùng tiết diện xa gối tựa.
  • Phát hiện 2 (Hiện tượng tập trung ứng suất và sai số ứng suất pháp): Tại vị trí góc liên kết giữa chân thân trụ và mặt trên của bệ cọc xuất hiện hiện tượng tập trung ứng suất cực kỳ lớn do hiệu ứng góc nhọn và sự thay đổi độ cứng đột ngột. Ứng suất pháp kéo uốn $\sigma_y$ tính theo phương pháp dầm thông thường thấp hơn từ 2.8 lần đến 3.2 lần so với kết quả giải tích chi tiết bằng phần tử khối 3D ANSYS. Từ dữ liệu thực nghiệm mô phỏng, tác giả đã đề xuất hệ số điều chỉnh chuẩn hóa là $\alpha = 3.0$ cho công thức tính ứng suất theo phương pháp dầm.
  • Phát hiện 3 (Quy luật độ cứng tương đối): Nghiên cứu xác lập quy luật tương quan nghịch giữa "độ cứng tương đối" ($\lambda$) của bệ cọc và trị số ứng suất cục bộ tại các vị trí nguy hiểm. Khi tăng chiều cao bệ cọc giúp độ cứng tương đối $\lambda$ tăng thêm 25%, mức độ tập trung ứng suất tại đáy bệ cọc giảm từ 18% đến 22%, đồng thời chuyển vị lệch giữa các đầu cọc biên giảm rõ rệt.
  • Phát hiện 4 (Đảm bảo điều kiện an toàn chịu nén tổng thể): Trên toàn bộ 9 kịch bản phân tích ở cả hai giai đoạn thi công đúc hẫng và khai thác toàn tải, giá trị ứng suất nén lớn nhất trong bệ cọc dao động từ 8.5 MPa đến 14.2 MPa, hoàn toàn nhỏ hơn giới hạn cường độ chịu nén tính toán của bê tông cấp độ bền B30 ($R_b = 17.5\text{ MPa}$).

Thảo luận kết quả

Sự chênh lệch đáng kể giữa mô hình dầm và mô hình phần tử khối 3D bắt nguồn từ bản chất cơ học của vùng D (Discontinuity regions). Phương pháp dầm Bernoulli thuần túy giả định ứng suất biến thiên tuyến tính trên toàn bộ chiều cao tiết diện. Tuy nhiên, đối với kết cấu khối lớn như bệ cọc móng cầu, sự truyền lực diễn ra theo cơ chế vòm không gian ba chiều với các luồng ứng suất nén dạng hình nón tỏa từ đáy trụ xuống từng đầu cọc.

Trong báo cáo kỹ thuật, dữ liệu so sánh được trực quan hóa sinh động thông qua biểu đồ phân bố đường đẳng ứng suất Von Mises và bảng đối chiếu sai số tương đối tại 7 điểm kiểm soát trên 2 mặt cắt nguy hiểm (Mặt cắt 1-1 tại tim trụ và Mặt cắt 2-2 tại mép thân trụ). Khi khoảng cách tim cọc tăng từ 2.5D lên 3.0D, mômen uốn cục bộ tại mặt cắt chân trụ tăng khoảng 14.8%, đòi hỏi diện tích cốt thép gia cường đáy đài phải tăng tương ứng để khống chế bề rộng vết nứt. Phát hiện này hoàn toàn phù hợp và phát triển sâu hơn các nhận định của các nghiên cứu quốc tế trước đây về tương tác kết cấu - nền đất, khẳng định tính tất yếu của việc đưa hệ số hiệu chỉnh vào các mô hình tính toán đơn giản.

Đề xuất và khuyến nghị

  • Chuẩn hóa hệ số điều chỉnh trong thiết kế sơ bộ: Đề xuất các đơn vị tư vấn thiết kế cầu đường áp dụng hệ số hiệu chỉnh $\alpha = 3.0$ khi sử dụng sơ đồ dầm consol để tính toán ứng suất pháp uốn và bố trí cốt thép chịu lực cho bệ cọc. Giải pháp này giúp nâng cao độ an toàn kết cấu đạt hệ số an toàn $K \ge 1.5$ ngay từ bước thiết kế cơ sở mà không làm tăng khối lượng vật liệu quá mức. Chủ thể thực hiện: Kỹ sư kết cấu mố trụ cầu. Thời gian triển khai: Áp dụng ngay trong quy trình tính toán sơ bộ (1 - 2 tuần).
  • Tối ưu hóa kích thước hình học đài cọc: Khuyến nghị lựa chọn chiều cao bệ cọc tối thiểu theo quan hệ kinh tế $H \ge 600\text{ mm}$ (áp dụng công thức Reynolds $H = 2h_p + 100\text{ mm}$) và duy trì khoảng cách cọc tối ưu trong khoảng từ 2.5D đến 2.75D. Việc kiểm soát khoảng cách này giúp triệt tiêu hoàn toàn sự chồng chéo ứng suất trong đất nền, đồng thời giúp tiết kiệm từ 12% đến 15% khối lượng bê tông bệ móng so với phương án bố trí khoảng cách 3.0D trở lên. Chủ thể thực hiện: Chủ nhiệm đồ án và kỹ sư thiết kế chính.
  • Quy chuẩn hóa mô hình chia lưới phần tử hữu hạn: Khi thực hiện kiểm toán chi tiết bằng phần mềm ANSYS, các kỹ sư cần thiết lập mật độ lưới phần tử SOLID186 mịn tại khu vực liên kết chân trụ với kích thước phần tử từ 100 mm đến 150 mm để đảm bảo sai số hội tụ ứng suất nhỏ hơn 5%. Chủ thể thực hiện: Bộ phận mô phỏng số và tính toán kết cấu chuyên sâu (CAE).
  • Xây dựng quy trình tự động hóa liên kết phần mềm: Kiến nghị các viện nghiên cứu và doanh nghiệp tư vấn giao thông phối hợp phát triển module tự động trích xuất phản lực đầu cọc từ phần mềm CDC PRO 1.0 sang mô hình khối ANSYS nhằm rút ngắn 50% thời gian phân tích kết cấu móng phức tạp. Timeline: Triển khai thử nghiệm và chuyển giao công nghệ trong vòng 6 tháng.

Đối tượng nên tham khảo luận văn

  • Kỹ sư thiết kế kết cấu cầu - đường bộ: Nắm vững phương pháp ứng dụng hệ số điều chỉnh $\alpha = 3.0$ và sơ đồ dầm công-xôn để tính toán nhanh, chính xác nội lực bệ cọc cho các công trình cầu nhịp lớn chịu tải trọng trên 5.000 kN.
  • Chuyên viên thẩm tra và quản lý dự án hạ tầng giao thông: Sử dụng các kết quả nghiên cứu làm căn cứ kỹ thuật để đánh giá, phản biện tính hợp lý của hồ sơ thiết kế móng cọc từ 6 đến 9 cọc, kiểm soát rủi ro nứt vỡ cục bộ đài cọc trong quá trình thi công đúc hẫng.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Công trình Giao thông / Địa kỹ thuật: Tiếp cận phương pháp luận nghiên cứu kết hợp giữa mô phỏng số phần tử hữu hạn 3D nâng cao, thuyết bền Von Mises và các tiêu chuẩn quốc tế như ACI 318-08, 22 TCN 272-05.
  • Giảng viên và cơ sở đào tạo kỹ sư xây dựng: Bổ sung tư liệu giảng dạy thực tế cho học phần Thiết kế Móng Cầu, Mô phỏng số trong Xây dựng và so sánh chuyên sâu giữa lý thuyết dầm Bernoulli với mô hình chống giằng Strut-and-Tie.

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

  • Tại sao phương pháp dầm truyền thống lại đánh giá thấp ứng suất tại bệ cọc?
    Phương pháp dầm truyền thống dựa trên giả thiết tiết diện phẳng của Bernoulli, vốn chỉ chính xác cho các cấu kiện uốn thanh mảnh (vùng B). Bệ cọc là kết cấu khối lớn có tỷ lệ chiều cao trên nhịp ngắn, chịu tác dụng tải trọng tập trung cao độ tạo ra vùng không liên tục (vùng D). Tại đây, ứng suất thực tế tập trung rất lớn ở chân trụ, cao gấp khoảng 3 lần so với tính toán tuyến tính thông thường.

  • Hệ số điều chỉnh $\alpha = 3.0$ được áp dụng trong phạm vi và điều kiện nào?
    Hệ số $\alpha = 3.0$ được áp dụng để nhân với giá trị ứng suất pháp kéo uốn tính toán theo sơ đồ dầm giản đơn khi thiết kế sơ bộ bệ cọc bê tông cốt thép đối xứng (từ 6 đến 9 cọc ma sát) có khoảng cách cọc từ 2.5D đến 3.0D, đảm bảo an toàn chịu lực trước khi kiểm toán chi tiết bằng mô hình phần tử hữu hạn 3D.

  • Khi nào bắt buộc phải sử dụng mô hình dàn ảo (Strut-and-Tie Model) thay cho mô hình dầm?
    Theo tiêu chuẩn ACI 318-08 và 22 TCN 272-05, mô hình dàn ảo bắt buộc phải áp dụng khi bệ cọc có tỷ lệ khoảng cách truyền lực trên chiều sâu hữu dụng $a_v/d < 2.0$ (kết cấu dầm cao/khối lớn) hoặc khi cấu tạo bệ cọc có hình học phi đối xứng, vị trí đặt cọc phức tạp không thể phân tách thành các dải dầm phẳng độc lập.

  • Việc sử dụng phần tử SOLID186 trong ANSYS mang lại ưu thế kỹ thuật gì?
    Phần tử SOLID186 là phần tử khối bậc hai đẳng tham số gồm 20 nút với 3 bậc tự do tịnh tiến trên mỗi nút. Phần tử này có khả năng mô phỏng chính xác trường biến dạng phi tuyến, hiện tượng tập trung ứng suất không gian 3 chiều và vật liệu đàn dẻo, loại bỏ hoàn toàn sai số do hiệu ứng khóa trượt cắt thường gặp ở phần tử bậc một 8 nút.

  • Khoảng cách bố trí giữa các tim cọc tối ưu được xác định như thế nào?
    Theo các tiêu chuẩn xây dựng quốc tế như BOCA và NBC, khoảng cách tối thiểu giữa các tim cọc ma sát phải đạt $s \ge 2.5D$ (không nhỏ hơn 760 mm). Khoảng cách tối ưu từ 2.5D đến 2.75D vừa đảm bảo loại trừ hiện tượng giao thoa ứng suất gây phá hoại cắt đất nền giữa các cọc, vừa tối ưu hóa kích thước đài móng và tiết kiệm chi phí vật liệu.

Kết luận

  • Luận văn đã phân tích toàn diện giới hạn của phương pháp dầm Bernoulli truyền thống và khẳng định tính tất yếu của việc sử dụng mô hình phần tử hữu hạn 3D trong phân tích ứng suất bệ cọc.
  • Xác lập sơ đồ dầm công-xôn (consol) là mô hình tính toán đơn giản có mức độ tương thích cao nhất phục vụ công tác thiết kế sơ bộ kết cấu móng mố trụ cầu.
  • Đề xuất thành công hệ số hiệu chỉnh chuẩn hóa $\alpha = 3.0$ cho bài toán tính ứng suất pháp uốn theo mô hình dầm giản đơn, khắc phục hiện tượng đánh giá thấp ứng suất cục bộ tại chân trụ.
  • Chứng minh quy luật tương quan nghịch giữa độ cứng tương đối $\lambda$ của bệ cọc và mức độ tập trung ứng suất tại các điểm xung yếu đáy đài móng.
  • Khẳng định 100% các trường hợp khảo sát móng trụ T8 cầu Sông Lũy đều đảm bảo an toàn tuyệt đối về khả năng chịu nén của vật liệu bê tông ($< 17.5\text{ MPa}$).

Công trình nghiên cứu đã đóng góp cơ sở lý luận và thực nghiệm số quan trọng, tạo tiền đề thúc đẩy chuyển đổi số và ứng dụng công nghệ mô phỏng CAE trong ngành xây dựng cầu đường Việt Nam. Trong định hướng giai đoạn 2026-2027, các nghiên cứu tiếp theo nên mở rộng phân tích tương tác phi tuyến đất - cọc dưới tác động của tải trọng động động đất và sóng thần. Hãy liên hệ ngay với các chuyên gia kỹ thuật kết cấu để nhận tài liệu hướng dẫn và tối ưu hóa giải pháp thiết kế bệ cọc cho công trình của bạn!