Tổng quan nghiên cứu

Nhu cầu xây dựng các công trình nhà cao tầng và siêu cao tầng từ 40 tầng đến trên 100 tầng, tiêu biểu như tòa tháp Taipei 101 cao 508m hay tháp đôi Petronas Towers cao 450m, đang gia tăng mạnh mẽ tại các đô thị hiện đại. Nếu chỉ sử dụng kết cấu bê tông cốt thép truyền thống cho một tòa nhà 40 tầng, kích thước tiết diện cột yêu cầu thường lên tới khoảng 1.5m x 1.5m, khiến công trình trở nên nặng nề, gia tăng tải trọng truyền xuống móng và làm giảm đáng kể không gian sử dụng hữu ích. Giải pháp kết cấu thép – bê tông cốt thép liên hợp giúp giảm kích thước tiết diện cột xuống còn khoảng 1.0m x 1.0m, tiết kiệm hơn 30% kích thước bao cảnh, giảm tải trọng bản thân và hạ giá thành xử lý nền móng từ 15% đến 25%.

Tuy nhiên, việc tính toán thiết kế khung liên hợp hiện nay còn nhiều bất cập. Các tiêu chuẩn quốc tế như Eurocode 4 hay AISC thường đơn giản hóa biến dạng trượt và phi tuyến vật liệu. Phương pháp khớp dẻo hay khớp dẻo hiệu chỉnh truyền thống chỉ giả định chảy dẻo cục bộ tại hai đầu phần tử, không phản ánh được sự lan truyền dẻo dọc theo chiều dài thanh và qua mặt cắt ngang. Đồng thời, ảnh hưởng của ứng suất dư chế tạo và phi tuyến hình học chuyển vị lớn chưa được đánh giá toàn diện.

Luận văn thạc sĩ chuyên ngành Xây dựng Dân dụng và Công nghiệp tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh (hoàn thành tháng 06/2012) tập trung giải quyết bài toán phân tích phi tuyến khung phẳng thép – bê tông liên hợp bằng phương pháp vùng dẻo. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi xây dựng mô hình số chuẩn xác, dự báo ứng xử phá hoại với sai số dưới 5% so với thực nghiệm, góp phần hoàn thiện cơ sở lý thuyết cho tiêu chuẩn xây dựng nhà cao tầng tại Việt Nam.

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 cơ học kết cấu phi tuyến, lý thuyết dầm Euler – Bernoulli mở rộng có xét trượt tiếp xúc và nguyên lý thế năng toàn phần dừng kết hợp phương pháp Rayleigh – Ritz. Các khái niệm và mô hình lý thuyết then chốt bao gồm:

Phương pháp vùng dẻo (Plastic Zone Method): Chia kết cấu thành nhiều phần tử hữu hạn và chia mặt cắt ngang thành nhiều thớ vi phân. Cách tiếp cận này cho phép theo dõi chính xác sự phát triển vùng dẻo, sự dịch chuyển trục trung hòa của lõi đàn hồi và ảnh hưởng của ứng suất dư trong thép hình.

Mô hình vật liệu phi tuyến: Vật liệu bê tông tuân theo mô hình Karayannis (1994) với biến dạng nén cực hạn đạt 0.002 và bỏ qua khả năng chịu kéo. Vật liệu thép hình và thép gia cường áp dụng mô hình đàn dẻo lý tưởng không tái bền.

Tương tác bán phần (Partial Interaction): Sử dụng mô hình phi tuyến của Ollgaard (1971) để mô tả quan hệ lực cắt – trượt tại bề mặt tiếp xúc giữa dầm thép và bản sàn bê tông, khắc phục giả định liên kết cứng tuyệt đối không thực tế.

Phi tuyến hình học: Sử dụng ten-xơ biến dạng Green trong hệ tọa độ Lagrangian nhằm mô tả chính xác hiệu ứng chuyển vị lớn và hiệu ứng P-Delta trong hệ khung chịu lực.

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

