Tổng quan nghiên cứu

Trong bối cảnh đô thị hóa nhanh chóng, nhu cầu xây dựng các tòa nhà cao tầng từ 40 tầng trở lên tại các đô thị lớn như Hà Nội và Thành phố Hồ Chí Minh đang gia tăng mạnh mẽ. Khi áp dụng các giải pháp kết cấu bê tông cốt thép truyền thống, tiết diện dầm cột thường phải tăng từ 25% đến 35%, khiến kết cấu trở nên nặng nề, tốn kém chi phí vật liệu và làm giảm đáng kể diện tích sử dụng hữu ích của công trình. Để khắc phục triệt để nhược điểm này, giải pháp dầm bê tông cốt cứng (Fully Encased Beam - FEB) – một dạng kết cấu liên hợp thép và bê tông đặc biệt với lõi thép hình được bao bọc hoàn toàn bên trong bê tông – đã trở thành hướng đi mang tính đột phá.

Đề tài luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng được thực hiện trong thời gian 06 tháng tại Trường Đại học Kinh Tế Công Nghiệp Long An dưới sự hướng dẫn khoa học của PGS.TS Trương Tích Thiện. Mục tiêu cốt lõi của nghiên cứu là xây dựng quy trình phân tích số mô phỏng chính xác ứng xử uốn, khả năng chịu lực giới hạn và cơ chế rạn nứt của dầm bê tông cốt cứng tiết diện 250 x 400 mm có cốt cứng chữ I bằng chương trình phần tử hữu hạn ANSYS.

Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu từ 60% đến 70% chi phí chế tạo mẫu và thử nghiệm phá hủy trong phòng thí nghiệm. Đồng thời, kết quả nghiên cứu cung cấp cơ sở dữ liệu tin cậy giúp các kỹ sư kết cấu rút ngắn 40% thời gian thiết kế, kiểm soát tối ưu độ võng và độ nứt công trình theo các tiêu chuẩn quốc tế hiện hành.

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 lý thuyết biến dạng phi tuyến của vật liệu liên hợp thép - bê tông, kết hợp các mô hình phá hoại và ứng xử đàn - dẻo tiên tiến:

  • Mô hình phá hủy bê tông Willam - Warnke (1974): Sử dụng bề mặt phá hủy 5 thông số để mô tả trạng thái ứng suất đa trục phức tạp của bê tông, kết hợp tiêu chuẩn chảy dẻo Von-Mises để xác định chính xác ứng xử nứt khi chịu kéo (ft = 1.5 MPa) và vỡ vụn khi chịu nén (fc' = 19.002 MPa).
  • Mô hình ứng suất - biến dạng Hognestad cải biên (1951) và Mander (1988): Mô tả quan hệ ứng suất - biến dạng phi tuyến của bê tông chịu nén một trục với nhánh parabol tăng dần đến biến dạng đỉnh xấp xỉ 0.002 và mô đun đàn hồi ban đầu Ec = 20875 MPa theo tiêu chuẩn ACI 318-08.
  • Mô hình đàn - dẻo song tuyến tính cho thép: Ứng dụng cho cả cốt thép dọc chịu lực (giới hạn chảy fys = 343.8 MPa, mô đun đàn hồi Es = 2.1 x 10^5 MPa) và thép hình chữ I (giới hạn chảy fyc = 280 MPa), cho phép mô tả chính xác giai đoạn biến cứng sau chảy dẻo.
  • Các khái niệm then chốt: Dầm liên hợp bọc bê tông toàn phần (FEB), trạng thái ứng suất thủy tĩnh, hệ số truyền lực cắt qua khe nứt (hệ số mở và hệ số đóng), và tương tác đồng chuyển vị giữa thép hình và bê tông.

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

Nghiên cứu kết hợp chặt chẽ giữa tổng quan lý thuyết, mô phỏng số học nâng cao và kiểm chứng thực nghiệm đối chứng:

  • Nguồn dữ liệu và cỡ mẫu: Dữ liệu hình học và cơ lý được thu thập từ mô hình thực nghiệm chuẩn dầm A1 trong nghiên cứu của Huỳnh Phúc Linh và Hồ Hữu Chỉnh (2011). Tiết diện dầm khảo sát gồm chiều rộng B = 250 mm, chiều cao D = 400 mm, cốt thép dọc gồm 4 thanh đường kính 14 mm, và lõi thép hình chữ I có diện tích Ays = 3440 mm2 (chiều cao ds = 200 mm, bề rộng cánh bf = 146 mm, chiều dày cánh tf = 8 mm, chiều dày bụng tw = 6 mm).
  • Phương pháp chọn mẫu và xây dựng mô hình: Lựa chọn phương pháp mô hình hóa phân tách rời rạc trong ANSYS. Bê tông được mô hình hóa bằng phần tử khối 3D SOLID65 (8 nút, 3 bậc tự do mỗi nút), cốt thép thanh dọc mô hình bằng phần tử LINK180 (2 nút chịu kéo/nén), và thép hình chữ I được mô hình bằng phần tử vỏ phẳng SHELL181 (4 nút, 6 bậc tự do mỗi nút).
  • Lý do lựa chọn phương pháp phân tích: Phương pháp phần tử hữu hạn trên ANSYS cho phép phân tích phi tuyến hình học và phi tuyến vật liệu ở mức độ chi tiết cao, đồng thời kỹ thuật liên kết nút (Couple DOF và Merge Nodes) giúp mô phỏng chính xác sự tương tác tiếp xúc hoàn toàn giữa cốt cứng và khối bê tông mà giải tích cổ điển không thể giải quyết được.
  • Timeline nghiên cứu: Triển khai tuần tự trong 06 tháng, bao gồm 02 tháng tổng quan lý thuyết và xác lập thông số vật liệu, 02 tháng lập trình code APDL và phân tích số trong ANSYS, 02 tháng so sánh đánh giá sai số với thực nghiệm và mở rộng bài toán tham số.

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

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

Quá trình mô phỏng và kiểm chứng dầm bê tông cốt cứng trên phần mềm ANSYS đã đem lại các phát hiện khoa học quan trọng:

  1. Độ võng giữa dầm phản ánh trung thực qua các cấp tải: Tại mức tải trọng 5000 N, độ võng đàn hồi giữa dầm ghi nhận được là rất nhỏ; khi tăng tải lên 10000 N và 13300 N, biến dạng uốn tăng dần theo quy luật phi tuyến rõ nét. Độ võng tính toán từ ANSYS bám sát đường cong thực nghiệm với độ sai lệch cực đại chỉ trong khoảng 5% đến 8%, khẳng định độ tin cậy của mô hình số.
  2. Quy luật phát triển vết nứt bê tông: Vết nứt kéo đầu tiên xuất hiện ở thớ đáy giữa dầm khi tải trọng đạt khoảng 30% đến 40% tải trọng giới hạn (ứng suất kéo vượt ngưỡng 1.5 MPa). Khi tải trọng tiến sát mức 13300 N, mật độ vết nứt mở rộng nhanh chóng và lan dần lên trục trung hòa, phù hợp với cơ chế phá hoại dẻo do uốn trong thực tế.
  3. Phân bố trường ứng suất Von-Mises: Trường ứng suất tương đương trong cốt thép dọc và thép hình chữ I phân bố tập trung cao nhất tại vùng giữa nhịp. Ứng suất trong cánh dưới của thép hình chữ I đạt ngưỡng 280 MPa và cốt thép kéo đạt 343.8 MPa tại tải trọng 13300 N, chứng minh lõi thép hình đã gánh từ 50% đến 65% tổng mômen uốn sau khi bê tông vùng kéo bị nứt hoàn toàn.
  4. Hiệu ứng gia cường của kích thước thép hình: Khi khảo sát 3 trường hợp thay đổi kích thước thép hình chữ I, việc tăng bề rộng cánh bf từ 100 mm lên 146 mm giúp mômen quán tính tương đương tăng lên đáng kể, làm giảm độ võng giữa dầm từ 18% đến 22% tại cùng một cấp tải trọng uốn.

Thảo luận kết quả

Sự tương đồng chặt chẽ giữa mô phỏng ANSYS và dữ liệu thực nghiệm bắt nguồn từ việc thiết lập chính xác các thông số truyền lực cắt. Khi đặt hệ số truyền lực cắt cho vết nứt mở ở mức 0.2 đến 0.3 và vết nứt đóng ở mức 1.0, bài toán phi tuyến đạt độ hội tụ ổn định tuyệt đối mà không xảy ra hiện tượng mất cân bằng số học. Đồng thời, việc gán giá trị -1 để tắt tính năng nén vỡ đơn trục trong SOLID65 theo khuyến nghị của các nhà nghiên cứu quốc tế đã loại bỏ hoàn toàn các lỗi phân kỳ toán học khi bê tông bước vào giai đoạn biến mềm.

Khi so sánh với các nghiên cứu của Mazen AL-Bdoor (2013) về cột cốt cứng và Eva Maria Regy (2017) về dầm liên hợp chịu nhiệt, kết quả của luận văn này củng cố thêm khẳng định: kết cấu bê tông cốt cứng có khả năng duy trì độ cứng chống uốn sau nứt cao hơn 30% so với dầm bê tông cốt thép thông thường có cùng lượng tiêu hao vật liệu. Toàn bộ quan hệ tải trọng - chuyển vị và sự chuyển dịch trục trung hòa có thể được trực quan hóa sinh động qua biểu đồ đường cong uốn P - delta và bảng phân tích trường ứng suất phân lớp theo chiều cao dầm.

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

Dựa trên các kết quả phân tích số và kiểm chứng thực nghiệm, luận văn đưa ra 4 giải pháp và khuyến nghị kỹ thuật có tính ứng dụng cao:

  • Ứng dụng quy trình mô phỏng số ANSYS vào thiết kế sơ bộ: Các công ty tư vấn thiết kế xây dựng cần đưa quy trình mô phỏng phần tử hữu hạn sử dụng phần tử SOLID65, LINK180 và SHELL181 vào giai đoạn thiết kế kỹ thuật cho các dự án cao tầng trên 20 tầng, mục tiêu rút ngắn 35% thời gian thẩm tra kết cấu và hoàn thành chuyển giao công nghệ trong giai đoạn 2026-2027.
  • Chuẩn hóa thiết lập thông số vật liệu nứt trong phần mềm: Kỹ sư kết cấu khi mô hình bê tông mác 300 cần cố định hệ số truyền lực cắt vết nứt mở trong khoảng 0.2 đến 0.3 và hệ số vết nứt đóng bằng 1.0 nhằm đảm bảo tỷ lệ hội tụ đạt 100% trong phân tích phi tuyến tải trọng tăng dần.
  • Tối ưu hóa tỷ lệ diện tích lõi thép hình trong dầm: Các chủ đầu tư và nhà thầu xây dựng nên áp dụng tỷ lệ diện tích thép hình chữ I trên diện tích tiết diện bê tông trong khoảng 3% đến 6% để tăng khả năng chịu tải trọng uốn thêm 25% đến 30%, áp dụng ngay cho các dự án khởi công trong vòng 12 tháng tới.
  • Biên soạn hướng dẫn kỹ thuật liên kết nút liên hợp: Bộ Xây dựng và các viện nghiên cứu chuyên ngành cần chủ trì xây dựng tài liệu hướng dẫn kỹ thuật chuẩn hóa phương pháp Couple DOF và Merge Nodes trong mô hình liên hợp thép - bê tông, dự kiến ban hành bổ sung vào tiêu chuẩn thiết kế kết cấu liên hợp quốc gia trước quý IV/2027.

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

Nội dung và dữ liệu thực nghiệm trong luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  • Kỹ sư kết cấu và chuyên viên R&D tại các công ty tư vấn thiết kế: Nắm vững phương pháp khai báo vật liệu phi tuyến, cách chia lưới phần tử khối và kỹ thuật xử lý tiếp xúc giữa bê tông và cốt thép hình chữ I trong ANSYS, giúp tối ưu hóa kích thước dầm cho các dự án nhịp lớn.
  • Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Công trình Xây dựng: Sử dụng toàn bộ cấu trúc code APDL trong phần phụ lục và khung lý thuyết Willam - Warnke làm tài liệu tham khảo chuẩn để phát triển các đề tài mô phỏng kết cấu phức tạp.
  • Giảng viên và nhà nghiên cứu tại các trường đại học kỹ thuật: Sử dụng làm tài liệu giảng dạy chuyên đề Phương pháp phần tử hữu hạn nâng cao, ứng dụng ANSYS trong cơ học kết cấu và phân tích công trình chịu tải trọng tĩnh/động.
  • Ban quản lý dự án và các nhà phát triển bất động sản: Có thêm căn cứ kỹ thuật để đánh giá tính khả thi kinh tế - kỹ thuật của giải pháp dầm cốt cứng, giúp tiết kiệm 15% chi phí kết cấu và tăng diện tích thông thủy cho các tòa nhà văn phòng, căn hộ cao cấp.

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

Dầm bê tông cốt cứng FEB khác gì so với dầm bê tông cốt thép thông thường?
Dầm bê tông cốt cứng đặt toàn bộ thép hình chữ I bên trong bê tông thay vì chỉ dùng các thanh cốt thép tròn. Cấu tạo này giúp tăng khả năng chịu lực uốn lên hơn 40%, giảm độ võng giữa nhịp từ 18% đến 22% và tăng khả năng kháng cháy vượt trội so với dầm thép hở ngoài.

Tại sao phần tử SOLID65 lại phù hợp nhất để mô hình bê tông trong ANSYS?
SOLID65 là phần tử khối 3D 8 nút chuyên dụng, tích hợp sẵn mô hình Willam - Warnke để mô phỏng chính xác sự hình thành vết nứt kéo theo 3 phương trực giao, sự biến cứng, biến mềm và khả năng chịu nén đa trục của bê tông dưới các cấp tải trọng từ 5000 N đến 13300 N.

Kỹ thuật Couple DOF và Merge Nodes trong ANSYS có vai trò gì?
Do phần tử bê tông SOLID65 có 3 bậc tự do tịnh tiến mỗi nút, trong khi phần tử thép hình SHELL181 có 6 bậc tự do mỗi nút, kỹ thuật Couple DOF hoặc Merge Nodes giúp ràng buộc chặt chẽ các bậc tự do chuyển vị trùng nhau, đảm bảo sự truyền lực dính bám hoàn hảo giữa thép và bê tông.

Vì sao phải tắt tính năng nén vỡ đơn trục bằng cách gán giá trị -1 trong ANSYS?
Khi kích hoạt tính năng nén vỡ đơn trục, phần mềm thường xuyên gặp lỗi phân kỳ toán học khi vật liệu bước vào miền phi tuyến cao. Gán giá trị -1 giúp chương trình duy trì sự hội tụ ổn định của thuật toán Newton - Raphson mà vẫn bảo toàn độ chính xác về độ võng và phân bố ứng suất.

Độ chính xác của mô hình ANSYS trong luận văn được kiểm chứng qua đâu?
Mô hình số được kiểm chứng trực tiếp với dữ liệu dầm thực nghiệm A1 chịu uốn 4 điểm. Kết quả chuyển vị độ võng tại các mức tải trọng 5000 N, 10000 N và 13300 N cho thấy sai số so với thực nghiệm chỉ dao động từ 5% đến 8%, chứng minh độ tin cậy tuyệt đối của phương pháp.

Kết luận

  • Luận văn đã xây dựng thành công quy trình mô phỏng phần tử hữu hạn hoàn chỉnh cho dầm bê tông cốt cứng FEB trên phần mềm ANSYS với độ chính xác cao, sai số so với thực nghiệm dưới 8%.
  • Xác lập bộ thông số chuẩn cho mô hình vật liệu bê tông mác 300 và thép hình chữ I, giải quyết triệt để bài toán hội tụ số bằng cách tinh chỉnh hệ số truyền lực cắt vết nứt mở ở mức 0.2 đến 0.3.
  • Chứng minh rõ nét hiệu quả của việc tăng kích thước cánh thép hình chữ I trong việc giảm 18% đến 22% độ võng uốn và tăng khả năng chịu mômen uốn của cấu kiện.
  • Đóng góp nguồn dữ liệu tham khảo quý giá cùng toàn bộ mã code APDL phục vụ công tác nghiên cứu, giảng dạy và ứng dụng thực tiễn trong ngành xây dựng Việt Nam.
  • Trong lộ trình 12 tháng tiếp theo, nghiên cứu cần được mở rộng cho bài toán dầm liên hợp chịu tải trọng động, tải trọng lặp tuần hoàn và tác động của nhiệt độ cao khi hỏa hoạn.

Hãy ứng dụng ngay giải pháp dầm bê tông cốt cứng và phương pháp phân tích số trên ANSYS để tối ưu hóa an toàn chịu lực và nâng cao hiệu quả kinh tế cho các dự án công trình cao tầng của bạn.