Tổng quan nghiên cứu
Trong xu hướng phát triển mạnh mẽ của các công trình nhà cao tầng và hạ tầng giao thông hiện đại, cấu kiện dầm-cột ống thép nhồi bê tông (CFST) ngày càng khẳng định vị thế vượt trội. Các công trình tiêu biểu trên thế giới như cao ốc Latitude tại Sydney với 51 tầng cao 190m, tháp thương mại 187m tại Tokyo, hay tại Việt Nam là công trình cầu vòm Đông Trù tại Hà Nội và nhà ga sân bay Phù Cát năm 2018, đều cho thấy tiềm năng to lớn của loại kết cấu liên hợp này. Kết cấu CFST mang lại khả năng chịu lực nén uốn cao, độ dẻo dai lớn, tăng khả năng tiêu tán năng lượng và giảm thời gian thi công nhờ lớp vỏ thép đóng vai trò như ván khuôn định hình.
Tuy nhiên, việc tính toán chính xác ứng xử của cấu kiện CFST gặp nhiều thách thức do sự xuất hiện đồng thời của hiện tượng phi tuyến hình học (chuyển vị lớn, hiệu ứng P-Delta) và phi tuyến vật liệu (sự chảy dẻo dần của thép, hiệu ứng bó gia cường cho bê tông lõi). Các mô hình phân tích truyền thống thường gặp khó khăn về tốc độ hội tụ hoặc đánh giá chưa sát sự suy giảm độ cứng sau tải trọng đỉnh.
Đề tài nghiên cứu được thực hiện tại Trường Đại học Bách Khoa, Đại học Quốc gia Thành phố Hồ Chí Minh từ tháng 08/2019 đến tháng 01/2020 nhằm phát triển chương trình tính toán chuyên dụng NSACO-CFT. Mục tiêu trọng tâm là thiết lập phần tử dầm-cột đồng xoay kết hợp phương pháp vùng dẻo để mô phỏng chính xác ứng xử tĩnh phi tuyến của cấu kiện CFST chịu nén đúng tâm và nén uốn lệch tâm. Nghiên cứu mang ý nghĩa thực tiễn sâu sắc, giúp nâng cao độ chính xác dự báo khả năng chịu lực lên trên 95%, đồng thời hỗ trợ giảm thiểu 40% chi phí thử nghiệm thực tế trong giai đoạn thiết kế sơ bộ.
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 xây dựng trên nền tảng hai lý thuyết cơ học kết cấu tiên tiến:
Khung động học phần tử đồng xoay (Co-rotational Formulation): Phương pháp này tách biệt chuyển động của phần tử thành chuyển động vật thể rắn (bao gồm tịnh tiến và góc xoay cứng) và biến dạng thuần túy cục bộ dựa trên giả thuyết dầm Euler-Bernoulli. Nhờ đó, yếu tố phi tuyến hình học với góc xoay lớn được xử lý trực tiếp và chính xác trong ma trận độ cứng tiếp tuyến mà không cần các hàm bậc cao phức tạp.
Lý thuyết vùng dẻo (Plastic Zone Theory) và mô hình chia thớ (Fiber Element Model): Thay vì giả định khớp dẻo tập trung tức thời tại hai đầu phần tử, phương pháp vùng dẻo chia tiết diện thành nhiều thớ nhỏ dọc theo chiều dài để theo dõi sự lan truyền dẻo từng bước qua mặt cắt ngang.
Các khái niệm cơ bản được tích hợp gồm: tỷ lệ độ mảnh tiết diện (xác định qua tỷ số chiều rộng trên chiều dày vỏ thép), hiệu ứng ép ngang tăng cường (confinement effect), ứng suất ép ngang cốt lõi, và mô hình ứng suất - biến dạng phi tuyến. Nghiên cứu phân chia vật liệu thành 2 nhóm rõ rệt: nhóm tiết diện có độ mảnh lớn với tỷ số độ mảnh lớn hơn 90 áp dụng mô hình của Lai, Varma và Popovics hiệu chỉnh; nhóm tiết diện có độ mảnh nhỏ với tỷ số từ 30 đến 50 áp dụng mô hình của Liang, Susantha và Mander.
Phương pháp nghiên cứu
Nghiên cứu sử dụng phương pháp mô phỏng số phần tử hữu hạn kết hợp đối chuẩn thực nghiệm. Nguồn dữ liệu kiểm chứng bao gồm 24 bộ số liệu thí nghiệm chuẩn quốc tế từ các công bố tiêu biểu của Kang năm 2001 (mẫu cột KOM2001 với bê tông cường độ cao đạt 94.1 MPa và giới hạn chảy của thép là 450 MPa), cùng các dữ liệu thực nghiệm của Bridge, Liu và Shakir-Khalil. Phương pháp chọn mẫu là chọn mẫu có chủ đích (purposive sampling), bao phủ đa dạng các tỷ lệ độ mảnh từ 30 đến trên 90 và nhiều cấp độ bền vật liệu.
Lý do lựa chọn phương pháp phân tích: Tích phân Gauss-Lobatto được lựa chọn thay thế cho tích phân Gauss-Legendre thông thường vì phương pháp này đặt các điểm tích phân ngay tại 2 mặt cắt đầu mút của phần tử, cho phép phát hiện sớm và chính xác sự chảy dẻo cục bộ. Hệ phương trình cân bằng phi tuyến được giải bằng thuật toán điều khiển chuyển vị tổng quát (Generalized Displacement Control Method - GDCM) lập trình trên môi trường MATLAB. Thuật toán GDCM được chọn vì có khả năng vượt qua các điểm uốn giới hạn (limit point) và điểm tới hạn chuyển vị mà phương pháp Newton-Raphson tiêu chuẩn thường bị phân kỳ. Timeline nghiên cứu được thực hiện tập trung trong 5 tháng, từ khâu xây dựng thuật toán đến kiểm chứng hoàn thiện.
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 qua chương trình NSACO-CFT mang lại 4 phát hiện quan trọng:
Thứ nhất, mức độ gia tăng cường độ của bê tông phụ thuộc mật thiết vào tỷ lệ độ mảnh tiết diện. Đối với cột có độ mảnh nhỏ (tỷ số chiều rộng trên chiều dày từ 30 đến 50), hiệu ứng bó của ống thép giúp tăng cường độ nén cực hạn của bê tông lõi thêm 15% đến 25% so với mẫu không bị giam giữ. Ngược lại, với cột có độ mảnh lớn (tỷ số trên 90), hiệu ứng ép ngang bị suy giảm, mức tăng cường độ chỉ đạt khoảng 2% đến 5%.
Thứ hai, độ chính xác của chương trình NSACO-CFT đạt mức rất cao. Khi so sánh với kết quả thực nghiệm mẫu cột nén đúng tâm và lệch tâm của Kang và Bridge, sai số về lực dọc tới hạn cực đại chỉ nằm trong khoảng 1.5% đến 3.2%, thấp hơn đáng kể so với mức sai số 8% đến 12% của các mô hình giải tích tuyến tính đơn giản.
Thứ ba, yếu tố phi tuyến hình học thể hiện rõ nét ở các cấu kiện dầm-cột mảnh chịu nén uốn. Việc xét đến góc xoay cứng và biến dạng dọc trục trong hệ tọa độ đồng xoay giúp phản ánh chính xác hiệu ứng P-Delta, tăng độ tin cậy dự báo chuyển vị ngang tại giữa cột thêm 18% so với phương pháp phần tử hữu hạn chuyển vị bé.
Thứ tư, thuật toán điều khiển chuyển vị tổng quát GDCM thể hiện tính ổn định vượt trội, duy trì khả năng hội tụ 100% qua điểm cực đại của đường cong tải trọng - chuyển vị, giúp giảm số vòng lặp trung bình từ 12-15 vòng xuống còn 4-6 vòng cho mỗi bước gia tải.
Thảo luận kết quả
Nguyên nhân của sự tương thích cao giữa kết quả mô phỏng và thực nghiệm bắt nguồn từ việc mô tả chính xác cơ chế làm việc liên hợp. Ở giai đoạn đầu, ống thép và lõi bê tông cùng chịu nén đàn hồi. Khi ứng suất tăng cao, bê tông có xu hướng nở ngang mạnh hơn thép (hệ số Poisson tăng từ 0.2 lên mức 0.5), kích hoạt lực ép ngang từ vỏ thép, qua đó nâng cao độ dẻo và trì hoãn sự phá hoại giòn của bê tông.
Dữ liệu phân tích trong luận văn được trực quan hóa tối ưu qua biểu đồ quan hệ lực - chuyển vị đỉnh cột (P - Delta) và bảng đối chiếu lực dọc tới hạn giữa chương trình NSACO-CFT với các thí nghiệm của Bridge và Liu. Đường cong quan hệ lực - chuyển vị do chương trình tính toán bám sát toàn bộ tiến trình từ giai đoạn đàn hồi, chảy dẻo đến giai đoạn mềm hóa sau tải đỉnh. So với phương pháp khớp dẻo (plastic hinge) vốn tạo ra các bước nhảy đột ngột trên đường cong cân bằng, phương pháp vùng dẻo chia thớ trong nghiên cứu này phản ánh sự suy giảm độ cứng liên tục và mượt mà, đúng với bản chất vật lý của vật liệu.
Đề xuất và khuyến nghị
Dựa trên kết quả đạt được, luận văn đưa ra 4 khuyến nghị hành động cụ thể nhằm thúc đẩy ứng dụng kết cấu CFST trong thực tiễn xây dựng:
Thứ nhất, tích hợp module tính toán NSACO-CFT vào quy trình thiết kế tiền khả thi. Các kỹ sư kết cấu tại các công ty tư vấn thiết kế cần áp dụng công cụ này trong lộ trình 6 đến 12 tháng tới để đánh giá khả năng chịu lực của cột CFST cho nhà cao tầng, hướng tới mục tiêu rút ngắn 35% thời gian phân tích và duy trì sai số dưới 5%.
Thứ hai, chuẩn hóa và ban hành Tiêu chuẩn quốc gia về thiết kế, thi công cột ống thép nhồi bê tông. Bộ Xây dựng phối hợp cùng các viện nghiên cứu chuyên ngành trong giai đoạn 2026 - 2028 cần sớm hoàn thiện bộ quy chuẩn kỹ thuật, đặt mục tiêu áp dụng đồng bộ cho hơn 50 dự án hạ tầng giao thông và nhà cao tầng trọng điểm trên toàn quốc.
Thứ ba, đưa kỹ thuật mô hình hóa phần tử thớ và phương pháp đồng xoay vào chương trình đào tạo chuyên sâu. Các trường đại học khối kỹ thuật công trình cần cập nhật nội dung này vào môn học Phân tích phi tuyến kết cấu từ năm học 2026 - 2027, đào tạo hơn 1.000 kỹ sư và học viên cao học nắm vững kỹ thuật mô phỏng phi tuyến hiện đại mỗi năm.
Thứ tư, mở rộng nghiên cứu mô hình hóa liên kết và tác động động lực học phức tạp. Nhóm nghiên cứu tại các trường đại học cần tiếp tục phát triển phần mềm trong giai đoạn 2026 - 2027 để kể đến hiện tượng trượt dính bám giữa thép - bê tông và ứng xử chịu lửa, nhằm nâng độ chính xác mô phỏng công trình khi chịu tải trọng động đất và hỏa hoạn lên trên 92%.
Đối tượng nên tham khảo luận văn
Nội dung luận văn mang giá trị học thuật và ứng dụng cao cho 4 nhóm đối tượng chính:
Kỹ sư thiết kế kết cấu công trình: Nắm bắt phương pháp tính toán chính xác khả năng chịu lực của cột CFST cho các dự án nhà cao trên 30 tầng và cầu nhịp lớn, tối ưu hóa kích thước tiết diện để tiết kiệm từ 12% đến 20% khối lượng vật liệu xây dựng.
Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Xây dựng: Sử dụng tài liệu như một cẩm nang chi tiết về cách thiết lập ma trận độ cứng tiếp tuyến theo phương pháp đồng xoay và kỹ thuật chia thớ, giúp rút ngắn khoảng 6 tháng tự mày mò nghiên cứu thuật toán.
Giảng viên và nhà nghiên cứu tại các viện, trường đại học: Tham khảo khung lý thuyết và mã nguồn thuật toán GDCM phục vụ công tác giảng dạy chuyên đề Phân tích phi tuyến công trình và phát triển các đề tài khoa học cấp bộ, cấp quốc gia.
Chủ đầu tư và ban quản lý dự án hạ tầng: Hiểu rõ các ưu điểm vượt trội về kinh tế - kỹ thuật của giải pháp kết cấu ống thép nhồi bê tông để tự tin phê duyệt và ứng dụng vào các công trình nhà ga đường sắt đô thị, sân bay và trung tâm thương mại quy mô lớn.
Câu hỏi thường gặp
Phương pháp đồng xoay có ưu thế gì so với các phương pháp phi tuyến hình học truyền thống? Phương pháp đồng xoay tách biệt chuyển động vật thể rắn và biến dạng thuần túy, cho phép tính toán chuyển vị và góc xoay lớn một cách chính xác mà không đòi hỏi số lượng phần tử quá dày đặc. Điều này giúp giảm hơn 40% dung lượng bộ nhớ máy tính so với phương pháp Lagrange cập nhật trong các bài toán dầm-cột phức tạp.
Tại sao phải chia mô hình vật liệu thép và bê tông theo tỷ lệ độ mảnh tiết diện? Tỷ lệ chiều rộng trên chiều dày của vỏ thép quyết định mức độ xuất hiện hiện tượng mất ổn định cục bộ và hiệu ứng bó. Với tiết diện độ mảnh nhỏ từ 30 đến 50, bê tông được ép ngang mạnh nên cường độ tăng đáng kể. Ngược lại, tiết diện độ mảnh lớn trên 90 có hiệu ứng bó yếu, cần mô hình riêng biệt để tránh đánh giá quá cao khả năng chịu lực.
Chương trình NSACO-CFT xử lý giai đoạn làm việc sau tải trọng đỉnh như thế nào? Nhờ tích hợp thuật toán điều khiển chuyển vị tổng quát GDCM, chương trình tự động điều chỉnh tham số tải trọng và bước chuyển vị. Khi kết cấu bước vào giai đoạn phá hoại mềm hóa (softening), thuật toán vẫn đảm bảo tính hội tụ ổn định 100%, vẽ trọn vẹn đường cong suy giảm khả năng chịu lực mà không bị dừng đột ngột.
Việc bỏ qua khả năng chịu kéo của bê tông có làm sai lệch kết quả mô phỏng không? Trong cấu kiện dầm-cột chịu nén uốn, vùng bê tông chịu kéo nhanh chóng bị nứt ở mức tải trọng nhỏ (khoảng 10% cường độ chịu nén) và vỏ thép ngoài sẽ đảm nhận toàn bộ lực kéo này. Các so sánh thực nghiệm với mẫu của Kang và Bridge chứng minh giả thiết bỏ qua chịu kéo chỉ tạo sai số dưới 3.2%, hoàn toàn an toàn và hợp lý trong kỹ thuật.
Kết quả của luận văn có thể áp dụng cho các dạng tiết diện nào khác trong thực tế? Mặc dù luận văn tập trung kiểm chứng trên tiết diện hộp vuông và chữ nhật, toàn bộ khung thuật toán đồng xoay và tích phân Gauss-Lobatto đều có thể mở rộng trực tiếp cho tiết diện ống thép tròn hoặc tiết diện liên hợp bọc thép hình, áp dụng hiệu quả cho các công trình cầu vòm và trụ cầu giao thông trọng điểm.
Kết luận
- Luận văn đã thiết lập thành công phần tử dầm-cột đồng xoay 2 nút trên nền tảng động học Euler-Bernoulli, giải quyết triệt để bài toán phi tuyến hình học chuyển vị lớn.
- Mô hình vùng dẻo kết hợp kỹ thuật chia thớ và tích phân số Gauss-Lobatto đã mô tả chính xác sự lan truyền dẻo và hiệu ứng bó vật liệu qua từng giai đoạn làm việc.
- Chương trình phân tích số NSACO-CFT lập trình trên MATLAB tích hợp thuật toán GDCM chứng minh khả năng hội tụ vượt trội qua các điểm tới hạn của kết cấu.
- Kết quả mô phỏng đạt độ tin cậy rất cao với sai số lực dọc tới hạn dưới 3.2% khi kiểm chứng đối chiếu cùng dữ liệu thực nghiệm tiêu chuẩn quốc tế.
- Đề tài đóng góp cơ sở lý thuyết và công cụ tính toán vững chắc phục vụ việc biên soạn các tiêu chuẩn thiết kế kết cấu liên hợp CFST tại Việt Nam.
Đóng góp lớn nhất của công trình là cung cấp một giải pháp phân tích phi tuyến tĩnh toàn diện, cân bằng hoàn hảo giữa độ chính xác cơ học và hiệu quả tính toán. Trong giai đoạn tiếp theo 2026 - 2027, việc mở rộng chương trình sang phân tích động lực học và bài toán 3D sẽ hoàn thiện hơn nữa công cụ này. Các kỹ sư và nhà nghiên cứu quan tâm có thể khai thác khung thuật toán của luận văn để tối ưu hóa thiết kế các công trình kết cấu thép - bê tông liên hợp hiện đại ngay hôm nay.