Nguồn dữ liệu và kiểm chứng: Luận văn sử dụng tập mẫu kiểm chứng gồm 7 bài toán thực nghiệm kinh điển quốc tế, bao gồm: dầm đơn giản E1 và U4 của Chapman & Balakrishnan (1964), dầm A3 của Haremza (2009), dầm liên tục hai nhịp CTB1 và CTB4 của Ansourian (1981), khung Portal của Liew & Chen (2001) và hệ khung 6 tầng 2 nhịp của Fang et al. (1999). Việc chọn mẫu thực nghiệm này đảm bảo tính đa dạng từ cấu kiện dầm đơn lẻ đến hệ kết cấu khung phức tạp.

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 vùng dẻo được lựa chọn vì đây là công cụ chính xác nhất hiện nay để mô tả trạng thái ứng xử không đàn hồi toàn diện, thay thế phương pháp khớp dẻo truyền thống vốn bỏ qua sự suy thoái độ cứng liên tục.

Quy trình thuật toán và công cụ: Tác giả phát triển hai phần tử hữu hạn mới gồm phần tử dầm – cột 6 bậc tự do (tương tác toàn phần) và phần tử dầm liên hợp 8 bậc tự do (tương tác bán phần). Thuật toán giải phi tuyến gia tăng Euler kết hợp kỹ thuật điều chỉnh công hằng được lập trình tự động hóa hoàn toàn bằng ngôn ngữ C++ qua hai chương trình SCCF và SCCB, thực hiện trong khung thời gian nghiên cứu 12 tháng (từ tháng 07/2011 đến tháng 06/2012).

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

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

Phần mềm C++ phát triển trong luận văn mô phỏng xuất sắc sự lan truyền dẻo: Chương trình SCCF và SCCB đã ghi nhận trực quan quá trình chuyển đổi từ trạng thái đàn hồi sang chảy dẻo từng phần của các thớ vi phân trên tiết diện dầm và cột, theo sát 100% bản chất vật lý của sự dịch chuyển trục trung hòa.

Ảnh hưởng rõ rệt của tương tác bán phần đến độ võng và độ cứng: Trên các mẫu dầm U4 và CTB4, sự trượt tại bề mặt tiếp xúc làm độ cứng uốn tổng thể của dầm giảm từ 10% đến 18%, đồng thời làm tăng độ võng ở giữa nhịp lên tới 25% so với trường hợp giả định tương tác toàn phần.

Tác động tiêu cực của ứng suất dư chế tạo: Đối với hệ khung liên hợp 6 tầng 2 nhịp, sự hiện diện của ứng suất dư trong thép hình làm giảm tải trọng tới hạn của hệ khung từ 6.0% đến 8.5%, đồng thời thúc đẩy quá trình chảy dẻo xuất hiện sớm hơn 15% so với cấu kiện không có ứng suất dư.

Xác định mật độ chia lưới tối ưu: Phân tích độ hội tụ chứng minh rằng việc chia 6 đến 8 phần tử trên mỗi cấu kiện thanh và chia từ 20 đến 40 thớ vi phân trên một mặt cắt ngang đảm bảo nghiệm số hội tụ ổn định với sai số tải trọng tới hạn nhỏ hơn 3% so với dữ liệu thực nghiệm của Ansourian (1981).

Thảo luận kết quả

Dữ liệu tính toán từ chương trình SCCF và SCCB được tổng hợp qua các bảng so sánh tải trọng tới hạn và trực quan hóa chi tiết bằng biểu đồ quan hệ tải trọng – chuyển vị đỉnh khung cũng như biểu đồ quan hệ lực cắt – trượt tại mặt phân giới thép và bê tông.

Khi so sánh với các nghiên cứu của Liew & Chen (2001) trên khung Portal và Fang et al. (1999) trên khung 6 tầng 2 nhịp, đường cong ứng xử phi tuyến do luận văn tính toán đạt mức độ tương đồng trên 95%. Nguyên nhân của sự suy giảm khả năng chịu lực trong giai đoạn sau đàn hồi được giải thích là do sự kết hợp đồng thời giữa hiện tượng nứt bê tông vùng kéo, chảy dẻo thép vùng nén và hiệu ứng suy giảm độ cứng hình học khi chuyển vị ngang tăng nhanh. Kết quả này khẳng định các phương pháp thiết kế truyền thống bỏ qua tương tác bán phần và ứng suất dư sẽ dẫn đến đánh giá thiếu an toàn về khả năng chịu lực cực hạn của công trình.

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

Áp dụng phần tử dầm liên hợp 8 bậc tự do trong tính toán độ võng thực tế: Các đơn vị tư vấn thiết kế công trình cần áp dụng mô hình phân tích có xét tương tác bán phần theo công thức Ollgaard nhằm kiểm soát độ trượt tiếp xúc dưới 2.0mm, bắt đầu áp dụng ngay trong giai đoạn thiết kế kỹ thuật từ năm 2026.

Bổ sung hệ số ảnh hưởng ứng suất dư vào quy chuẩn thiết kế: Đề nghị Viện Khoa học Công nghệ Xây dựng và các cơ quan quản lý chuyên ngành đưa quy định chiết giảm từ 5% đến 8% khả năng chịu lực của cột thép hình cán nóng khi phân tích phi tuyến vào tiêu chuẩn thiết kế nhà cao tầng giai đoạn 2026 – 2028.

Chuẩn hóa quy trình chia lưới phần tử hữu hạn vùng dẻo: Các kỹ sư kết cấu nên thiết lập mật độ chia lưới chuẩn gồm tối thiểu 6 phần tử dọc trục mỗi thanh và 30 thớ vi phân trên tiết diện ngang để vừa đảm bảo độ chính xác trên 97%, vừa khống chế thời gian chạy mô phỏng dưới 15 phút trên máy tính thông thường.

Mở rộng mô hình phân tích sang không gian ba chiều và tải trọng động: Các nhóm nghiên cứu chuyên sâu tại các trường đại học kỹ thuật cần tiếp tục nâng cấp mã nguồn C++ của chương trình SCCF để tích hợp bài toán khung không gian 3D chịu tải trọng động đất và cháy nổ trong lộ trình 24 tháng tới.

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

Học viên cao học và nghiên cứu sinh ngành Xây dựng: Tài liệu cung cấp cơ sở lý thuyết chuẩn mực về giải thuật phần tử hữu hạn phi tuyến, phương pháp năng lượng Rayleigh – Ritz và mã nguồn lập trình C++, phục vụ đắc lực cho các đề tài nghiên cứu chuyên sâu.

Kỹ sư thiết kế kết cấu công trình cao tầng: Giúp nắm bắt bản chất làm việc thực của hệ liên hợp thép – bê tông, từ đó tối ưu hóa kích thước tiết diện dầm cột, giúp tiết kiệm từ 10% đến 15% khối lượng vật liệu cho các tòa nhà thương mại trên 20 tầng.

Giảng viên và nhà nghiên cứu cơ học công trình: Sử dụng làm tài liệu tham khảo chất lượng cao trong giảng dạy các học phần Cơ học kết cấu nâng cao, Phân tích phi tuyến công trình và Thiết kế kết cấu liên hợp.

Chủ đầu tư và nhà quản lý dự án xây dựng: Có thêm căn cứ khoa học để phê duyệt giải pháp kết cấu liên hợp, rút ngắn tiến độ thi công khoảng 20% và tối đa hóa diện tích sàn thương phẩm hữu ích cho dự án.

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

Phương pháp vùng dẻo có ưu điểm gì vượt trội so với phương pháp khớp dẻo? Phương pháp khớp dẻo chỉ xem xét chảy dẻo cục bộ tại hai đầu mút và coi đoạn giữa hoàn toàn đàn hồi. Ngược lại, phương pháp vùng dẻo chia nhỏ tiết diện thành nhiều thớ vi phân, mô phỏng chính xác sự phát triển dẻo dọc thanh và qua tiết diện, giúp dự báo tải trọng giới hạn chính xác hơn 15%.

Tại sao bắt buộc phải xét đến hiện tượng tương tác bán phần trong dầm liên hợp? Trong thực tế, các chốt neo chống cắt luôn có biến dạng trượt đàn dẻo. Nếu bỏ qua hiện tượng này và coi liên kết cứng hoàn toàn, độ cứng uốn sẽ bị ước lượng cao hơn từ 10% đến 20%, dẫn đến tính thiếu độ võng thực tế của sàn từ 15% đến 25%, gây mất an toàn cho công trình.

Ứng suất dư ảnh hưởng như thế nào đến khả năng chịu lực của khung liên hợp? Ứng suất dư hình thành do quá trình nguội không đều khi cán thép hình. Khi chịu tải trọng lớn, ứng suất dư làm các thớ ngoài cùng bị chảy dẻo sớm hơn dự kiến, làm giảm tải trọng tới hạn của hệ khung liên hợp từ 6.0% đến 8.5% so với kết cấu lý tưởng.

Thuật toán điều chỉnh công hằng có vai trò gì trong giải bài toán phi tuyến? Khi kết cấu bước vào giai đoạn biến dạng dẻo lớn hoặc mất ổn định sau điểm cực hạn, ma trận độ cứng dễ bị suy biến. Thuật toán Euler kết hợp kỹ thuật điều chỉnh công hằng giúp tự động điều chỉnh bước tải, vượt qua các điểm giới hạn uốn cong và tránh hiện tượng phân kỳ nghiệm số.

Chương trình SCCF và SCCB viết bằng C++ có thể mở rộng phát triển như thế nào? Hai chương trình SCCF và SCCB được viết theo kiến trúc mô-đun hướng đối tượng rõ ràng cho phần tử 6 và 8 bậc tự do. Trong tương lai, mã nguồn hoàn toàn có thể tích hợp thêm ảnh hưởng của từ biến co ngót bê tông, tải trọng động đất theo lịch sử thời gian và mô hình liên kết nút nửa cứng.

Kết luận

  • Xây dựng thành công hai mô hình phần tử hữu hạn phi tuyến: phần tử dầm – cột 6 bậc tự do (tương tác toàn phần) và phần tử dầm liên hợp 8 bậc tự do (tương tác bán phần).
  • Phát triển toàn diện phương pháp vùng dẻo có xét đồng thời phi tuyến vật liệu, phi tuyến hình học Lagrangian và ứng suất dư chế tạo trong cấu kiện thép hình.
  • Lập trình thành công hai phần mềm chuyên dụng SCCF và SCCB bằng ngôn ngữ C++ với thuật toán gia tăng Euler và kỹ thuật điều chỉnh công hằng ổn định.
  • Thẩm tra độ chính xác thông qua 7 mô hình thực nghiệm quốc tế uy tín với độ tin cậy cao và sai số dự báo tải trọng tới hạn dưới 5%.
  • Khẳng định giải pháp kết cấu khung liên hợp giúp giảm trên 30% kích thước tiết diện cột so với bê tông cốt thép, mang lại hiệu quả kinh tế và thẩm mỹ vượt trội.

Đóng góp chính của luận văn là đã hoàn thiện công cụ tính toán số học thuật chính xác cao, lấp đầy khoảng trống nghiên cứu về phân tích vùng dẻo khung liên hợp tại Việt Nam.

Hướng nghiên cứu tiếp theo sẽ tập trung mở rộng thuật toán cho hệ khung không gian 3D chịu tải trọng động đất và phân tích ứng xử cháy trong lộ trình 12 đến 24 tháng tới.

Các kỹ sư và chuyên gia kết cấu quan tâm có thể khai thác phương pháp nghiên cứu này để tối ưu hóa thiết kế các công trình cao tầng hiện đại ngay hôm nay.